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XXXVI. A New Method of Determining the Charge and Energy of Heavy 
Nuclei in the Cosmic Radiation. 


By A. D. Datyton, P. H. Fowxer and D. W. Kenr, 
H. H. Wills Physical Laboratory, University of Bristol*, 


[Received February 5, 1951.] 


SUMMARY. 


The charge and energy of a heavy particle in cosmic radiation may be 
determined from the frequency of the accompanying delta rays and the 
multiple scattering in its path through an emulsion. The method is 
illustrated by the charge spectrum of 226 tracks in emulsions exposed at 
95,000 feet. 


§1. INTRODUCTION. 


Since the discovery of the presence of heavy nuclei in the cosmic radiation 
(Freier et. al. 1948), several experimenters have employed the photographic 
method to determine the magnitude of the charges of the individual 
particles. The problem is of great importance because a knowledge of the 
“ charge spectrum ” of the primary radiation is essential for any discussion 
of its origin (Bradt and Peters 1950 a and b). 

The method employed hitherto has been based on the measurement of, 
(a), the frequency of occurrence, N;, per unit length of the track, of 
S-rays with energies in excess of a certain minimum value ; together with, 
(b), the ranges of the particles in the assembly of glass and emulsion 
constituting the photographic plates employed to record the tracks. At 
geomagnetic latitudes less than 30°, most of the primary particles are 

‘moving at relativistic velocities. For such particles, N; varies as Z* and 
this relation is independent of the precise value of the energy. The charge 
can therefore be deduced from the observed value of N, without an exact 
knowledge of the range, except in the few cases where an observable 
change of N, with range indicates that a particle is of lower energy. 

Bradt and Peters (1950 b) have recently refined this method by using 
insensitive emulsions. The grain-density in the tracks of multiply 
charged particles is then so reduced that it can be determined by grain 
counts. Such observations give a measure of the ionization of the particle 
with a much greater statistical weight than a determination of N; deduced 
from an equal length of track. 

For particles of lower energy, which occur frequently in the primary 
radiation at geomagnetic latitudes greater than 45°, the above methods 
suffer from the fact that, in general, the charge and the energy can be 
accurately determined only if the particle happens to reach the end of its | 
range in one of the emulsions of the assembly. It is therefore important to 
develop other modes of attack, and for this purpose we have found it 
possible to modify the methods employed in this laboratory in determining 
the energy of fast protons and other particles of charge |e]. 


* Communicated by Prof. C. F. Powell, F.R.S. 
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In the observation on the tracks of fast protons, the energy of a particle 
can be determined, and its identity established in favourable cases, by 
measuring, (a), the grain density, g, in the track ; and, (6), the average 
angular deviation in the direction of motion, a, per unit length of the track, 
due to multiple coulomb scattering. In extending this method to the 
heavy particles, we have made similar observations on the scattering, but 
instead of the grain-density in the track, we determine N, the frequency 
of occurrence of 5-rays. By these observations it is possible to determine 
the charge and energy of a particle if its track in a single emulsion is long 
enough to allow measurement of « and N, to be made with sufficient 
statistical weight. A knowledge of the range of the particle is unnecessary. 
Further, the method appears to be sufficiently precise to allow elements of 
low atomic number to be identified with little ambiguity, so that the charge 
spectrum shows clearly resolved peaks due to the different elements. 


Fig ls 


50 MICRONS 


GRAINS per 


fe) 2 “4 “6 8 ye) 
DELTA RAYS per 100 MICRONS 
Relation between grain-density, g, and delta-ray density, N3, in the trajectory 
of particles of charge |e]. Measurements on 7-mesons are shown thus: @ ; 
on protons and deuterons, thus: ©; and on fast «-particles, thus: ©. 


§2. RELATION BETWEEN N, AND VELOCITY FOR PARTICLES OF 
DIFFERENT CHARGE. 


In order to determine the relationship between the 6-ray density in the 
track of a particle of known charge and its velocity, we have made 
observations on particles of charge |e]. For these particles, the relation- 
ship between residual range and velocity, and between grain-density and 
velocity, is now well established, for the type of emulsion in which the 
present tracks are recorded. 
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In the first instance, we made measurements on the tracks of a+ or > 
particles of which the range in the emulsion was about 8 mm., and which 
were identified by the secondary effects observed at the end of their range ; 
viz., either -decay, or the nuclear absorption of the particle and the 
formation of a ‘‘ star’ with two or more branches. The observed values 
of N, were plotted against the corresponding values of grain-density in the 
track, for various values of the residual ran ge, and the results are shown in 
fig. 1. In making the observations, only 8-rays giving rise to a track 
containing four or more grains were accepted. 

Similar observations of greater statistical weight were made on the tracks 
of fifty protons and deuterons emitted from nuclear explosions occurring in 
the emulsion, particles of which the identity had been established by the 


DELTA RAYS per 100 MICRONS 


RANGE in ems, 


i i t bon and of a 
iati f delta-ray density with range for tracks of a car 
“Eee Sos shat reach the end of their range in the emulsion, 

The nuclei were identified by observations of Ng and range. 


methods previously described (Fowler 1950). Tt will be seen from ae 1 
that for particles of charge |e|, moving at velocities for which the Bea 
density is less than five times the minimum value (g <59min)s the nes ot 
N, are proportional to the grain-density. On the basis of these obsery a- 
tions, curves can be drawn showing the relation between Me and ae 
scattering parameter «—for particles with any value of the ghar Ze a » d 
mass, A—using the relation between « and g established for z-particles 


(Camerini et al. 1950). 
Z2 


320 A. D. Dainton et al. on a New Method of Determining the 


As a check on equation (1), values of N, were determined by observations 
on long tracks of magnesium and carbon nuclei which came to rest in the 
emulsion, and which could be identified by the original method based on 
determinations of N,and range. The observations showed that 


Nockee uo) ae 


where R is the residual range; see fig. 2. Using this result and the 
established relation between range and energy for particles moving in a 
photographic emulsion, 

Bice byes, (2) 


the relation between N, and 8 was computed*. The result thus obtained 
was found to be consistent, within the limits of experimental error, with 
the calculated family of curves. 


§3. OBSERVATIONS ON THE SCATTERING OF Heavy NUCLEI. 


The measurements of the deviations in the tracks of heavy nuclei due to 
multiple coulomb scattering present no particular difficulties and have 
been carried out in the manner described by Fowler (1950). The thickness 
of the track gives rise to larger reading errors than in the case of similar 
observations with fast protons or «-particles, but these are not important. 
In order to evaluate the mean angle of scattering from the readings taken 
on the microscope, we have corrected for reading errors, stage “ noise ”’ and 
distortion, by the procedures to be described in a later communication. 
If A is the mass number and Z the charge number of the nucleus, its energy, 
corresponding to a given value of «, can be calculated from the relation :— 

a= Zh Apv, sa. own be vie 
where p is the momentum of the nucleus per nucleon. With the units 
commonly employed, the scattering constant k is taken as 26, a value 
indicated by a preliminary determination made in this laboratory. With 
the present exposures, “noise”? and distortion prevent any accurate 
measurement of energies greater than 3 BeV./nucleon. It may, however, 
be possible to refine the processing technique in order to diminish the 
distortion errors and thus to improve the precision of the measurements. 


§4. CHARGE SPECTRUM OF Heavy NUCLEI. 


The present methods were employed in a study of the charge distribution 
of the heavy nuclei recorded in an assembly of Ilford G5 plates exposed by 
means of free balloons at an altitude of 95,000 ft., 40 c.c. of the emulsion of 
these plates were searched for tracks of which the lengthin any one emulsion 
was greater than 6 mm., and in which the value of N; was greater than 
1:0 per 1004. For each of these tracks the values of « and N, were 
determined. ‘The results thus obtained are represented in fig. 3. 


aa naar n UNI connunsetCSSLIBLS ssl eee ie es ee 

* Equation (3) is not seriously invalidated, for particles with values of the 
energy with which we are here dealing, by the processes of ‘ capture and loss ” 
of electrons occurring at the end of the range of the heavy particles. 


ENERGY . per NUCLEON 
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The family of full lines in fig. 3 represent the expected relation between 
N; and « for different values of the charge number Z. They have been 
calculated on the assumption that the mass number of a nucleus, A, is’ 
equal to 2Z. This assumption is approximately justified for stable nuclei 


whilst the dashed lines are 


e figure. 


le plate (see fig. 2), are shown 


h 


lines of constant range. An energy scale is included at the top of t 


Observations on tracks found in a 
Measurements on two tracks 


systematic survey of 40c.c. of emulsion are shown thus :—O. 


per IOO MICRONS 


x 


(100) for the tracks of 230 particles. 


a 


which terminate in the emulsion after travelling over 4 cm. in a sing 
thus :— +--+. The full lines represent lines of constant charge, 


Relation between N; and 


Ng NUMBER of 5-RAYS per 1OO MICRONS. 


up to Z~20; any nuclei of greater charge seldom occur in the primary 
radiation. Further, it seems very unlikely that any but stable or very 
long-lived isotopes will have survived the long periods believed to he 
involved in the acceleration process which produces the primary particles. 


322 A. D. Dainton et al. on a New Method of Determining the 


The approximation is therefore adequate, with the present state of precision 
of the mass measurements, in the case of the primary particles. The 
greatest errors are likely to occur in the case of nuclei of lithium and 
beryllium, but they can hardly be serious. 

From the results, the value of the charge for each individual particle 
can be determined and the resulting charge “ spectrum” is shown in 
fig. 4.. The spectrum includes results due to slow «-particles, many of 
which are products of the fragmentation of heavier nuclei as a result of 
nuclear collisions. Lithium nuclei are also present with energies as low 
as 60 MeV./nucleon. Energetic lithium nuclei are sometimes emitted 
from nuclear explosions caused by protons, neutrons, or «-particles 


a ari 


le ie ee 


Fig. 4. 
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Number of Particles. 


Charge Number Z. 
Charge spectrum of 226 of the tracks represented in fig. 3. The four other tracks 
observed in the systematic survey were of charge Z greater than 12; three of 
them were identified as due to silicon, and the fourth due to acalcium nucleus. 


(Sérensen 1951), but they occur very infrequently. Most of the particles 
of charge Z=3 observed in the present experiment must therefore have 
been either produced by the fragmentation of heavier nuclei, or they are 
true primary particles. 

The charge spectrum shown in fig. 4 has not been corrected for loss of : 
particles by ionization—an effect which would tend to reduce the numbers 
of more massive nuclei observed—or for the general degradation in charge 
due to nuclear collisions in the atmosphere above the plates. Even after 
making allowance for these factors, however, it is difficult to account for 
the observed numbers of beryllium and boron nuclei if these elements 
are absent from the primary radiation at this latitude. Such a result, if 


— 
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it can be manitained by observations of greater statistical weight, would 
appear to be at variance with those of other observers working at a 
geomagnetic latitude of 30° N. (Bradt and Peters 1950 b). 


§5. APPRAISAL OF THE METHOD. 


The most important advantage of the present method of determining 
the charge of fast atomic nuclei over those previously employed is that, at 
high geomagnetic latitudes, if the track of a particle is of sufficient length, 
the observations on it can be confined to a single emulsion. Further, for 
values of Z between three and six inclusive, the accuracy of the measure- 
ments appears to be greater than that previously attained (Bradt 
and Peters 1950 a, Hoang 1950, Dainton and Kent 1950). It is difficult, 
for example, to evaluate the charge of such particles, even when they are of 
low velocity, from the changes in the values of N, as the particle passes 
from one emulsion to another in traversing the stack. 

In addition to the determination of the charge, the method gives the 
energy of the particles as a result of the same limited observations and a 
scale of energy per nucleon is included in fig. 3. It would be premature 
to discuss the detailed features of the observations at this stage of the 
investigation. 
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XXXVI. An Interferometric Study of Cleavage Surfaces of Artificially 
Grown Orystals of Sodium Chloride. 


By 8. AMELINCEX, 
University of Ghent, Belgium*. 


[Received January 30, 1951.] 
[Plate XIV.] 


ABSTRACT. 


Cleavage surfaces of an artificially grown crystal of sodium chloride are 
examined by multiple-beam interferometry. Two main features are 
noted, namely long narrow, more or less parallel strips, and superposed 
on this is a secondary fine grain structure. Silver films on sodium chloride 
are found to be unstable and apparently photosensitive. The cleavage 
mechanism is discussed. 


CLEAVAGE faces of an artificially grown crystal of NaCl have been examined 
by means of optical interference methods, using both multiple-beam 
Fizeau fringes and fringes of equal chromatic order (Tolansky 1948). 
The silvered specimen was matched against a high-grade silvered optical 
flat accurate to within A/40. 

Cleavage was induced in a large clear crystal (grown by Hilger) by - 
striking perpendicular to a face with a razor blade. Rocksalt shows three 
cleavage directions parallel to (001), (010) and (100), which are considered 
as crystallographically equivalent. This was confirmed since no specific 
difference could be detected between them. 

Although the cleavage is generally described as “* perfect ’’, interfero- 
metric examination revealed a rather complex structure on these faces. 
The Fizeau fringes (A 5461) for a typical low magnification of a well-cleaved 
specimen is shown on fig. la. This shows the presence of narrow strips, 
generally more or less parallel to each other and radiating out from the 
point where the cleavage was induced. There is a close resemblance 
between this picture and that of the calcite cleavage plane reported by 
Tolansky and Khamsavi (1946). 

The strip widths vary from 1/100 to 1/10 mm. and their lengths are of 
the order of | mm. The planes of the strips are not always parallel but 
are sometimes inclined one to another as in the case of selenite (Tolansky 
1945). In fig. la the strip AB is inclined to the strip CD at an angle of 
1-33 x 10-4 radians. ; 

Another fact of interest is clearly visible. The part of the cleavage face 
in the right bottom corner (below the line A’B’) is inclined to the adjacent 


* Communicated by Professor Tolansky. 
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region above A’B’ at an angle of 2-11<10-* radians. The part of the 
crystal face to the left of the strip CD consists of strips all very nearl 
parallel to the strip CD; whereas the part to the right of AB menace 
of strips all very nearly parallel to AB. 

The strips CD and AB themselves are inclined one to another. The 
part of the crystal face shown here consists consequently of tinge main 
regions which are inclined one to another over very small angles. 


Fig. 2 


N 


Distance along 
Surface Eres 


os 0:5 ma. Y '0mm. 


Fig. 1 (Pl. XIV.) shows some other characteristics of the rocksalt cleavage 
structure rarely seen on other cleavage structures investigated such as 
mica (Tolansky 1945, Tolansky and Morris 1947), selenite and calcite 
(loc. cit.). Besides well defined steps the change in level sometimes takes 
place gradually as can be seen along the lines ED and FG on fig. la. 
Fig. 2 represents a cross-section along a line XY, overlapping order 
allocations being ascertained by using several wavelengths of light. The 
gradual changes in level are probably due to the non-resolution of very 
small steps. The elementary distance, perpendicular to the face as given 
by X-ray diffraction, is dy9).=5°63 A. which is beyond the resolving power 


of the method. 


326 S. Amelinckx on an Interferometric Study of Cleavage 


Although in general within a narrow strip cleavage must probably 
proceed along a lattice plane, in some cases curved strips were seen (fig. 1 5). 
Some radii of curvature calculated by using the approximation : 
R=1/2A{d,d,(d, +d,)|/[d, —d.], where d, and d, are the distances between 
three adjacent fringes, are for example, R=1-58, 2-02 and 3-29 metres 
respectively. 

Uniform areas are rare and very small and in general high dispersion 
gives ragged fringes, indicating a complicated fine structure as in 
figs. le and 1d. This fine structure, probably also due to non-resolved 
steps, corresponds to level changes up to 200 A. 

Fig. 16 shows another striking feature of the rocksalt cleavage. The 
surface is striated in two directions: one parallel to (100) and the other 
parallel to the cleavage lines. The striations lead to the dotting of the 
fringes. 

The following proposed cleavage mechanism is suggested by the fringe 
patterns. When a knock is applied to a crystal face, cleavage starts in 
the plane under the razor blade’s edge containing the largest number of 
dislocations or defects. Cleavage will proceed along this plane until 
reaching the immediate neighbourhood of some other plane which contains 
a higher concentration of dislocations. Cleavage then Jumps and continues 
along that new plane. The result is a repeated change of level and the 
production of steps. This mechanism would probably account for the 
smaller steps only. The higher steps could be associated with inner 
cracks existing previously in the crystal or with the boundaries of the 
separate blocks which build up a coarse mosaic structure of the crystal. 
Evidence in favour of this last argument may be found in the fact that 
the large step of fig. 2 occurs in the immediate neighbourhood of the 
border between two inclined crystal-blocks. | 

According to this view the quality of a cleavage would depend on the 
concentration of lattice defects and on their distribution in the crystal 
and consequently on growth conditions. This view is supported by 
fig. 2 of Tolansky and Khamsavi's paper on calcite (1946). 

This investigation of the rocksalt cleavage involved certain difficulties. 
As rocksalt is hygroscopic it was necessary to silver the specimens 
immediately after cleavage. Keeping the silvered specimens in a 
desiccator filled with P,O, was not sufficient to prevent the immediate 
attack on the rocksalt. One hour only is needed to make this effect 
visible. It results in a decrease of reflectivity and an increase of the 
absorption of the silver layer. A specimen is completely spoiled and 
unusable for examination some twelve hours after silvering. 

The break-up of the film is much more visible with reflection fringes 
than with transmission, as the local loss of reflectivity is then seen as a 
dark region on a bright background. In transmission.on the contrary 
the effect is seen only as an interruption in the bright fringes. |The 
disintegration could therefore be studied much easier in reflection. There 
is some indication that it starts in regions of maximum illumination 
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which, as Holden (1949) has shown is asymmetrically on one side of the 
dark reflection fringe. It appears to spread out in the direction of the 
illuminated region and stops abruptly when a cleavage line is reached 
(fig. le). 

Although the silver break-up is undoubtedly partly due to the hygro- 
scopic character of rocksalt, it may be that a photosensitive effect 
contributes if it be supposed that the silver in contact with the sodium 
chloride forms a certain amount of silver chloride, which is photosensitive. 

Further evidence in favour of this hypothesis is this: the rocksalt 
specimens when subjected to a glow discharge turned yellow, in accordance 
with known expectations. It is found that yellowish specimens were 
spoiled much more quickly than the others. Fig. le referring to such a 
specimen, the picture being taken approximately half an hour after the 
exposure to the glow discharge followed by silvering. 

In terms of Mott’s explanation (1940) of Pohl’s experiments the 
formation of silver chloride can be explained in the following way. The 
glow discharge results in a stochiometric excess of the metal. As a conse- 
quence a certain amount of chlorine becomes available. When afterwards 
the silver is deposited some silver chloride will be formed, which has a 
low photosensitivity. If exposed only to low intensity daylight, the 
darkening takes place slowly, but under the influence of the very intense 
Hg source of illumination, used in interferometry, the silver chloride 
shows the “ print-out ” effect which causes the dark stains in the picture. 
This picture was taken a very short time after the specimen was silvered, 
it was completely clean when putting it in the light beam. Once in the 
interferometer the fringe system was kept fixed. As a result one can see 
‘the preferential darkening on one side of the dark fringes, i.e. on the 
side at which the light has a reflection maximum. [If this is indeed the 
correct explanation then the effect is highly critical since the maximum 
on the side of the fringe is but a few per cent above that of the intensity 
between fringes. 

This work has been carried out whilst holding a Fellowship from the 
“Centre National Belge de Chimie Physique Moleculaire ”’, and a visit 
to the laboratory at Royal Holloway College, London, was made possible 
with aid from the ‘‘ Fonds pour les relations culturelles Anglo-Belges a 
I am very much indebted to Professor S. Tolansky ; without his constant 
advice and encouragement this work would have been impossible. 
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[Received January 30, 1951.] 


ABSTRACT. 


The total cross-sections of carbon and hydrogen have been measured 
at neutron energies of 39, 64:5 and 97 MeV., using a recoil proton telescope 
as a neutron detector. At these three energies the hydrogen cross-sections 
are 223+7:6, 12643 and 73-9-+3-0x10-2? cm.? respectively. The 
corresponding carbon cross-sections are 1100+20, 784 +5 and 
508 +5 x 10-2? cm.?. 


INTRODUCTION. 


NEvTRON beams, produced by proton bombardment of internal beryllium 
targets at various radii in the 110-inch Harwell cyclotron, have been used 
to measure the n—p total cross-section at 39 MeV., 64:5 MeV., and 97 MeV. 
The method employed was similar to that used at 156 MeV. (Taylor, 
Pickavance, Cassels, and Randle 1950, 1951). 

The collimation of the neutron beam and geometrical layout of the 
apparatus differed little from the previous arrangements, but at each new 
energy a change was needed in the polythene radiator and carbon absorber 
of the triple coincidence telescope. The thicknesses of these are included 
in the summary of the results. 


EFFECTIVE ENERGY. 


As before, only the high energy end of the neutron spectrum was used 
(see Table), and the limits of detection were now such that not less than 
98 per cent of the recoiling protons counted came from the hydrogen 
content of the polythene radiator. The effective energy was determined 
in a subsidiary experiment at each energy by counting the number of 
recoiling protons as a function of absorber thickness in the telescope and 
hence as a function of energy. The maximum neutron energy determined 
this way was compared with the value expected knowing the energy of 
the internal proton beam and the geometrical arrangement. These 
results were used to check the range energy relation, and so an estimate 
of the error in the effective energy was obtained. 

To test for “ hardening ” of the neutron beam in the attenuator, the 
effective energy after passing through the carbon was determined in the 
course of the 64-5 MeV. experiment. Since this differed from the effective 
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energy of the unattenuated beam by less than 1 MeV., no appreciable 
hardening was taking place. The mean of the two values was taken as 
the effective energy. 


ATTENUATORS. 


The carbon and polythene 24-inch diameter bars were of such lengths 
that the polythene attenuator contained the same number of carbon 
atoms as the carbon attenuator and the attenuation produced by the 
hydrogen content varied between exp (0-5) and exp (0-8). Details of 
the attenuators are given in the Table. 


RESULTS. 


The results, including those at 156 MeV., are presented in the Table. 
The hydrogen total cross-section, oy, is plotted as a function of energy 
in the figure (p. 330). 

TABLE. 


Effective energy (MeV.) 39 
Estimated error in effective energy (MeV.) cen! 
Upper energy limit of neutrons in the beam 
(MeV.) 
Lower energy limit of detection (MeV.) 
Absorber in telescope (gm.cm.~? of carbon 
equivalent) ; 3-01 


Radiator thickness (gm.cm.~? of polythene) : 0-63 
Carbon attenuator (gm.cm.~*) 23° 70-76 
Polythene attenuator (gm.cm.~*) . 82-41 
Hydrogen total cross-section x 10°? (cm.”) 23 126 

; +3-0 
Carbon total cross-section x 10?” (cm.?) 784 


‘DISCUSSION. 


The value obtained for the hydrogen total cross-section at 97 MeV. is 
in good agreement with the most accurate experiment at Berkeley, in 
which the effective energy was 95 MeV. (DeJuren and Knable 1950). 
However, the result at 39 MeV. is greater than that expected from earlier 
Berkeley measurements at 40 MeV. (Hadley, Kelly, Leith, Segre, Wiegand, 
- and York 1949), but in good accord with the more recent determination 
at 42 MeV. (Hildebrand and Leith 1950). 

It may be of interest that a simple empirical relation can be used to 
fit the results within the experimental errors quoted, namely, 


ple TS Eeeccd CMH ADE 


On 


where oy is in units of 107°’ cm.2 and Ein MeV. This relation also gives 
agreement with the results at 14 MeV. (Salent and Ramsey 1940), 25 MeV. 
(Sherr 1945), and 280 MeV. (Fox, Leith, Wouters, and MacKenzie 1950). 
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ABSTRACT. 


The rate of production in ice of fast neutrons is given. It is shown 
that the number of fast neutrons produced at any point in the ice is 
approximately proportional to the number of cosmic ray stars at that 
point. The results are compared with those of other investigations, and 
it is pointed out that many previous results are too low due to the misuse 
of a formula given by Bethe et al. (1940). 


INTRODUCTION. 


DurRINnG recent years there have been many experiments to find the 
absolute rate of production of neutrons by cosmic rays, and the variation 
of this rate with the height of production. In general the results of 
these experiments have been in agreement only as far as the height 
variation is concerned, the values found for the rates of production varying 
by a factor as large as five. It has been found that the variation with 
height over the lower part of the atmosphere is exponential, the rate of 
production at a depth # gm./cm.? below the top of the atmosphere being 
proportional to exp (—a/A) where A is 140-150 gm./em.?. This height 
variation is the same as that found for stars observed in photographic 
emulsions. ° 

For measurements on slow neutrons, the use of boron loaded nuclear 
emulsion has a great advantage over counter or ionization chamber 
tech ique, in that the boron disintegrations 

B+ > Lij+ Hes 

can actually be observed and there is therefore no doubt as to what is 
being measured. In the electronic devices, on the other hand, elaborate 
precautions must be taken to eliminate background counts, and one 
can never be certain that no spurious events are counted, or that all 
true counts are recorded. The photographic plate, however, requires 
long exposures, of the order of 50 days, and cannot, therefore, be used to 
examine short term fluctuations, such as have recently been observed by 
Adams (1950) and others, using electronic devices, after intense solar 
activity. 


* Communicated by Sir George Thomson. 
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EXPERIMENTAL. 


The plates used were Ilford 100 C.2 boron loaded emulsion and were 
exposed vertically, at depths of 0, 0-5, 1, 2 and 3 metres below the surface 
of the glacier ice at 3600 m. elevation at the Jungfraujoch. The geometry 
was excellent, the plates being contained in an iron pipe, internal diameter 
21 inch., wall thickness } inch., which had been sunk into the ice. The 
glacier ice was covered with approximately 15 gm./cm.? of snow. 

The results are shown in fig. 1, where the numbers of boron disinte- 
grations observed (in arbitrary units) are plotted logarithmically against 
the depth below the glacier ice. As can be seen, the results are consistent 


Wag. io 


ie 
1-0 Ss 


log,, (No. of B. disintegrations). 


0 100 200 300 
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with an exponential absorption. The best values for the range in ice 
and the number of disintegrations at the surface of the ice were calculated 
by the method of least squares, and were found to be : 

Range in ice=126 gm./em.?. 

No. of boron disintegrations at the surface of the ice =893/c.c. of 
emulsion/day. 


INTERPRETATION OF RESULTS. 


In the following paragraphs we shall show how the number of neutrons 
captured at any point, which is the quantity we measure, is related 
to the number of neutrons produced at that point. 


& 


¥ 
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The boron disintegrations are produced by thermal and near-thermal 
neutrons only, since it is only in this region that the cross section is large. 
The neutrons will be slowed down to thermal energies by elastic collisions 
with the hydrogen in the ice. It seems probable that the neutrons will 
be produced in cosmic ray stars in the ice. We shall assume that (1) the 
neutrons are produced in stars, and that the mean number of neutrons/star 
is independent ofthe depth in theice. In this case, the neutron production 
will decrease exponentially with depth in the ice, with a range of 
130 gm./em.? (Harding 1950). (2) the neutrons can be divided into two 
groups, one group consisting of those neutrons that are evaporated from 
stars, the second of those neutrons that correspond to the high energy 
tail (Perkins 1947) observed in the proton spectrum. The evaporated 
neutrons will have an energy spectrum given by the formula 


P(E) dE=A aad Cx Ey bah ee ee LY 
where P(E) is the number of neutrons emitted with energies between 
EK and E+dH#H, A is the total number of neutrons emitted, and T is a 
constant, with the dimensions of an energy. From a study of the protons 
it has been found that for stars of three or more visible tracks, T is equal 
to 4MeV. It has further been found that the evaporated protons are 
emitted isotropically, (Harding et al. 1949), and we shall suppose that 
this is also true for the evaporated neutrons. Now, from formula (1), 
we can calculate that 96 per cent of the evaporated neutrons will have 
energies less than 20 MeV. We shall therefore assume that all those 
neutrons with energies below 20 MeV. are evaporated, and have an 

isotropic angular distribution, while none of those with energies greater 
than 20 MeV. are evaporated, but they are all strongly collimated 
downwards. We shall deal with these two groups separately. 


EVAPORATED NEUTRONS. 


Bethe et al. (1940) have shown that the mean square distance r 
travelled by a neutron in losing energy from E, to E, is given by _ 


‘ B, dE 
aM | Tha) gee eae Me 2 at?) 
E, A 


where M is the atomic weight of the moderating material and 1(E) is the 
mean free path at an energy E. If we put E, and Hy equal to 20 Mev. 
and leV. (thermal energies) and neglect the scattering with the oxygen, 


we find 
F 72=120 (gm./em.?)?. 
This value will be greater than the true value because of the neglect of 
the oxygen scattering. We should note that the value of 7? decreases 
as KH, decreases. 
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Bethe ef al. have also shown that the distribution, after the slowing 
down process of neutrons produced at a depth xp, is given by 


N(x) dv=n-!L-1 exp [—(a—a,)*/L?7]dz, . . . . (3) 
where N(x) da is the fraction of the neutrons that will be found between 
x and w-+da as thermal neutrons and L is given by 


oN 
L= (5) = 7 gm./cm.* for 20 MeV. neutrons. 


It must be emphasized that formula (3) holds only if the neutrons are 
emitted isotropically. 

Now, the number n’(a) of evaporated neutrons produced at a depth x 
in the ice can be written 

n'(x)=n'(o) exp (—a/A)/gm./sec. .. . 9. . . (4) 
where A=130 gm./cm.?. 

Using equations (3) and (4), we find that B’(x)dx, the number of neutrons 
of the first group, assumed all to have energies of 20 MeV., reaching thermal 
energies/sec. between depths # and x-+-dz is given by 

‘ ‘(0) exp (—a/A) /” 
B(x) dua=dau pee { exp (y/A) exp (—y?/L*) dy 
n'(o) exp (—2/A) 


=r aay 


bd 1 
exp (L¢/4a%) | exp (—y"/L*) dy. 
| Relief 
Now, for L=7 gm./cm.*, we find that the upper limit of the integral in 
equation (5) can be replaced by oo with less than 3 per cent error for all 
values of x>10 gm./cm.?. 
Then we find 
B’(w)=n'(o) exp(—w/A) exp (L?/4A?), 
and for L=7 gm./em.*, A=130 gm./cm.?, exp (L?7/4A?) is very nearly 
unity (0-999), so therefore we can write 
B'(x%)=n'(0) exp (—a/A)=n'(a)-.-. . . . (6) 
We have derived equation (6) by considering all the evaporated neutrons 
to be emitted with an energy of 20 MeV. In fact, the majority of the 
neutrons will have energies of only three or four MeV. For these lower 
energy neutrons, 7°, and therefore L, will be reduced. This will make 
the contribution between y= +10 gm./cm.* to the integral in equation (5) 
even larger than 97 per cent, and will make exp (L?/4\?) nearer unity, 
therefore equation (6) will still hold. : 


HicH ENERGY NEUTRONS. 


The scattering mean free path in ice for neutrons of energy greater 
than 20 MeV. is largely governed by the oxygen cross section, and further 
varies only slowly with the energy. For a 50 MeV. neutron the mean 
free path is of the order of 20 gm./em.”. Now, since the oxygen nucleus 
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is Massive compared to the neutron the scattering is very large, and 
therefore, after one collision the neutrons will have an isotropic angular 
distribution. We will take the mean energy of the fast neutrons as 
50 MeV. and suppose that all the neutrons have this energy. Let the 
number of fast neutrons produced at a depth x be given by 


m (zee (oyexp(—%/A)). . «9. > (7) 


where A is again equal to 130 gm./em.?. 

Then the number of fast neutrons produced isotropically at a depth x 
will be approximately equal to the number of fast neutrons produced 
from stars a distance 20 g¢m./cm.? above x, since these fast isotropic 
neutrons are produced after one collision from the fast collimated neutrons. 

We can now apply equations (2), (3) and (5) to these isotropic neutrons, 
and we then find that 


72 = 6,500 (gm./em.2)2, 


L=30 gm./em.? 
and that 
BAG) =e 20), 9 Mie. 22). 28) 


' provided «x is greater than 50 gm./cm.’. | Fe 
B’’(x) dx is the number of neutrons, which were originally fast collimated 
neutrons, that reach thermal energies between 2 and x+dz. For values 
of x less than 50 gm./cm.”, B’’() will not have the full contribution from 
the neutrons produced in the ice, but this will be compensated for by 
the entry of fast neutrons from the air. 
From this analysis, we can say that, for depths of more than 50 gm. /om.? 
of ice, the total number of neutrons reaching thermal energies/sec. at a 
depth x is given by 
B'() +B" (x)=n'(w) +n (w—20), 


and ithat for values of x between 15 and 50 gm./em.”, we will’ have the 
full contribution to B’(x), but a reduced contribution to B’’(x) from 
neutrons produced in the ice, which will be compensated to some gz tent 
by neutrons produced in the air. This view 1s strengthened by a sid) 
of fig. 1, which shows that all the experimental points lie on mens 
line, with a slope corresponding very closely to that expected for a value 
“Sa pos B(a)+B’ (x) for a=15 gm./em.’. Bethe et al. show 
that if p(x) is the number of boron: disintegrations/c.c./day at a depth 


x gm./em.”, then 
ple) = =~ —? (Bz) + BY (a))p, 


‘and D are the number of boron atoms/c.c. of emulsion and the 
ey th 1 : toms/c.c. of ice respectively, pj and py the atomic 
number of hydrogen atoms/c.c. es 


ZAZ 
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capture cross sections of boron and hydrogen respectively, and p the 
density of ice. It is assumed that no capture occurs in the oxygen. 
Putting 


ee 115 borons, 
a R. K. Adair (1950) 
o)=9:049 borons, 


we find 
B’(15)+ B’’(15)=20/gm./day in ice. 


Fig. 2. 


log,) (No. of neutrons/gm./day). 
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Pressures in atmospheres. 
This figure shows some results obtained for the rate of neutron production as a 
function of the atmosphere depth, both before (dashed lines) and after 
(full lines) the corrections given in this paper have been applied. 
D—Davis (1950). 
K—Korff and Hamermesh 1946. 
L—Lattimore (present results). 


M—Montgomery and Tobey 1949 (no corrections required). 
Y—Yuan 1946. 


If we assume that only the oxygen in the ice produces neutrons, and 
also correct for the absorption in the 15 gm./em.? of snow by multiplying 
by exp (15/130), we find | 

B’(0)+B"(0)=25/em. of oxygen/day 
at an elevation of 3600 m. Now, 
B(x) +B" (a4)=n' (x) +n’ (a—20) 
=n'(x)+1-15n'" (x), 
since n'’(x—20)=exp (20/130) n(x). (See equation (7).) 


Rate of Production of Neutrons in Ice by Cosmic Rays 337 


If we assume that the neutron and proton spectra are similar, 7’’(x)/n'(a) 
is of the order of 4, and we can therefore write 


n'(x)+1-15n'"(~%)=1-04(n'(a)+n'(a)), 
B’(0)+B’(0)=1-04(n'(0)+-n'"(0)), 
and therefore the rate of production of fast neutrons in oxygen at 3600 m. 
is given by 25 divided by 1-04 which is equal to 24/gm./day. 
Rate of production of fast neutrons =24/gm. of oxygen/day at 3600 m. 


DISCUSSIONS OF RESULTS. 


The value obtained for the absorption of the neutron producing primaries, 
namely 126 gm./cm.? of ice, is in excellent agreement with that found by 
Harding (1950) for the star producing radiation. The rate of production 
also agrees well with that found in carbon at the same elevation by 
Montgomery and Tobey (1949), who used a similar method to that used 
here, that is to say, they surrounded an ionization chamber with large 
masses of paraffin wax. These two values are, however, greater by a 
factor of about two than those given for air by several other investigators. 
As had been pointed out in another paper (Lattimore 1951), this is 
probably due to the misuse of a formula given by Bethe et al. (1940) 
which involves the nitrogen cross section/molecule. In many cases the 
cross section/atom has been used, which gives a rate of production of 
only half the true value. In fig. 2 is shown some corrected values, and it 
can be seen that the agreement between these and the results presented 
here is much improved. 
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ABSTRACT. 


The’ binding energies of the nuclei H®, He* and He? are calculated, using 
the standard variation method with a simple two-parameter wave- 
function of a new form. Purely central interactions of the Yukawa and 
exponential types are considered. It is found that, when the nuclear 
parameters are chosen to fit the low-energy scattering and deuteron data, 
reasonable values are obtained for the three-particle nuclei, but the 
a-particle binding energy is too high. ‘This result is in agreement with the 


calculations of Rarita and Present for the exponential potential, and those — 


of Svartholm for the Yukawa potential. In contrast with Svartholm’s 
method of calculation, the potential and kinetic energies of the nuclei are 
obtained explicitly. In the case of the «particle, the Yukawa potential 
yields exceptionally high values for these quantities, and results also in a 
very concentrated nucleus. The exponential potential, on the other hand, 
does not give rise to such high values for the potential and kinetic energies 


nor does it concentrate the nucleus to the same extent as the former 
potential f. 


§ 1. INTRODUCTION. 


NUMEROUS investigations have been made on the binding energies of the 


three- and four-particle nuclei, using a central two-body interaction with | 


various forms of potential well. The methods of calculation adopted in 
these investigations and the conclusions arising therefrom have been 


summarized by Rosenfeld (1948). The so-called consistency problem may 
be stated as follows :— 


(1) A set of parameters compatible with the deuteron and low-energy 
scattering data, and with the binding energy of the triton, yields too large 
a binding energy of the «-particle. 

(2) A set of parameters yielding the correct binding energy of H?, He* 
and He involves a longer range than allowed by the p—p data and gives 
rise to a bound !8 state of the deuteron. 

In this paper the problem is approached from the point of view of the 
first statement. 

The standard variation method is used with a new two-parameter wave- 
function of simple form for both the Yukawa and exponential ‘central 
interactions. ‘The wave-functions chosen give rise to relatively simple 


* Communicated by the Author. 
+ See note added in proof. 
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mathematical analysis, no approximations being necessary in the evalua- 
tion of the integrals involved. This is to be contrasted with the more 
laborious calculations of Rarita and Present (1937) for the exponential 
potential, using the variation and equivalent two-body problem tech- 
niques ; and the work of Svartholm (1945), using the much more tedious, 
but admittedly very accurate variation-iteration method. The conclusion 
reached from the calculations carried out here is given in statement (1), 
and is thus in agreement with that of Rarita and Present, and Svartholm. 
It might be remarked, however, that the Yukawa potential yields results 
for the «-particle, which, though agreeing with the latter’s conclusion, 
differ from those of Frohlich et al. (1947) for the longer range of the 
forces—corresponding to a meson mass of 210-220 m,—which they use. 
The value obtained here for the binding energy for this particular range is 
approximately twice as large as that of Frohlich. 

The total kinetic and potential energies for each nucleus are given 
explicitly in contrast with Svartholm’s calculations, where the nuclear 
parameters are evaluated to fit the experimentally known binding energies. 
It is found in the case of the «-particle that the kinetic energy for the 
Yukawa potential is exceptionally high; the excess binding energy is 
large ; and, moreover, the nucleus is very concentrated. The exponential 
potential, on the other hand, gives rise to a much smaller excess binding 
energy, more reasonable values of the kinetic energy and a not so concen- 
trated nucleus *. 

The calculations of Gerjuoy and Schwinger (1942) have indicated, of 
course, that the introduction of a non-central interaction of the tensor 
type would considerably reduce the excess binding energy of the x-particle. 
In fact, their investigations with a square well potential have yielded too 
small binding energies for the triton and the «-particle. However, if one 
considers the exceptionally large binding energy obtained tor the Yukawa 
potential with purely central forces, it is feasible that the introduction of 
tensor forces might reduce it sufficiently so that the experimental value 
might be obtained when the normal ranges of the forces used to fit the 
two-body data are considered. This problem is now being investigated 
with the wave-functions used in the present calculations. 

The wave-functions introduced in this paper are of a form which might 
be extended for systems of more than four particles, provided the space- 
part is so modified as to have the correct symmetry properties for the 
particles in p-states. 


§2. Brinpinc ENERGY FORMULA) AND WAVE-FUNCTIONS. 

The two-particle central interaction with a mixture of Wigner, Bartlett, 
Heisenberg and Majorana forces may be written in the form (Rosenfeld, 
1948) a | 

Voaner= V (112) {4p + 4,6 PE +a,7 07 +-4,,0 D627 72) 
Hue. La 7 7 


or 


Vana= V ("1 2) {Ww BP, BE, MEL. Is 


* See note added in proof. 
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where V(r,,) is the two-particle triplet potential, 7,, being the distance 
between the particles, and the other quantities have their usual meanings. 
Now the ground states of the nuclei H®, He® and He‘ are S-states, so that 
an interaction consisting ofa mixture of Majorana and Heisenberg forces— 
so-called exchange forces—yields the same binding energy as an equi- 
valent mixture of Wigner and Bartlett forces. For a wave-function 
symmetrical in space coordinates the formula for the total energy of the 
nucleus has the well-known form (Rosenfeld 1948), 


EO —(h?/2M) f pAb dv®+ (3/2). +) J V(rre)? dv Econ 


(1) 
for the three-particle system, and 


BO —(12/2M) [ 64S dv+3(1 +9) | V(rr9)d? de + Bou 
: (2) 
for the four- pore system, Bow q is the ratio of singlet to triplet well 
depth, 40= 5 A, and 4®%= 2 A;, where 4; is the Laplace operator. 


i=1 =1 
acting on the coordinates of particle whose position vector is r;; % and ¢ 
are the space parts of the wave-functions for the three-and four- -particle 


systems respectively ; M is the mass of the nucleon. Eg¢,,,; denotes the - 


Coulomb energy of the protons in the nucleus and will be treated in the 
usual way as a perturbation. 

The problem of the calculation of the binding energy involves the choice 
of a good trial wave-function with one or more parameters, subsequent 
evaluation of the integrals, and minimization with respect to the variation 
parameters. The Gaussian wave-function has been frequently used in 
calculations in nuclear problems, because of the great simplification it 
usually gives in the mathematical analysis. It has, however, an incorrect 
asymptotic behaviour, and is a very poor approximation to the correct 
wave-function for a three- or four-particle system, when a Yukawa 
potential is used. The following new wave-functions with an exponential 
type of asymptotic behaviour are introduced in the subsequent calculations 
and result in fairly simple mathematical analysis. 

(a) Three-particle system. 
Let 
p=BY exp {—a(rig tris tres)? h/{rigtrig trea}. s -. (3) 
be the wave-function for the three-particle system, where B12 is the 
normalization coefficient ; « and are parameters ; 71, 713 and 743 are the 
distances between the particles | and 2, 1 and 3, 2 and 3 respectively. 

The wave-function is of course a function of the relative coordinates of 
the particles and it may be simplified by a transformation. 

Let 

e=—r,+(r.+r5)/2, (1/3/2)(r,—r5), R=(r,+r,+45)/3, 


(4) 


Seat Qee 
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r,, r, and rz being the position vectors of the particles 1, 2 and 3 
respectively. Relation (3) then assumes the form 

pNy" exp {—+/2a(p%+r2)VMi(p2tr2yn, ss, (5) 
N7? being the normalization coefficient in the new coordinate system and 
p=lel|,7=|r|. 

In this system of coordinates it is easily shown that the operator 
—(#?/2M)4® becomes —(h?/2m)(A,+ Ap)—(h?/6M) Ag, where J,, Ap and 
A, act on the coordinates r, gp and R respectively ; andm=2M/3. NowR 
is the position vector of the centre of mass of the three particles, so that, 
since # or 7’ is a function of the relative coordinates, the expectation value 
of the kinetic energy of the mass-centre is zero. 

Equation (1) may then be written in the form 


HO — (h2/2m) J p'(4,+ Ap)’ dr dp+(3/2)(1-+q) J V(ry2)h? dr dp+Egou: 
(6) 
The Coulomb energy for two protons, that is, the He® nucleus, may be 
written in the form 
| Eoou=@ J (Urb drde, =... (7) 
é being the electronic charge. 
(b) Four-particle system. 

The form of the wave-function taken for this case is similar to that of ¢ 
in equation (3). 

Let . 

p=C? exp {—a(rig tris tia tra t1aat 734)? 3/ 


Bee Wee eo 1832 i 2 ie 
(Tot 1ig +14 +193 +124 +1 34}"> oe ES Sa ae} 


where C? is the normalization coefficient ; « and n are parameters; 719, 
etc. are the distances between the particles 1 and 2, etc. The trans- 
formation 


u=(r,+1,)/2—(r3+4,)/2, R=(r,+r,+4r3+r,)/4, | 
v=(r,—r,)//2, w=(t3—F,4)/V/2, 
where Yr, fy, f, and r, are the position vectors of the particles 1, 2, 3 and 4 
respectively, applied to (8) gives 
d’ =N#? exp {—2a(u2-+v?+w*)1?}/{u? +0? + w? }”, eto) 
N3” being the normalization coefficient in the new coordinate system, and 
u=|ul, v=|v], w=|w]. It is easily shown that equation (2) becomes 
E@= — (h?/2M) f 6'(4,+4,+4,)¢' du dv dw 
| +3(1+9) { V(rie)$? du dv dw+Hooy + (11) 
in the new coordinate system, the expectation value of the kinetic energy 
of the mass-centre vanishing as in the case of the three-particle system. 
Au, Av and Aw are the Laplace operators acting on the coordinates 
u, v and w respectively. 


(9) 
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The Coulomb energy for the two protons may be written in the form 
Eoou=e” J (1/712) du dv dw, ig TO CEP a) (12) 


e being the electronic charge. 


§ 3. Evayuation or Brypine Enercy or H® and He’. 


The normalization coefficient N¥? for the wave function (5) is easily 
determined. 


J JP dp de =N,(4n)? |" [exp {—lp*+19)"2}/(p*+r2)*"Jp¥r* dp dr= 1, 


where B=2,1/2«. 
Using the transformation p=R cos %, r=R sin ¢ it is found that 
N= (08-4 , 29-8n)/f78(5—4n) 13. 2 we (EB) 
For the above-mentioned wave-function, the kinetic energy of the nucleons, 


‘corresponding to the first term on the right-hand side of equation (6), is 
given by 


—(H2/2m)N, f fexp {—B(p?+r2)12}/(p2-+ 12/20] 
x [202 (p? +r?) — +/2u(5—4n)(p?-+7?)!?2—n(8—4n)] dp dr. 
This is simply evaluated, using the transformation p=R cos #, r=R sin ys, 
so that the kinetic energy is 
3 (5—n) he 


2° (i—n)(6—4n) ° M” 
or 
(5—n) ey? 14 
ql wean) OMe og Be ome 
where c’= /6a/«, (1/«) being the range of the nuclear forces. 

When the wave-function (5) is substituted in the expression (7) for the 
Coulomb energy, then, since from symmetry considerations 7,5 may be 
replaced by 7y, 

16 y/ 6e?a 16e?« 
EK ‘ =—_—_ eer o_O 
Onul = Si 4nja 1 a(S oSan imo ee 


where ¢=k/4/6a. 
(a) Yukawa potential well: V(r\,)=—A exp (—k?y9)/(K?' 49). 
The second term on the right-hand side of equation (6) represents the 


potential energy of the three-particle system apart from the Coulomb 
energy. Now, when the wave-function given by (5) is used, 


J Vers? dr dp=—AN, J exp {—B(p?-+r2)!2—xer9}/ {eer ya(2-+72)"} de de. 


Since from symmetry considerations 7,, may be replaced by 723, that is, 
27/4/3, the integral may be written in the form 


7/2 
(AN /«’) | cos® ys sin yf up| R**" exp (—yR) dR, 
0 
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where y=$+«' sin, x’=2x/\/3 and the transformation p=R cos ¢, 
r=R sin % has been applied. Hence, by elementary integration, the 
potential energy has the form 


72 — cos? % sin xs 


— ((24ln\(1+qAc'/(5—4n)} | (teain pm W - . (16) 


where c’=1/c= /6a/k. 
Thus, for the H® nucleus, the total energy is given by 


(5—n) h” on 2A +) Ac’ (7? cos? ys sin yp 
4(1—n)(5—4n) ° Ma m(5—4n) a (1-Fe sin ps) 


9 
4 


dep. 


(17) 


The total energy for the He® nucleus is obtained by the addition of the 
Coulomb energy term given in equation (15). 

The expression (17) is now minimized with respect to the parameters c 
and n for different values of the range (1/«), with the corresponding values 
of A and q fixed by the low-energy »—p scattering and deuteron data. 
The integral appearing in the expression is easily evaluated for the values 
n=), 1/4, 1/2 and 3/4. These are given in § 1 of the Appendix. Minimi- 
zation with respect to ¢ is carried out for each of the values of n mentioned. 
It should be noted that n must be less than or equal to 3/4 for convergence 
of the integral representing the kinetic energy. 


TABLE I, 


Minima of E for H? and He? with respect to ¢ for each value of n. 


K é E \2 A K ee Hain Ecou He?/Emnin 
we E (MeV.) (MeV.) (MeV.) (MeV.) | (MeV.) 
0) 0-38 52-3 —60-2 =7-9 1-10 —6-8 

1/4 0-46 56-5 — 65-0 85 1-13 —7-4 
1/2 0-70 46-3 —53-3 —7:0 0-99 —6-0 
3/4 1-60 25-1 =27:2 = 2-2 0-65 —1:5 


Table I. gives the minimum values of E for H® for each value of n 
considered, for the range (1/«)=1:17< 10-8 em., which corresponds to a 
meson mass of 326m,, A and q having the values 67-3 MeV. and 0-69 
respectively. The corresponding values of the kinetic and potential 
energies, the Coulomb energies for the He® nucleus and the resultant 
values of the binding energy for the latter nucleus are also tabulated. 

The values given in Table I. indicate that the maximum value of | E | of 
approximately 85 MeV. arises for n~1/4, and c=0-46, that is, 
I/x=1-318X10-% em. The value of (1/x) gives an indication of the 
“ size’ of the nucleus. The maximum value of | E| for He® is 7-4 MeV. 
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The experimental value of the mass-difference of the nuclei H® and He’, 
that is, the Coulomb energy of the protons, is accurately known, being 
77146 keV. (Tollestrup et al. 1949). Comparison of the theoretical 
value with this gives a check on the accuracy of the binding energy 
calculations. From Table I., it is seen that the theoretical values of the 
Coulomb energy are rather high, though the actual values of the binding 
energies of H® and He® are in good agreement with the experimental values 
of (8:51 -+0-091) MeV. and (7:739 0-1) MeV. respectively. 

Svartholm, using the variation-iteration method, has calculated the value 
of the parameter A to fit a triton binding energy of 8-33 MeV. for the 
ranges 1-185 10- em. and 1-835 10-% em. Brown (1939) has also 
fitted the binding energy of the triton for the latter range, using the 
variation method with a large number of parameters. Table II. gives the 
binding energies for these ranges with the new wave-functions, using the 
values of g and A corresponding to Brown and Svartholm’s calculations. 
Considering the simplicity of the present calculations the results obtained 
are reasonable. 

TABLE II. 


Values of the binding energy of H® with new wave-functions corresponding 


to different sets of nuclear parameters. 


1/k Binding-energy 
(10-28 em.) ; of H? (MeV.) 


1-185 dl: 0:75 6-1 


(Svartholm) { 1-835 29. 0-58 8-4 


(Brown) .. “S3E 32-( 0-58 
(Assumed) . 37° 0-69 


(b) Haponential well: V(r,.)=—A exp (—Kryo). 
The analysis for this potential is of the same form as that for the Yukawa 
case. The potential energy reduces to the expression 
| 7/2 cos? ys sin? yb 
—(24/7)(1 A ————— dy 
(24/7)(1+-q) 9 (fesin yea y. 2 2), 4, eee 
The integral in (18) may be evaluated by elementary methods for n=0, 
1/4, 1/2, 3/4. [The values for n=0 and n=1/4 are given in §2 of the 
Appendix.| Thus, for the H*® nucleus, the total energy is given by 
“7/2 cos? x sin? ys d 
o (1+e sin f)e4" Ms 
: (19) 
The values taken for the nuclear paramerers are: (1/x)=0-865 x 10-Bem., 
A=123-3 MeV., q=0-6. The minimum value of E in this case occurs at 
approximately n=0, c=0-69, the kinetic and potential energies having the 


(5—n) 


: ih Aan | 
i= 4(1—n)(5—4n) 5 M KC “— (24/7) (1 +a)A | 
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values 29-1 MeV. and —36-7 MeV. respectively. Thus the minimum value 
for EK of —7-6 MeV. for H? is reasenably good. The Coulomb term given 
in the relation (15) has the value 1:02 MeV. at this minimum, so that the 
binding energy of the He? nucleus is approximately 6-5 MeV. The value 
of (1/x), corresponding to c=0-69, is 1-462 10- em., indicating the 
“size” of the nucleus. It is to be noted that there is a considerable 
difference in the values of the kinetic energies for the Yukawa and expo- 
nential potentials at the minimum of E*. 

Rarita and Present (1937) have fitted the value of the binding energy of 
the triton, using the exponential potential and the same values of the 
nuclear parameters as above, but a much more complicated trial wave- 
function than that used here has to be introduced. 


§4. EVALUATION OF THE BINDING ENERGY oF He’. 


The normalization coefficient Nj? for the wave-function 4’, given by 
equation (10), is easily determined 


[ 6? du dv dw=N,(47)? i [exp {—B(u?+-v?+ w?) 12} /(w? +2 + w?)2"] 
0 


x wv? du dv dw=1, 


where B=4«. 
Using the transformation 
u=R sin 6 cos ¢, v=R sin @ sin 4, t= Kicos'G, ~~ (20) 


it is found that | 
Np=(105 . 213-8" | 9-4n)/fr4(8—4n) 1}. 2... (21) 


For the above mentioned wave-function, the kinetic energy of the nucleons, 
corresponding to the first term on the right-hand side of equation (11), is 
given by 
— (h2/2M) f [exp {—B(u? +0? + w2)12}/(u2 + v2 w2)n tt] 

— & [402(w? + v2? + w*) —20(8—4n)(u?+ v?+ w?)!/2—n(14—4n)] du dv dw. 


This is simply evaluated, using the transformation (20), so that the kinetic 
energy iS 


GS / 2 
ee si Cs (22) 
4(4—2n)(7—4n) © M 
where c=2,4/2a/k. , 
When the wave-function (10) is substituted in the expression (12) for the 
Coulomb energy, this may be evaluated, using (20), to give 


354/ 2e7%a 35e?cK 


yp ee Ue (23) 
Conl ~ 32(2—n)  64(2—n) 


E 


* See note added in proof, 
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(a) Yukawa potential well: V(r,3)=—A exp (kj 2)/(K 12): 
The second term on the right-hand side of equation (11) represents the 
potential energy of the four-particle system apart from the Coulomb 
energy. Now when the wave-function ¢’ given by (10) is used, 


{ V(7r,2)6? du dv dw 

=— AN, | exp {—B(u2+02-++-w?)!2—xryo}/ {ar a(u? +0? +0?)* } du dv dw. 
Since from symmetry considerations r,, may be replaced by r,4, that is, 
\/2w, the integral may, with the help of the transformation (20), be 
written in the form ; 

~71/ _ 

0 


n/2 2 me é 

—(AN,|«’)(477)? | sin? 6 cos? dd | sin°d cos6 dé R74" exp (—yR) dR, 
0 0 

where «’=1/2« and y=fB-+.«’ cos#. Hence, by elementary integration, 


the potential energy has the form 
315(1+¢q)Ac*” is u(l—u?)? 
~~ 82(2—n) Jy (wFoyea 
where w=cos 0 and c=B/k’=(24/2a/k). 
Thus for the He* nucleus the total energy is given by 
ee en) he 
~ 4(4—2n)(7—4n) ° M 
315(1+q)Ac? 4” 7! w(1l—u?)* 
32(2—n) ih (utes 
Kcom being given by (23). 

The minimization is carried out in the same way as in § 3 (a) for the H® 
nucleus, the Coulomb term being omitted in this stage, and then added 
later as a perturbation. The values of n=0, 1/4, 1/2, 1, 5/4, 3/2 are 
considered ; convergence of the integral for the kinetic energy imposes the 
condition that n<3/2. The values of the integral appearing in the poten- 
tial energy term in equation (25) are given in § 3 of the Appendix. The 
results for the values (1/«)=1-17 x 10- em., A=67-3 MeV., g=0-69 are 
given in Table III. 


dase! Sika Ee Teme 


K 


R202 


bith 


TABLE III. 


Minima of E for He* with respect to ¢ for each value of n. 


n C K.E Pole Eyin Eeou eee 

(MeV.) (MeV.) (MeV.) (MeV.) (MeV) 
() 5-0 188-9 — 243 +5 —54:6 1-97 52-6 
1/4 45 200-4 —255-4 —55-0 1-73 53:3 
1/2 3-7 176-6 — 230-8 — 54-2 1-66 52-5 
l 2-3 173°5 — 220-0 —46-5 1:55 44-9 
5/4 1-5 139-0 = 1744) —35-0 1-33 33-7 
3/2 0-63 75-0 — 88-4 —13-4 0-85 12-5 
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The values given in Table III. show that the parameter n for values from 
0 to 1/2 does not influence the value of | E| very much. The maximum 
value of the binding energy is approximately 53-3 MeV. for n=1/4, c=4-5, 
that is (1/«)=0-735 x 10-% cm. Thus a very large excess binding energy 
is obtained, proving the inadequacy of the central Yukawa interaction. 
The kinetic energy of the nucleons is exceptionally high and the nucleus 
itself is very concentrated. The Coulomb energy too is rather large. 

The minimum values obtained for the binding energy, for the nuclear 
parameters of Svartholm and Brown, are given in Table IV. In addition, 
the binding energy is calculated for the parameter (1/)=1-185 x 101 em., 
q=0-75, A=56-0, assuming this to be the approximate convergence limit 
for A for the iteration method. Svartholm gives only the value obtained 
by the first iteration. It is also shown that the very long range of 
1-9 10~¥ em. has to be taken to obtain a value, which is fairly near the 
experimental binding energy of 28-2 MeV. It might be mentioned that 
the value obtained by Frohlich et al. (1947) for a comparable range is only 
half of the experimental.value. 


TABLE [V. 


Values of the binding energy of He4 with new wave-functions corresponding 
to different sets of nuclear parameters. 


lik A Binding energy 
(10-13 cm.) q : (MeV.) 

61-01 * 43 

(Svartholm) 1-185 0-75 
60-034 + 40) 
x5 1-185 0-75 56 (assumed) 27 
33°80 * 39 

ex 1-835 0-58 
29-69 23 
(Brown) 1-835 0-58 32-06 * 32 
(Assumed) 1-9 0-57 29-3 § 27 
(Assumed) [7 0-69 67:3 § 53 


‘ i] 
* Depth of well obtained by variation-iteration method to fit binding energy 
of H?. Apetiels 
+ Depth of well obtained in one iteration to fit binding energy of Her 
+ Depth of well obtained by two iterations to fit binding energy oLHe-, 
§ Depth of well to fit deuteron and low-energy scattering data. 


Thus the simple wave-functions introduced in this investigation yield 
results in close agreement with the calculations of Svartholm for the 
Yukawa potential. 

(b) Haponential potential well: V(712)= —A exp (—krj9). 

The analysis for this potential is of the same form as that for the Yukawa 

case. The potential energy reduces to the expression 


a 


—(315/8)(1+q)Ac?-*™” | {u2(1—u?)?/(ut+e)?*"} du: ee 26) 
0 
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The integral in (26) is of a similar form to that in (25). Its values for 
n=0, 1/4, 1/2 are given in § 4 of the Appendix. The total energy of the 
He? nucleus is then given by 


ae (28—2n) h? Be 
1 w(1—w)? “ 
— (315/8)(1-+q)Ac?-™” » wpeman MU Bou ‘Lyte wine) 


Eoou being given by (23). 

The values taken for the nuclear parameters are: (1/«)=0-865 x 10-¥em., 
A=123:3 MeV., g=0-6. The maximum value of the binding energy 
occurs at approximately »—0, c=2-25, the kinetic, potential, and Coulomb 
energies having the values 70 MeV., —102-9 MeV. and | MeV. respectively. 
_ Thus the maximum binding energy is approximately 31-9 MeV., an excess 
of 3-7 MeV. over the experimental value of 28-2 MeV. The value of (1/«) 
corresponding to c=2:25 is 1-087 x 10-1 cm. Thus there is a considerable 
difference in the effect of the Yukawa and exponential central interactions 
on the binding energy of the «-particle. The latter potential gives a much 
smaller excess binding energy, a much lower kinetic energy, and does not 
concentrate the nucleus so much as the former potential *. 

The conclusion that the exponential potential leads to too much binding 
energy for the «-particle is in agreement with the calculations of Rarita 
and Present, though their method results in a slightly larger value for the 
excess binding energy. 

Since the above calculations confirm that central Yukawa and expo- 
nential interactions are inadequate for a description of the binding energy 
of the «-particle, consistent with the binding energy of the three-particle 
system, the deuteron and the low-energy scattering data—the same 
conclusions have been reached by Svartholm, and Rarita and Present for 
the Gaussian interaction—the question of non-central forces should be 
considered. A tensor force is of course necessary to explain the quadrupole 
moment of the deuteron. Gerjuoy and Schwinger, using a square- -well 
potential with a mixture of central and tensor forces, have come to the 
conclusion that tensor forces are not very effective in binding the triton 
and «-particle. More recent calculations on the binding energy of the 
triton (Rarita and Feshbach, 1949, Clapp, 1949, Hu and Hsu, 1950) have 
shown that it is possible to obtain fairly good agreement with the experi- 
mental binding energy with a mixture of central and tensor forces when a 
Yukawa potential is used. It is still an open question whether such a 
mixture can account for the binding energy of the «-particle when the 
parameters are chosen to fit the deuteron and low-energy scattering data. 
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* See note added in proof. 
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APPENDIX” 


The values of the integrals listed in the following sections are used in 
the text :— 


7/2 cos? ab sin ys 
§ ib. I, (1+e sin p)>-™ dis. 
0) (1/24) {(8 + 9e2—2¢4)/(1 —c?)®—15¢ cos} ¢/(1—c?)72}, Gaal. 
n=1/4: (1/6){(2+c?)/(1—c?)?—3c cos-! e/(1—c?) 5/2}, Oma Ee 
n=1/2: —(1/2c*){—c/(1—c?)—(2—3c?) cos-1 CHE eye c<l. 
m=3/4: (1/c*){m—2c—(2—c?2) cos-1 ¢/(1 Sia ba Cal. 
(1/c?) {7 —2c—(2—c?) log {e+ (c?—1)¥?}/(c21) 12}, Cel. 
72 cos? x sin? y 
§ 2. I. eden pee 
m=0: (1/120){15(1+ 6c?) cos-1 e/(1—e?)°? —¢(81 +28¢?—4c4)/(1 —c2)4}, 
Osll.. 
n=1/4: (1/24){3(1+4c?) cos-} o/(1—e?)*? —c(13-+-2c?)/(1—c?)3}, C= 1) 


a 
$3. i mL —a2) (abe) 4 dau 


m=0: (35¢3+47c?+25c+5)/{210c8(1+-c)5}. 

n=1/4: (5c?+4c+1)/{30c5(1+c)*}. 

n=1/2: log {(1+¢)/c}—(15c?-+36c+23)/{15(1-+-c¢)?}—(2c2—1)/(20c?4). 
m=1: .(10c?—2) log {(1+c)/c}+5+1/6c?—10c—4/{3(1--c)}. 
n=5/4: 10c?+1/2c—5ce—8/3—(10c?—6c) log {(1+e)/c}. 

n=3/2: (5c*—6c2+1) log {(1 }c)/c}+-5c?/2+13¢/3—5c2— 7/4. 


: 9—4n 
§ 4. Ip u?(1—u?)?/(u+c)>-™ du. 


n—0: (64c?+ 69¢c?-++30c-+5)/{840c%(1-+c)*}, 
n=1/4: (8c?-+5c+1)/{105c5(1+c)>}, 
n=1/2: log {(1+e)/e}—(4c?—1)/(60c*) 
— (30c?+103c2-++122¢-+51)/{30(1-+c)#}. 
Note added in proof.—Since this paper was written, the author has 
investigated the effect of changing the range in the case of the exponential 


potential to the value given by a recent analysis of the low-energy p—p 
scattering data by Jackson and Blatt (1950). The values for the nuclear 


parameters are then: | 
(1/<)=0-709 x 10-8 cm., A=173-0MeV., g=0-64. 
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The method of calculation of the triplet well depth is given by Blatt and 
Jackson (1949). 

For H8 the maximum value of | E| is approximately 9-5 MeV., corre- 
sponding to n~0 and c=0-72, that is, 1/c<=1-25x10-%cem. The binding 
energy of He is then 8-5 MeV. Thus an excess binding energy is obtained 
in this case for the three-particle nuclei. The kinetic energy of H? is 
39-7 MeV. For He‘a binding energy of 47 MeV. is obtained, corresponding 
to n~0, c==2°32, that is, 1/x=—0-864 x 10-13 cm., the Coulomb energy of 
1-3 MeV. being taken into acount. Thus the shorter range exponential 
potential of Jackson and Blatt yields a much greater excess binding 
energy for the «-particle than the range used by Rarita and Present. 
The excess is in fact comparable to that obtained for the Yukawa 
potential. The kinetic energy of 106-2 MeV. obtained for the «-particle 
for this case is much more reasonable than the value given by the Yukawa 
potential. 
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ABSTRACT. 


A method is given for finding the equilibrium positions of a set of like 
dislocations in a common slip-plane under the influence of a given applied 
stress. Their positions are given by the roots of a certain set of orthogonal 
polynomials. The case of a set of free dislocations piled up against a fixed 
dislocation by a constant applied stress is discussed in detail and the 
resulting stress-distribution is compared with that produced by a crack 
with freely slipping surfaces. 


$1. INTRODUCTION. 


THE following problem (Frank 1951, Kuhlmann 1951, Nabarro 1951) 
arises in interpreting the plastic behaviour of solids in terms of dislocations. 
A set of identical dislocations lie in the same slip-plane. What positions 
will they take up under the combined action of their mutual repulsions and 
the force exerted on them by a given applied shear stress, in general a 
function of position along the plane ? 

Since the dislocations (assumed to be of infinite length and parallel to one 
another) repel each other inversely as the distance between them the 
problem is unchanged if they are replaced by a set of line-charges and the 
applied stress is replaced by an electric field. This electrostatic problem 
was used by Stieltjes (1885) (cf. also Szegé 1939) to illustrate the properties 
of the zeros of orthogonal polynomials. In the present paper we show how 
the properties of the classical orthogonal polynomials may be used to 
discuss the dislocation problem. 

Stieltjes solved the problem by minimizing the potential energy of the 
charges. We start from the idea of the force acting on a dislocation. 


The two methods are equivalent. 


$2. 

Suppose that there is an infinite straight dislocation parallel to the 
z-axis, having the plane y=0 as slip plane, and passing through the point 
a=a, Consider the shear stresses p,,, P,, in the slip-plane. If the 

* Communicated by the Authors. 
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dislocation is of pure screw type pz,=0, P,,40. If it is of pure edge type 
Pry, Py,=9. In either case the non-vanishing stress component which 


it produces at the point x in the slip-plane is 
A 


L—2; 

For a screw dislocation A=A,=b/27 and for an edge dislocation 
A=A,=pb/27(1—o) 

(Burgers 1939). Here » is the shear modulus, o is Poisson’s ratio 

and 6 is the magnitude of the Burgers vector b which gives the change 

of the displacement vector on passing once round the dislocation line. 

For an edge dislocation in an anisotropic material A=Kb/2z7 where K is a 

certain function of the elastic constants (Eshelby 1949). 

A screw dislocation is acted on by a force bp,, per unit length, where p,, 
is the total stress acting at its centre, excluding that produced by itself. 
Similarly an edge dislocation is acted on by a force bp,,. If there are n 
dislocations at the points x, v2... . x, all of the same type and each one is 
in equilibrium, the equations 


(1) 


p= 


a eee 
vaaa (2,0, J=1f2/ tne 
wFD 


must hold. P(x) is the appropriate component of the applied stress at the 
point z, i.e. the xy component for edge and the yz component for screw 
dislocations. 

Equations (1) and (2) apply, more generally, to any set of parallel like 
dislocations lying in their common slip plane. If is the angle between 
the line and Burgers vector of each dislocation, we then have 


A=A, cos /-++A, sin x, 
P(x)=P,.(«) cos $+P,,,(x) sin yf, 


since the screw and edge components interact independently with one 
another. 

In all cases, therefore, the following condition must be satisfied at each 
dislocation : the component in the direction of b of the traction on the 
plane y=0 due to the other dislocations and to the applied stress must 
vanish. 

It will be convenient to choose the unit of stress so that A=1. The 
condition of equilibrium is then 

x 1 


igj Uj—%; 


te )=0, j==1)2.) 2.7 ee et a 
We may regard the a; as the zeros of the polynomial 


f= Il (x—a;). eee <seeera err ( 4) 


14= 


i 
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f has the convenient property that its logarithmic derivative is equal to the 
stress due to all the dislocations : 


u n 
ar ee lew 
fo gn Ba; 
The stress when the jth dislocation is missing is 
feel 
fo &— 3 
The value of this expression when «=a; is 


(w—ay) f'(#)—f(@) 


lim Baill li "(2 ) 
toa em aefe). Fey)’ 
by double differentiation of numerator and denominator. 
The conditions (2) can thus be written 


f(x;)=0 


f"(x;) ee ee 27) 
iejuMat Wee 


Consider the differential equation 


f(a) +2P(a)f"(a)+4(n, x)f(z)=0. - . 6 ee + (8) 


Suppose that we can choose qg so that this equation has a polynomial 
solution of the nth degree all of whose roots are real and distinct. Then if 
q does not have a pole at any of the roots the conditions (7) are satisfied and 
the problem is solved. 

In the physical problem we may be interested in cases where certain 
of the dislocations are “locked ” in fixed positions (Cottrell 1948). If 
a dislocation is locked at «=, we imply that it lies in a stress field 
p(x)=const.8’(~—x,) where 8’ is the derivative of Dirac’s 6-function. 
The dislocation is in equilibrium at =x, where p(x)=0. Also, since p’(x,) is 
negatively infinite, the equilibrium is stable and an added stress will 
shift the position of equilibrium only infinitesimally. The other 
dislocations are unaffected by p(x). It is then convenient to omit p(2) 
from the expression for the applied stress and to regard one of the 
dislocations as “locked ” at the point 2,. 

Suppose, then, that in addition to the n “ free ”’ dislocations there are 
locked dislocations at 41, %p42-+--y. The equilibrium positions of 
the free dislocations are found by regarding the stresses produced by the 
locked dislocations as forming part of the applied stress, so that (8) 


becomes 


(5) 


(6) 


pr+2{P(e)+ ees tan, 2/0, eect 42 1(9) 


a=n+1U—% 
In this case we expect q to have a pole at each x, in order that f'/fshall not 
vanish there. 


f 
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The solution will be of the form (4) and f’/f will give the stress due to the 
free dislocations alone. If we make the substitution 


v 


F(a)=f(z) I (x—«,) 


(9) becomes Sais 
SPP + re, —F=0, a Lee 
where 
Tera) HP ae, — 2 aaa 


and «, 8 take the values n+1,n+2...v. At the zeros %, 7%... 2%, the 
conditions (7) are satisfied (free dislocations) whilst at the zeros 2,41, 
Xn, ---+Xy the last term in (10) has a finite value owing to the vanishing 
of a factor in the denominator, and the second of the conditions (7) (with 
F replacing f) is not satisfied (locked dislocations). [F’/F is the stress 
produced by locked and free dislocations. 

If a particular problem has been solved in this way it remains to verify 
that the arrangement of dislocations is stable. The condition for this is 
that the «-derivative of the total stress at the position of a free dislocation 
(excluding its own stress) shall be shes es d. e. that 
d E'(x) | 
After some calculation this reduces to the “eondisen that 

I(x) =Q—P?—P’>0 «$2 uenet oa’ an an 


for =a, %,...2,. I(x) is the “ invariant ” of equation (10) which can 
be reduced to the form 


vo" +1(x)v=0 ee ea rs) Ee 
by the substitution 


v=F exp [Paz). Oo prs 


v has the same zeros as F with possibly additional ones. A region in which 
I<0 can usually be dealt with by the following argument. Suppose that v 
has a zero (which may be at infinity) in the region. Since v” has the same 
sign as v the curve v=v(«) is convex to the a-axis and there can be no other 
zero in the region. In particular, if <0 for 2 >2’ and v0 as aco then 
F can have no zero for 2 >2’. 

It is easily seen that v'/v gives the total stress due to free and locked 
dislocations and the applied stress. 


§3. 
The general method will be illustrated by some cases of physical interest. 


(i.) A row of n dislocations under zero applied stress, the outer two being 
locked at z=+L. If L—1, (9) becomes 


fafa + aah ream a) fo, 


ee 
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or if we take 
n(n—1)—2 
l—a? ”’ 


(1—a*)f”—4af" + {n(n—1) —2}f=0 


which is satisfied by f=P),_,(x), the first derivative of the (n—1)th 
Legendre polynomial. Hence 


fHP%_a(e/L), F=(L?—22) f=P%,_,(a/L)a/(L?—2). 
In ordinary stress units the stress on one of the locked dislocations is 
+AF"(1)/F’(1)= +4n(n—1)A/L. 


The roots of f up to n—1=8 have been given by Kopal (1946). 

(ii.) A row of n free dislocations with an applied stress proportional to x. 
Such a uniformly varying stress exists near the centre of a beam supported 
at the ends and uniformly loaded. We have P=a(dP/dx) with dP/dx 
constant. The resulting force on a dislocation can be regarded as arising 
from a potential energy V=—4 ba? (dP/dx). “ Potential troughs” of 
this type may exist in real crystals owing to the presence of various 
defects (Mott and Nabarro 1948). If we choose 1/4/(dP/da) as the unit of 
length the equation for F(=/) is 


F’—2aF’+QF=0. 


If Q=2n, (10) is satisfied by the nth Hermite polynomial, H,,(z). The 
invariant I=2n-++ 1—2? is positive if |7|<4/(2n-+1) and negative otherwise. 
The reduced equation is satisfied by v=H,,(x) exp (—42"), so that by the 
argument following equation (13) all dislocations lie in the region 
|a|<v/(2n+1). (cf. also Appendix, § A2). 

In ordinary units of length and stress 


ee) 


and all dislocations lie in the region 


lel f ara | 


(iii.) A set of n—1 free dislocations along the positive x-axis in equilibrium 
under the combined effect of a locked dislocation at the origin and a 
uniform stress tending to move them in the negative direction. This case 
is discussed in detail in § 4. 

(iv.) Two locked dislocations as in (i.) with an applied stress as in (ii.). 
This leads to an equation occurring in the theory of the hydrogen molecular 
ion. There is an extensive literature (Wilson 1928, Teller 1930, Hylleraas 
1937) but no general discussion of the polynomial solutions. As we should 
expect there are n polynomial solutions of order n, since we can put 0, 1, 
2 _..n dislocations in the interval —L<a#<L and the remainder in the 


interval L<z. 
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For case (iii.) above we have P(x) constant and equal to say —7p. 
Equation (8) becomes 


f’+2(a—19) f’ +af=9, 
If we take 1/27, as the unit of length and put g=(n—1)/a this becomes 
af’+-(2—2)f'+(n—1)f=0 .. .: . | (14) 


which is satisfied by the first derivative of the nth Laguerre polynomial, 
L/(z). Thus in ordinary units of length and stress 


fHLi(279x/A), F=(279%/A) Ly, (2797/A). 


The spacing of the dislocations is the same as that of the radial nodes of a 
hydrogen atom in the ns state. 

The invariant of (14), I=nx1—l is positive if x<4n. The reduced 
equation is satisfied by v=aL/(x) exp (—4) so that all zeros of f must 
occur for x<4n. A lower bound for the greatest root is obtained in the 
Appendix, §A.1. Since 


n—1 n! (—a)* 
Sd witch ETP TERR WUT PE GY, 
there are no negative roots. 


When z is small compared with n we may neglect } in comparison with 
n/x and the reduced equation v”-+ _Ilv=0 becomes 


v’ +(n/x)v=0 


0= 47m Ji {4/(4na)} . ee 


so that near the origin the position of the ith dislocation is given by 
«=j;/4n where j; is the ith zero of the Bessel function J,. The position of 
the free dislocation nearest to the locked dislocation is 2,~3-67/n. Since 
j;=77 for large i the number, 7, of dislocations between the origin and the 
point x is 


(15) 


with the solution 


t= (2/7) ia. 1<i<n. 


Returning to ordinary units of stress and length in which 2 becomes 
27)x/A we have the following approximate results for large n : 
L, the length of the slip plane occupied by dislocations 


2A fry mis os Sh a ee TT) 


the number of dislocations lying between 0 and x 


==) (7). ee ees 


d, the distance between the locked and nearest free dislocation 


saLiBAUA ery) atk i, <x Foetal, 7019) 
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The difference between the L of equation (17) and the true distance 


between the locked and furthest free dislocation, L’ say, can be estimated 
from the inequality 


1-47n-28—0-55n-58 < eel <2-56n-28, 


obtained by combining the result of § A/1 with an inequality given by Szegé 
(1939, equation 6.32.6). In what follows we shall refer to L (and not the 
accurate L’) as the length of the row of dislocations. Equation (18) is 
true in the sense that if 1<i<n the ratio of the right-hand side of (18) to 
the correct i approaches unity as n tends to infinity. Actually (18) is in 
error by a factor of only 4/7 even if i=n and in error by only one or two 
units for small ¢ provided n is large. d must be larger than a few atomic 
distances if our results are to be valid, since the expression (1) is inaccurate 
very close to a dislocation. (Nabarro 1940). 

If we move along the z-axis the discontinuity 4 in the displacement on 
crossing the slip plane changes by b each time a dislocation is passed. The 
curve 4— A(x) is thus stepped, but according to (17) and (18) it will 
approximate to the smooth curve 


2br 
A= — V (La), eee. te.” (20) 


provided v is large and we do not go too far from the origin. 

From the coefficients of x" and x”? in (16) we find (27)/A)X’%,=n(n—1) 
so that the centre of gravity of the free dislocations is at the 
point <=nA/27,=}L whilst the centre of gravity of the n dislocations is at 
x=(n—1)A/27,. The coefficients of x and x give AXa,1=(n—1)r». 
According to (1) this is the stress at the locked dislocation produced by the 
free dislocations. The total stress on the locked dislocation is thus 77». 
This also follows from (14) which gives f’(0)/f(0)=3(n—1) which must be 
multiplied by 27, to convert to ordinary units. Cottrell (1949) has already 
obtained this result by a virtual work argument. 

The behaviour of the stress, 7, beyond the locked dislocation (x <0) is 
also of physical interest (Frank 1951). To discuss it, it is convenient to 
reverse the sign of x so that the dislocations lie on the negative a-axis. The 
most important region is that in which « is of the order of L, or, in the 
reduced units in which 27,/A=1, a~n. In these units the solution (16) 
with the sign of « changed is 


v=const. y I, (y), y=v/(4nx) 


and this is a good approximation when «<4n. The total stress is given by 


v/v, for which we find 
v! m\ In(y) 
= 4) (=) Tw)’ u<4n. 
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From known properties of Bessel functions I,/I,—>l as y> © and 
I,/1, > 2/y as y> 0 so that v'/u=/(n/x) if 1/4n<a<4n and v'/v=1/x if 
x<1/4n or in ordinary units 
[t= V(L/x), a<a<L ree a ok ob 
7T=A/a, pd <x. 


Fig. 1. 


fey) 


om Ct 


0 0-5 oe oF 40 | 
Shear stress in front of a row of held-up dislocations. Upper and lower bounds 
from equations (22). © exact values calculated from the polynomial (15) 


for n=10. For comparison the stress due to a crack of length 4L is given, 
according to equation (24). 


The last equation implies that the effect of the locked dislocation is 
dominant at points much closer to it than is the nearest free dislocation. 
It is shown in the Appendix, (§ A3) that 


t= LOD OV: 2 | ay 


and also TAD it To. 


These relations enable us to set bounds to the departure of 7/7) from the 
simple expression 1/(L/z) for a given a; see figs. 1 and 2. 
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For large x we have 


TT +NA/e% or t/tp>=1+3Ljx, e<L, . . . (23) 
as can be seen from the inequality 
nA A A nA 


Pea al alae 


which follows from the facts that 0<]2;[<L. Equation (23) shows that at 
large distances the stress is the sum of the applied stress and that due to all 
the dislocations regarded as congregated at the origin. 


Fig. 2. 


1-25 


0 0-04 0:08 afl 012 0-14 


Shear stress in front of a row of ten held-up dislocations. a, 6 upper and lower 
bounds from equations (22). c, d upper and lower bounds from 
equation (a9) (Appendix, §A.3). ©, accurate values calculated from 
the polynomial (15). The stress falls off monotonically with x: the 
convexity arises from plotting 7/79/(L/x) instead of 7/79. 


Zener (1948) has compared a slip-band both to a crack and to an array of 
dislocations and it is interesting to compare the results of this section with 
those of Starr (1928) for a crack. He finds the following expressions for an 
infinitely narrow two-dimensional crack of length / extending along the 
x-axis from z—0 to z=—l when the body containing it is subjected to a 
uniform applied shear stress p,,—=7 . The stress in the plane y=0 beyond 
the tip of the crack is given by the exact formula 


T ca eee l a le t ‘6 
on TN (ze) (7 sete) ANG (24) 
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If the crack is not infinitely narrow but has a radius of curvature p at the 
end (24) is only correct if z>p. 7/7) then reaches a maximum 1n the 
neighbourhood of «=p and falls to zero at ~=0. We have 


afto=l+ dle), fel ee ee 


far from the crack. 

The discontinuity in the displacement in crossing the x-axis, 7. e. the 
relative displacement of the two faces of the crack, is given by the exact 
expression 


A= 2 Vid—l2pl2b, —t<«<0 
reducing to 


bt 
— ef & “Gi > S26 
4=—vil|z)), x<0, |x|<l (26) 


if we do not go too far from the tip of the crack. 

The elastic constants have been expressed in terms of b and A (¢f. 
equation (1)) to aid comparison, though of course the problem of the crack 
does not involve the Burgers vector. A is to be given its value for edge 
dislocations, bu/27(1—c). It can be shown that the last four equations 
also apply to a crack in an anisotropic medium if we give A the value 
bK/2z7. They also cover the case of a crack in an isotropic body in a state 
of anti-plane strain with 7) and 7 representing the yz stress component 
provided we put A=bu/27, the value appropriate to screw dislocations. 
Consequently these equations can be used to compare a crack under plane 
strain or anti-plane strain with an array of edge or screw dislocations 
respectively. 

Comparison of (22) and (24) or of (20) and (26) shows that over a certain 
range of x the row of dislocations of length L simulates a crack of length 
4L (cf. fig. 1). 

It will be noted that near the origin the stress produced by the 
dislocations and the crack with a finite curvature at the tip diverge in 
opposite directions from the value t)\/(L/2). 

From (23) and (25) it will be seen that for the dislocations 7—7, falls off 
as 2! at large distances, but for the crack as a”, that is, in the same way 
as for a group of positive dislocations near a group of an equal number of 
negative ones. It is also clear that the array of dislocations simulates only 
one end of the crack. These facts suggest that a better imitation of a 
crack would be given by the following sequence of dislocations : a locked 
positive dislocation, n free positive dislocations, n free negative dislocations, 
a locked negative dislocation. The free dislocations are prevented from 
coalescing by the applied stress and the locked dislocations prevent them 
from spreading and define the length of the equivalent crack (cf. Zener 
1948, fig. 4). The method of the present paper is not capable of finding 
their equilibrium positions. 
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§5. 

The polynomial F also has a physical interpretation for complex values 
of its argument. Let Z=ax+iy. Then for a set of screw dislocations the 
displacement w and stresses at the point (x, y) are given in the isotropic 
case by 


: b 
o+iw= 5, 08 F(Z) , 


: pb F(Z) 
Pye TP x= In F(Z) : 


¢ is the electrostatic potential in the corresponding Stieltjes line-charge 
problem. These results follow from the fact that for a single screw 
dislocation at the origin w=(b/27) tan“ (y/a). 

For an array of edge dislocations the results are less simple since problems 
in plane stress or strain cannot be solved in terms of a single complex 
variable. We find for the isotropic case 


(Dee +p) +2ipa/( 1—o)=— 2iG! 
Digi 44(Daw Dy =o iy" 


where 

G(Z)=A log F(Z), A=pb/27(1—c) 
and w is the rotation. The results follow from the value y=—Ay log r 
of the Airy function for an edge dislocation at the origin (Koehler 1941). 
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APPLE EN DIX. 


§A.1. 
The reduced equation for the problem of §4 is 


0+ (3-5) =0. Ree ee rte (Al) 


Let X=2,,_1 be the greatest root ofv. We can arrange that v>0 for x Se, 
Choose a constant €>X. We can write (al) in the form 
Dale AUG (Bie we lore ymke ve se 8 (a2) 


where 


362 J. D. Eshelby, F. C. Frank and F. R. N. Nabarro on the 


The solution of the inhomogeneous equation (a2) which satisfies v(X)=0 is 
v(x) =A sin A (w2—X)+A7 uli sin A (X—2’) g (x’) dx 
where A is an arbitrary constant. Hence 
X-+a/A ure 
o(X-fa[A=r7 | sin A (a’—X) g (x") dx 
ax 


Since the left-hand side is positive g(x’) must be positive for at least part 
of the range of integration and so 


é<X+n/A. 


So far is arbitrary. We now try to choose it to give as precise information 
as possible about X, bearing in mind that Ais a function of €. Put 


E=4n/(1+¢"), A=2¢. 
Then 


4n 2a 
XT gs a wits © 56th 0a ee ee 


and a fortiori 
X >4n(1—d?—7w/2nd). . . . « «© « © (a4) 
The right-hand side of (a4) has a minimum when ¢=(7/4n)"8 so that 
X >4n[1 —3(7/4)28 n-28], 


A slightly more precise lower limit could be found by minimizing the right- 
hand side of (a3). 


§A.2. 


The same method applied to the Hermite polynomials occurring in 
§ 3 (iii.) gives for the greatest root X 


X>v/(2n+1) [1 — (27) (2n-+1)-94], 


§ A.3. 
With the sign of x changed (a1) becomes 


n 
vy” — iat a) e e e é e e (a5) 


which has the exact solution 


v=(n—1) « L, (—«x) exp $x" , uh ae « (a6) 


with a convenient choice of the arbitrary multiplying factor. If the 


right-hand side of (a5) were zero it would have the two independent 
solutions 


%=y 1, (y), =yK, (y), where Yy=vV/ (Anz). 


ST 
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It is easily verified that v and v, agree in value and first derivative at the 
origin. Solving (a5) as an inhomogeneous equation in the usual way, 
making use of the relation 


Di ets Wg a= lies cote oss «ew (&7) 
we have 


v=y L(y) +y (lily) 4% —Ky(y)F}/8n . 6 1. -. (28) 
with the abbreviations 


I= jy Ly’) o(2') da’, H= [7K v(x’) da’, y'=»/(4nz’). 


The term y{ }/8n vanishes with its first derivative at the origin so that 
solving the integral equation (a8) would give the exact solution (a6). 
Differentiating (a8) we have 
| vy’ =2n I+ (Ip % +KyZ). 
If we write (a8) in the form 
v=yl,[1—y (1, +Ky9)/8n0] 

we find, using (a7) again, 

v=. 2nIy rhe a 

ae asa.) 
By finding upper and lower limits for.4/v we can give upper and lower 


limits for v'/v. The simplest approximation is the following. Since v(«) is 
a monotonically increasing function of x 


0<sjv<| y’ I, (y’') da’=y* 1,(y)/2n 
0 
and so, using the relation y IL=y Ip—2,, 
Ty hb y? y 
Le tne ae 1+ Té@8 pg as Gas cree (a9) 
This estimate could, of course, be improved. Using the recurrence 


relations for Bessel functions and the fact that, if m<n, I,,/I,, decreases 
monotonically to unity as y increases we find the following inequalities : 
IG Fe Dae 
oa ly = Saw Qe “tel. 
I, Gee lie Y 
With the help of these we can derive from (a9) this weaker relation free 
from Bessel functions : 


9 


2 y v ee Se 
max (. 1) SOE ala ie 1672" 


The notation on the left implies that we are to choose the greater of 2/y 
and 1. Equation (22) in the text follows from the fact that 


y=4nv/(x/L), yv' 2nv=7V/4|T)1/L 


in ordinary units. 
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SUMMARY. 


The behaviour of plasticine when permanently deformed is compared 
with that of an ideal metal, and the necessary conditions for similarity 
between their flow patterns are examined for flow in two and three 
dimensions. A brief account of the mechanical properties of plasticine 
includes a description of stress-strain curves obtained by compressing 
cylinders between lubricated plates at constant rates of loading. Plane 
strain experiments, using a plasticine model marked with a grid on a 
cross-section parallel to the plane of flow, are described, and photographs 
of the deformed grids are compared with theory. The processes thus 
examined were: extrusion through square and wedge-shaped dies ; 
indentation of a plane surface by a single wedge and by a row of wedges ; 
and compression of a narrow wedge by a flat die. The profile of a deep 
narrow block after a small indentation by a flat punch is also compared 
with theory. All these plane strain plasticine experiments agree well with 
theory. One experiment using a metal was performed—the compression 
of a narrow copper wedge by a smooth flat die. This copper wedge deformed 
in a similar way to the plasticine wedge. 


§1. [INrRoDUCTION. 


Many metal working processes can be simulated by plasticine model 
experiments. The models are cheap and easy to make, and only a small 
amount of apparatusis needed. Inthe present paper an attempt is made to 
examine the mechanical properties that plasticine (or, indeed, any 
substance) should possess, in order that its manner of deformation should 
approximate to that of an ideal metal, for flow in two and three 
dimensions. That plasticine does in fact deform in a similar manner to 
an ideal metal, under conditions of plane strain, is shown by the close 
agreement between the model experiments described below and calcula- 
tions made on the basis of existing theory for an ideal plastic—rigid 
non-hardening metal. There is no evidence, as yet, that this agreement 
will extend to flow in three dimensions. However, a model experiment 
should give a qualitative idea of the flow, provided the hydrostatic 
component of stress is compressive. 


* Communicated by R. Hill. 
SER. 7, VOL. 42, NO. 327.—APRIL 1951 2 
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§2. PROPERTIES OF PLASTICINE. 


Plasticine is a mixture of an inert mineral oil, probably petroleum 
jelly, and a mineral filler, mainly calcite. Sometimes an organic dye is 
added. Its compressive stress-strain curves at three different rates of 
loading were found by compressing cylindrical specimens, }-in high and 
l-in. diameter, between the parallel plates of a beam test apparatus. 
During the course of compression the sides of the specimens became 
slightly concave. This effect was probably due to extrusion of the 
grease lubricant between the ends of the specimens and the compression 
plates. In all the stress-strain curves there is a sharp bend, after which 
the material deforms at an approximately constant stress. Unfortunately, 
the curves cannot be directly compared, because the temperatures at 
which the tests were carried out differed by 1-2° C. and the yield stress 
is very temperature dependent. A rise in temperature from 17°C. to 


Fig. 1. 


0 0, 0-2 03 0-4 


0:5 
sTRAIN (5% 
Stress-strain curve for plasticine under uniaxial compression. 

Rate of loading 21-15 lbs.min.-1. Temp. 21° C. 


21° C. caused a decrease in yield stress of about 30 per cent. The curve 
for the rate of loading 21-15 Ib./min.~+ at 21° C. is shown in fig. 1. The 
curves for the two lower rates are similar, except that for strains greater 
than about 0-1, the small constant slope of the curve is greater. For 
the rates 21-15, 15-66, and 11-93 lb./min.~, the slopes are, approximately, 
3-0, 4:0 and 4:5 Ib./sq. in.-1 respectively. The rates of loading were 
chosen to produce rates of strain comparable to those which occurred in 
the model experiments. If, at any stage of a compression test, loading 
is stopped and the load is kept constant, deformation continues at a 
rapidly decreasing rate becoming negligible after a few minutes. For 
example, at the end of a test where the rate of loading was 21-15 lb.min.~1, 
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the load was kept constant, and after one minute the strain had increased 
from 0-51 to 0-53; after a further ten minutes it was 0-54 and over the 
next sixteen hours it increased to 0:55. Provided air bubbles are removed. 
by working the material, its permanent compressibility is negligible. 
This was checked by weighing cylindrical specimens in water before and 
after 50 per cent compression. The change in volume was always less 
than 0-2 per cent. 

There appears to be no generally accepted explanation of the 
stress-strain properties of plasticine. It is possible that friction between 
the mineral particles embedded in the petroleum jelly prevents 
continuous flow until a critical shear stress is reached. There may be, 
in addition, electrostatic forces between the particles, as has been 
suggested for clay and paints (Houwink 1937). 


§3. COMPARISON OF PLASTICINE AND METAL. 


It might appear surprising that close comparison is possible between 
the deformation of plasticine and that of a metal. It will now be shown 
why this comparison is possible under conditions of plane strain. 
Consider plastic flow in plane strain, where plastic strains are large, and 
the freedom of movement of the material in plastic flow is not constrained 
completely by material still elastic. Under these conditions of pure 
shear, plasticine behaves like a non-hardening isotropic metal, provided 
the flow is slow, so that it can be treated as quasi-static. This is because 
any element of either material deforms in shear under an approximately 
constant stress; both are virtually incompressible ; and their principal 
axes of stress and strain-increment coincide. The last property must be 
true of any isotropic body for which the ratios of the plastic strain- 
increment components depend only on the ratios of the current stress 
components, and not on the previous stress or strain history. This can 
be expected on physical grounds, if, after deformation, the new random 
arrangement of the mineral particles is indistinguishable from the original 
one, so that there is no ‘memory ”’ of past deformations. This appears 
to be true for plasticine judging by the agreement between theory and 
experiment. The precise degree of similarity between the plasticine and 
the metal flow patterns depends, of course, on how closely the shapes of 
the respective shear stress-strain curves correspond. 

The equations governing the plane plastic strain of an ideal non- 
hardening metal were originally derived by Hencky (1923), (stress 
equations), and Geiringer (1930), (velocity equations). The theoretical 
flow patterns given in this paper were calculated on the basis of this 
theory. The flow patterns for plasticine can be expected to differ in the 
region of small strains from the calculated patterns, because its yield 
point is not quite sharp, nor, of course, is it rigid before yielding. 

Tt is not known how closely plasticine will correspond to a non-hardening 
metal in three-dimensional plastic flow, since no experiments have yet 
been made. If good agreement is to be obtained, the yield criterion for 
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plasticine must approximate closely to that of the metal; for many 
metals, the Mises criterion is a good approximation. In addition, the 
relation between Lode’s variables for stress and strain-rate must be the 
same for the two materials. For a Newtonian viscous fluid Lode’s 
variables are equal; this is also a good approximation for many metals 


(Hill 1950). 
§4, EXPERIMENTAL DETAILS. 


Stone-coloured plasticine, containing no dye, was used for most of the 
experiments. After carefully working the plasticine by hand to eliminate 
air bubbles, the specimens were cut to shape with a fine wire. Each 
specimen was made in two halves, }-in. thick. On the inner face of one 
half a +-in. square grid was lightly stamped in ink. French chalk was 
used to prevent the two halves sticking when placed together. After 
deformation between parallel lubricated constraints, the two halves could be 
easily separated and the deformed grid photographed. If it isinconvenient 
to have the specimen in two halves, either the grid must be ruled on 
an outer face, or the model may be built up in layers, using two different 
colours of plasticine (see expt. 5 below). Unfortunately it is difficult to 
find different colours with exactly the same mechanical properties. 
Medium grease was used for the lubrication of external surfaces. 


§5. COMPARISON OF THEORY AND EXPERIMENT. 


With one exception, the following experiments were carried out under 
conditions of plane strain, using plasticine models. The exception is 
the compression of a narrow wedge of annealed copper by a smooth flat 
die (see expt. 6). 

(1) Direct extrusion through a square die with 66-7 per cent reduction 
(figs. 2a and b, Pl. XV.). 


(a) Full friction at the container walls. 

The container was lined with emery paper to ensure that the friction 
was sufficiently high to cause local yielding of the plasticine at the walls. 
The distortion of a square grid was calculated, using the slip line field 
given by Hill (1948 a) for direct extrusion. It is particularly interesting 
that the observed region of dead plasticine agrees well with that asswmed 
in the theory. 


(b) Smooth container walls and die face. 

The slip-line field is an extreme case of the one given by Hill (1948 a) 
for extrusion through a smooth wedge-shaped die (the die-angle here 
being 90°). In the theory it is assumed that there is no dead region, but 
in practice lubrication was not quite perfect, and it was impossible to 
avoid a small dead region in the corner, which modified the flow pattern. 
The cusps on the centre line in the theoretical pattern are a result of the 
tangential velocity discontinuities across the entry and exit slip lines. 
In practice, since the yield point of plasticine is not sharp, there is a region 
of rapid shear instead of a velocity discontinuity, and the cusps and sharp 
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corners in the flow lines are rounded off. The experimental and 
theoretical grids are out of step, because the actual moment of comparison 
was not the same. 

(2) Direct extrusion through a smooth wedge-shaped die; semi-angle 
of die 15°; 34 per cent reduction (fig. 3, Pl. XV.). 

The slip-line field for this problem (Hill 1948 a) is the same as that 
found for the plane strain problem of sheet drawing. The latter problem 
is dealt with in a paper by Hill and Tupper (1948), and the theoretical 
distortion of a square grid is taken from that paper. The theoretical 
region of uniform motion is well demonstrated in the experimental 
pattern. The squares become parallelograms, in this special configuration, 
due to a single line of shear at entry ; the preservation of the parallelo- 
grams indicates uniform flow. The boundary of the region of large plastic 


(a) Diagram of the experimental arrangement for wedge indentation and coining. 
(b) Segment of the coining process. 


strains, defined by the sharp bends in the longitudinal flow lines, agrees 
well with theory. There is again a rounding off of cusps and sharp 
corners. The slope of the straight sections of the deformed transverse 
lines is greater than in theory, probably due to friction at the die face. 

(3) Indentation of a plane surface by a smooth wedge of semi-angle 50° 
(fig. 4, Pl. XV.). 

The distortion of a square grid was calculated according to the theory 
given by Hill, Lee and Tupper (1947), who computed one for a wedge of 
semi-angle 30°. The wider semi-angle of 50° was chosen for this experiment 
because the region of non-uniform motion is larger, and a clearer comparison 
can be made with the theory. It is difficult in practice to obtain good 
lubrication using a full wedge. Instead, therefore, a rectangular block, 
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constrained between two parallel lubricated guides was compressed 
against an oblique face (fig. 8a). This face represented half a wedge, and 
one of the guides replaced the line of symmetry through the apex of the 
wedge. It is assumed that the slip-line field and velocity solution are 
the same as for one half of the field of a full wedge, so long as there still 
remains a part of the original surface undeformed. The stresses in the 
plastically deforming region differ only by the addition of a hydrostatic 
pressure due to the lateral constraint of the guides. The oblique face 
was removed and re-lubricated intermittently, but there still remained 
an appreciable amount of friction due to newly sheared plasticine being 
continually presented to the face, both at the tip of the half wedges and 
the apex of the plasticine lip. Calculations by Hill (1950, unpublished) 
show that the effect of friction is to lessen the slope of the lip which, 
however, remains straight if the friction is uniform. By comparing 
measurements on the experimental block with the theory, the mean 
coefficient of friction was deduced to be approximately 0-3. This high 
value appears reasonable in view of the breakdown of lubrication which 
occurs for the reason already explained. The curvature of the surface 
near the apex of the lip is, no doubt, partly due to uneven friction, and 
partly to strain hardening. 

(4) Indentation of a plane surface by a row of smooth wedges of 
semi-angle 45° (fig. 5, Pl. XVI.). 


EXPLANATION OF THE PLATES. 


PLATE XV. 
Fig. 2. 
Direct extrusion of plasticine through a square die with 66-7 per cent 
reduction. (a) Perfect friction ; (6) Smooth. 


Fig. 3. 
Direct extrusion of plasticine through a smooth wedge-shaped die 
with 34 per cent reduction ; total die angle 30°. 
Fig. 4. 
Indentation of a plane plasticine surface by a smooth 100° wedge. 


Pratt XVI. 
Fig. 5. 
Three stages in a plasticine coining process using smooth 90° 
wedges : Ratios of depths of indentation to width of specimen. 
(a) 0-3, (b) 0-42, (c) 0-5. 
Fig. 6. 
Compression of a 40° plasticine wedge by a smooth flat die. 
Fig. 7. 
Compression of a 34° annealed copper wedge by a smooth flat die. 
(a) Vertical section. 
(6) View from above showing compressed surface. 


Phil. Mag. Ser. 7. Vol. 42. Pl. XV. 


A. P. GREEN 


Ni 


Fig. 3. 


Fig. 4. 


Ee (SHS ER 


a 
call | ) 
4 
oy T 
Bees 
ae S| 
$ —t-— 
pr 


To face page 370 


Phil. Mag. Ser. 7. Vol. 42. Pl. XVI. 


A. P. GREEN 


Models to Simulate the Plastic Flow of Metals 371 


This may be regarded as a process of coining. A block is constrained 
on all faces, except the surface which is to be indented. The flow pattern 
can be split up into exactly similar segments between the planes of 
symmetry through the angles of the wedges (fig. 8b). Each of these 
segments is identical with that for half-wedge indentation, described in 
the previous section, and the experiment can, therefore, be carried out in 
a similar manner. In this case, however, compression was continued, 
until the specimen was completely “coined”. After the region of 
deformation has reached the line of symmetry between adjacent wedges, 
(2. €. the outer guide), the configuration is no longer geometrically similar. 
The slip-line field is then continually changing and no theoretical solution 
has yet been found. It is interesting to see that there is little deformation 
below the tip of the wedge. 

(5) Compression of a 40° wedge by a smooth flat die (tis. "6, PI eV Jey 

The model was made of layers of two different colours. The mode of 
deformation shown in fig. 6 (Pl. XVI.) is one of the ways in which a narrow 
wedge can deform. The observed deformation suggested the following theory. 
It appears that a simple shear takes place along a single moving slip-line, 
BC, which makes an angle of 45° with the smooth die. Any element of 
the wedge remains stationary until reached by the moving slip-line. 
It then shears instantaneously, and continues to move as part of the 
rigid body ABC, parallel to BC. If angle ACB were less than 45°, the 
material in ABC would be unable to support the calculated forces on BC. 
This restricts the range of validity of the solution to wedge angles less 
than 45°. The experimental pattern agrees well with this theory. The 
bends in the deformed layers are sharp and lie on a straight line which 
makes an angle of 43° with the die face. 

Other experiments were attempted with wider angle wedges but it 
was difficult to obtain consistent results because of varying friction 
conditions at the die face. Two possible symmetrical solutions have 
been suggested by Hill (1948 b). In one of these a steadily increasing 
wedge-shaped region of material is carried down with the die. This type 
of solution tends to occur when lubrication is poor, since it involves no 
relative motion between the die and the plasticine surface in contact 
with it. Itis valid for wedge angles greater than 28° 4’. When lubrication 
is good the other type of solution, which is valid for wedge-angles greater 
than 53° 12’, is sometimes observed. The wedge gradually divides down 
the centre ; newly sheared plasticine is continually presented to the die 
face from the centre and is then displaced outwards towards the sides. 
Sometimes a wedge will change its mode of deformation as frictional 
conditions alter. 

(6) Compression of a 34° annealed copper wedge by a smooth flat die 
(figs. 7a and b, Pl. XVI). 

The wedge was not constrained at its ends, but towards the centre 
conditions approximated to plane strain, since the ratios of depth of 
compression to the dimensions of the wedge were small. Initially, the 


372 A. P. Green on the Use of Plasticine 


wedge was 0-75 in. high and 1-36 in. long. It was compressed through 
0-1 in., in six stages, between polished plates made of sintered hard 
metal, which were lubricated intermittently with calcium oleate. There 
was no deformation at the base, but the length of the top of the wedge 
increased by 0-06 in. After compression the wedge was sectioned 
vertically (i.e. along a plane perpendicular to its length), 0-35 in. from 
one end (fig. 7a, Pl. XVI.). It can be seen that its mode of deformation was 
similar to that shown in fig. 6 Although work-hardening broadened the 
single slip line, BC, into a region of rapid shear, the final shape of the 
section, except for the rounded corners, agrees well with the theory for 
an ideal metal. Initially, some vertical sections of the wedge commenced 
to shear to one side and some to the other, while some sections even 
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Indentation of a narrow plasticine block by a smooth flat punch. 
(a) Diagram showing nature of initial deformation (exaggerated) ; 
(b) Profile of block after indentation. 


deformed symmetrically. Viewed from above, the compressed surface 
zig-zagged along its length. As compression continued, the longer 
unsymmetrical sections began to carry the other sections with them. 
Finally, half the wedge was shearing to one side and half to the other 
(fg. 7 Oe elo Vee 

(7) Indentation of a deep narrow block by a smooth flat punch*. 


ae Sc a a nial ha 
*'The theory has been worked out by R. Hill (1950) and I quote, by kind 
permission, from his as yet unpublished results. , 
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With a sufficiently narrow block, the plastic region spreads out to the 
sides as the load is increased, and the initial deformation of the corners 
is a rotation outwards (fig. 9a). The rotation on each side takes place 
along a curved slip line which meets the side of the block at a point, A, 
whose position can be calculated. A plasticine block, ?-in. wide, was 
indented to a depth of 4-in. by a flat lubricated punch }-in. wide. After 
indentation the shapes of the lateral faces of the block were plotted with 
a travelling microscope. The occurrence of a rigid-body rotation and the 
position of A agreed well with theory (fig. 9 6). 


§6. CONCLUSION. 


Under conditions of plane strain, agreement with theory has been shown 
to be encouragingly good. It is planned to extend these experiments to 
other processes, and also to measure the applied loads. This should 
provide an additional test of the similarity between plasticine and an 
ideal metal. Similar experiments using a metal such as pre-strained 
copper are desirable for further comparison with theory and with 
plasticine. An experiment in which lead blocks were indented by steel 
wedges coated with vaseline is described in the paper on wedge indentation 
(Hill, Lee and Tupper 1947). When the theory is not known, a plasticine 
experiment may assist in finding a solution (for example, see expt. 5). 
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SUMMARY. 


Calculations similar to those of Stoner (1938, 1939) are carried out for an 
energy band for which the energy density of states is constant. These 
have a bearing on calculations of the actual band form for nickel (Fletcher 
and Wohlfarth 1951). After brief reference to previous work and a state- 
ment of the basic assumptions of the treatment, relations are derived 
giving the dependence on temperature and the interchange interaction 
energy of the spontaneous magnetization below and the susceptibility 
above the Curie point. Various explicit limiting relations are deduced. 
The numerical results of the calculations are given in Table I. and fig. 1 (§2). 
Formule are also derived which give the internal energy and specific heat 
over the whole temperature range. Relations are given for the most 
important experimentally observable thermal characteristics, such as the 
specific heat discontinuity at the Curie point. The results are shown in 
figs. 2 and 3 (§3). Experimental results which are briefly discussed (§4) 
include those on the magnetic properties of nickel and its alloys and the 
thermal properties of nickel, including the energy and specific heat 
characteristics, the temperature variation of specific heat and the magneto- 
caloric effect. The high temperature specific heat of paramagnetics is also 
discussed. In an Appendix certain results for an energy band of more 
general form are given. 


$1. INTRODUCTION. 


THE collective electron treatment of ferromagnetism was originally 
developed by Stoner (1938 a, 1939, 1948, 1951) for an energy band of 
standard form, for which the energy density of states is proportional to the 
square root of the energy. In the present paper similar calculations are 
carried out for an energy band for which the density of states is constant. 
Consideration of this problem arose mainly as a result of calculations, 
which are still proceeding, aiming at a determination of the actual band 
form for nickel. Preliminary results of these calculations (Fletcher and 
Wohlfarth 1951) indicate that the actual band form near the top of the band 
is almost certainly intermediate between that assumed by Stoner and that 
considered here. A comparison of the two sets of theoretical results 
serves to show to what extent discrepancies between theory and experiment 
are due to the assumption of a particular band form and to what extent 
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they result from other assumptions of the treatment. As compared with 
the parabolic band, the band with constant density of states has the 
practical advantage that the main formule may be obtained in closed 
algebraical form; for the parabolic band extensive numerical computation 
is often required. 

Rectangular bands have been considered previously by Néel (1938) and 

Watanabe (1948). Néel obtains expressions for spontaneous magnetization 
and susceptibility equivalent to relations (2.7) and (2.9) of the present paper. 
Watanabe considers susceptibility in some detail, but attempts to represent 
his results by a Curie-Weiss law, which has no immediate significance in 
the treatment. This type of representation underlies his analysis of the 
experimental results for iron, cobalt and nickel, and little significance can 
be attached to the estimated values of the specifying parameters for these 
metals. Watanabe also discusses the experimental results for some nickel 
alloys, but does not seem to have appreciated the primary effect of alloying 
on the magnetic properties (cf. §4). It is believed, however, that 
Watanabe’s procedure in attributing some of the details of the observed 
susceptibility curves to band overlap effects is essentially correct. In the 
same paper the temperature dependence of specific heat is investigated 
and the calculated curve for nickel is similar in shape to that given here 
(fig. 2). 
- In the present paper calculations of spontaneous magnetization and 
susceptibility, appropriate to the basic assumptions (§2), are followed by 
calculations of the internal energy and specific heat ($3). Use is made of 
the recently published table of the Fermi—Dirac function of order 1 
(Rhodes 1950). The bearing of the theoretical results on experimental 
data is discussed in §4. In an Appendix certain results are derived for 
an energy band of more general form. 


§2. Magnetic PROPERTIES. 


The calculations are based on the following assumptions: (i.) The energy 
band contains N particles, in actual ferromagnetics holes in the 3d band. 
At absolute zero and in the absence of any interchange interaction the 
holes occupy the band in balanced pairs over an energy range between the 
top of the band, «,, and the top of the Fermi distribution, «,—¢p. The 
energy density of states for particles of either spin, v(e), is constant and is 
given by 

v(e)=N/2e . ie) ees es eon?) 


(ii.) The interchange interaction effects give rise to a term in the total 


energy of amount 
LOSE ee (22) 


where ¢ is the relative magnetization, related to the magnetic moment, M, 
by (=M/Nup, and 0’ is a constant parameter. 
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Fundamental relations. With the above assumptions a development of 
the statistical treatment (Stoner 1938a) leads to the following expressions 
for Cas a function of temperature, T, and applied field, H: 


t=3 a {Fo(n+B-+P’)—Fo(n—B—B')}; 
a birch cep gei Ge} 
i=) (Rint B+ 84+ Falr—B-P)}. | 


Here 7 is related to the chemical potential, €, by y=€/kT and B=k6'C/kT, 
f’=pH/kT. The function F, is the FermiDirac integral of zero order 
which may be evaluated exactly, 
is dx 
F.(n)= | Cree eras 
07) o exp (c—n)+1 . 

In the absence of an applied field manipulation of (2.3) and (2.4) gives an 
implicit relation for spontaneous magnetization in the form 


ke’ or) fexp((Q+0/rI-1 
@ THE expt} vf fale seta eh 


where t=kT/ep. Two useful alternative relations may be obtained: 


In\(T ety ee ce ey 


= 4=F tanh {tanh (¢/27)/tanh (1/27) }, Sieh i) 
1 Fsinh-! fexp (—1/r) sinh B}, St) TURE (207) 


where B=(k6'/<9) ¢/7. The Curie temperature, 6, is given by 


eo tec 
hare (1 ) 
or inversely «age ple tae 
S. € 
gg exp (en) | 

The inverse of the susceptibility above the Curie point, y, is given by (2.3) 
and (2.4) as 

pris Ps ON se aa 

wee orig Toe ee ee (—1/7)] —kO’/éo, . . ° (2.9) 
if wH/kT and ¢ are small compared with 1. 

Limiting relations. (i.) The limit k@/eg> oo. Relation (2.8) shows that in 
this limit 0’-@ and hence, from (2.6), ¢=tanh (C6/T), the classical Weiss 
relation for spontaneous magnetization. Relation (2.9) shows that in the 
same limit y=Ny?/k(T—8@), the Curie-Weiss law. 

(ii.) Condition for ferromagnetism. Relative magnetization at absolute zero. 
Relation (2.8) shows the condition for ferromagnetism (finite spontaneous 
magnetization at absolute zero) to be k0"/e,>1 and (2.5) shows that as soon 
as ferromagnetism occurs at all the relative magnetization at absolute zero, 
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€o, is unity. As shown in the Appendix (relation (A.8)), however, in a band 
of general form there is usually a range of values of k6’/e) over which 
although ¢, is finite, it may be less than 1. For parabolic bands, Py 
example, this range was shown by Stoner (1938 a, relation (5.10)) to be 
3 <k0'/e) <2-1/8, the lower limit giving the condition for ferromagnetism in 
that case. In this respect, therefore, the rectangular band form with 
strictly constant density of states is a very special case. 

(iul.) Approach to saturation, kT/e)>0. In this limit €+¢,=1 and relation 
(2.5) gives the law of approach to saturation as 


eet ai Qe, /k0" 
t=1——exp{—Ft(——1)}, e 6 . 6 (2.10) 


analogous to Stoner’s relation (1938a, relation (5.21)) for parabolic bands in 
the range of values of k6’/e, for which {;=1. For parabolic bands, how- 
ever, the law of approach takes the form (¢/f))?=1—«(£,) (T/@)? when 
Co <1, where « is in general close to 1 (Stoner 1938 a, relation (5.26)) 

(iv.) Magnetization near the Curie point, T/6>1. Relation (2.6) may be 
used to show that in the neighbourhood of the Curie point 


(2=6h6/«, tanh (€/2k0)(1—T/0), 
which, on transforming (2.8) into coth (¢)/2k0’)=2k6'/e)—1, gives 


61hO/€5 


ryan rire ae eee. (2,11) 


Relation (2.11) gives correctly the classical dependence of magnetization on 
temperature near the Curie point, ¢?=3(1—T/@) for k@/e>> 0, but shows 
that the slope of the ¢?, T/ curve decreases as k6/e, decreases. In the limit 
kO/e>>0, kO’/eg>1, (2.11) becomes (?=6k6/e, (1—T/O). 

Computations. Relation (2.5) was used to compute values of ¢ for 
T/0=0-0 (0-1) 0-9 (0-05) 1-0 and for 46’/e,=1-01, 1-04, and 1-10, using a 
balancing method of solution. Since (.=1, ¢ equals the usual reduced 
specific magnetization, c/o). The range of values of £6’/ey covered is 
indicated by the susceptibility results to include practically all ferro- 
magnetic metals and alloys. The susceptibility above the Curie point was 
calculated using relation (2.9) for T/9=1-0 (0-2) 3-0 and the above 3 
values of £6'/e9, in the form ) 

lpHe 1pH | pool 

Gey kD CKO kB x’ 
which gives the classical curve as the straight line (u9/k0) 1/x=T/O—1. 
The final results are given in Table I. and fig. 1. 

Fig. 1 also includes the values calculated by Stoner (1938 a) for parabolic 
bands with £6’/ep5=2-1/8. The most marked difference between the results 
for ¢ calculated for rectangular and for parabolic bands is that for the 
former the ¢, T/@ curves move downwards as k6’/e, decreases. For the 
latter (cf. Stoner 1938a, fig. 5) the curves giving ¢/¢) move downwards as 
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k0’/<, decreases to the lower limiting value (2~1/*) for ¢)=1, but they 
actually move very slightly upwards as /6’/e) decreases further to the lower 
critical value 2, mainly as a result of the simultaneous decrease of Co. 

The susceptibility curves, on the other hand, have similar characteristics 
for both types of band form. The reduced curves move downwards as 
k0/«, decreases, coinciding with the T/@ axis as k6/e,+0 (cf. also Wohlfarth 
1949a, fig. 9). 


TasLe TI. Spontaneous magnetization and susceptibility. 


k8'/eo 1-01 1-04 1-10 o-) 
kO/<o 0-2167 0-3069 0-4170 ro) 
T/é Relative magnetization. 

0-0 1-0000 1-0000 1-0000 1-0000 
0-1 0-9889 0-9976 0-9996 1-0000 
0-2 0-9553 0-9799 0-9919 0-9999 
0-3 0-9083 0-9468 0-9702 0-9974 
0-4 0-8505 0-9003 0-9341 0-9856 
0-5 0-7826 0-8410 0-8834 0-9575 
0-6 0-7041 0-7679 0-8161 0-9073 
0-7 0-6122 0-6780 0-7286 0-8286 
0-8 0-5015 0-5638 0-6123 0-7104 
0-9 ~ 0-3556 0-4058 0-4448 0-5254 
0-95 0:2517 0-2897 0-3189 0-3795 
1-0 0-0000 0-0000 0-0000 0-0000 

Reduced inverse susceptibility. 

1-2 0-0546 0-1006 0-1363 0-2000 
1-4 0-1312 0:2219 §$0-2879 0-4000 
1-6 0-2271 0-3587 0-4501 0-6000 
1:8 0-3388 0-5071 0-6199 0-8000 
2-0 0:-4638 0-6646 0-7952 1-0000 
2-2 0-5995 0-8287 0-9748 1-2000 
2-4 0-7439 0-9983 1-1577 1-4000 
2-6 0-8955 1-1560 1-3429 1-6000 
2:8 1-0533 1-34.96 1-5303 1-8000 
3-0 1-2157 1-5299 1-7193 2-0000 


’ 

T/#, reduced temperature. For T/@<1 the figures give ¢, the relative 
magnetization. For T/@>1 they give (j0/k6)(1/x), the reduced inverse mass 
susceptibility. The ¢ values for k6/eg=co are taken from Stoner (1938 a 
Table I.) ; the susceptibility values in this limit are given by ; 


(09/k6)(1/x)=T/O—1. 


§3. THERMAL PROPERTIES. 
Fundamental relations. With the basic assumptions as in §2 and Stoner’s 
nomenclature (1939) the expressions for the reduced electronic and 
‘«‘ molecular field ’’ energies in the absence of an external field become 


Ey 1/kT\2 | 
m= a(-) Fig +B)+Fi(n—B)}, 2 6 6 © (3.1) 
By 1 


Ne, 5 & ap PE oS, al Pe ist ae: (3.2) 
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and the total internal energy is given by the sum E=E,+Ey. Here 
Fy(n)= [_ xde/lexp (en) +0), 
0 


the Fermi—Dirac function of order 1 (Rhodes 1950), and 7 and B were 
defined in §2. The reduced specific heat, C/Nk, may be calculated, with 
Stoner, as the sum of three contributory terms 
C,/Nk= {0(Eg/Neq)/07}., 
C,/Nk={0(Eg/Ne,)/0¢}, de/dr, 
Cyy/Nb={d(Byg/Neq)/dl} dC/dr, 
where t=kT/e). These may be calculated directly from (3.1) and (3.2) if 


use is made of the following auxiliary relations, derived from (2.3), (2.4) 
and the general relation F/(n)=sF,_,(7) (Rhodes 1950), 


dF, 
pe =F,(y+8)=In[1+exp (n+ A)J=(146)/7- . (3.3) 


It is found that 


C 2E il 
NEA No eae Oe ee 


Coat dl 
NETH +9X,-1-OX IE, Pe 
Oy a, MU 
Nk I One Ss c dt ? 
where X ,={l—exp[—(1+¢)/7]}-1. 


The functions X may be expressed in more convenient form using (2.5) 
and (2.7), giving 
sinh nZ/r 
<= "sinh Er CoP [F (n— 1) Git eee) 


where n=k6'/e). The relations for the contributory specific heat terms 
then become 


6) lL W2k, inh T 
xi =- Ne, = ee {(1+ 2?) cosh (n—1)¢/7—2¢ sinh m—1)¢/7} | ? 


C, inh Td 
Ni = {6 cosh (n—1)e/r—sinh (n—1)E/7} ee oe 
Neon ear 


- (3.6) 

The value of y, required for the evaluation of Ey/Ne, and C,/Nk, may be 

calculated from the relation, derived from (3.3) and (3.5) 
Ral sinh n¢/7 

)— a —In ‘sinh ¢/7 > . . . . . . (3.7) 


’ 
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and the temperature derivative dl/dz from (2.7), giving 


es = g + ={(a—1) coth (een —n coth = pbielner 358) 


T 


Energy and specific heat characteristics. Inmiting values. Relations 
giving the internal energy and specific heat at low temperatures and at the 
Curie point as functions of the specifying parameter k6/e, are useful in 
briefly summarizing the thermal properties of a wide range of ferro- 
magnetics. The relations may be used immediately in an analysis of 
experimental results (§4). With Stoner (1939) calculations are made of 
the internal energy change between absolute zero and the Curie point and 
the specific heat at low temperatures and just below and just above the 
Curie point. The relations are obtainable in algebraical form by direct 
manipulation of the expressions derived above. In the derivation use is 
made of the limiting relations for F,(y), viz. 


Fi()=37 +477—O(E"), 4>~, 
Fy (n)=e"—O(e?”), 1>— 2. 
(i.) Internal energy change between 0°K and the Curie point. 
It follows from (3.1), (3.2) and (3.7) that 
E,—E,_ ké Eo k@’ eo (,__ 1 ko" 
(ii.) Specific heat in the limit KT/«,—> 0. 
In this limit C, and Cy, tend to zero since d¢/dr->0 (cf. (2.10)), so that 
C=C,. Also from (2.10), (3.3) and (3.7) 


2 2 (k0" 


€0 
where t=kT/e). Hence from (3.1) 


ee ae 
Ne, Vee utc 
Oh oie 
so that Nie ohiee 
C/Nk yap 0 
: Sees St i ee eee I) 
Peis ( T/0 ieee 67 Ep ( ) 


(iii.) Specific heat just below Curie point. 
The three contributory terms become from (2.11), (3.1), (3.6) and ((3.7) 


C, Fat ko tin) Mi 
(si), 2() Ge) 7 
CoN ae U5 kO coralieee, gamit len 
(Si), 89 amet aT 


One = kO"/e€o 
NkJoo | 240'/eg—1° 
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(3.11) 
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(iv.) Specific heat just above Curie point. Specific heat discontinuity. 
It is easily seen that (C/Nk))49=(C,/Nk)o_o, given in (3.11), and hence the 
specific heat discontinuity at the Curie point becomes 


30’ k6’/es—1 
ACINK=(C/Nk)o-o—(C/NEoso= F 3EB7e—1° a Ee) 5 


(v.) Temperature variation of specific heat of paramagnetics. 

An immediate application of the relation for (C/Nzk),_o in (3.11) is to the 
calculation of the electronic specific heat of paramagnetics. This may be 
derived from, (3.11) by replacing k6/e, by t=kT/e) and consequently 
k6’/e, by {1—exp (—1/7)}“ (cf. (2.8)). The resulting expression is 


(rE), =2rF, {In [exp (1/7)—1]}—{7[1—exp (—1/7)]}-“%. . (3.13) 


An alternative derivation of (3.13) is by direct differentiation of the 
expression for the internal energy of paramagnetics 


(we), =F, {In [exp (1/r)—1]}, 


which may be derived from (3.1) and (3.3) by putting B=0. 
Relation (3.13) shows that in the limit 7 0 


(xz) =47"7, wi ns, de ip het eck ie mean eke 
p 


twice as large as the low temperature specific heat of ferromagnetics (cf. 
(3.10)). This change arises from the fact that in a paramagnetic the states 
near the top of the Fermi distribution are doubly occupied at low tempera- 
tures, while in a ferromagnetic, with saturation always complete in the case 
of rectangular bands, they are singly occupied. For parabolic bands 
saturation may not be complete at 0°K. if 40’/e) is small enough and the 
corresponding increase of specific heat as k0’/e) decreases is more gradual 
(cf. Stoner 1939, relation (3.15); also Appendix). 


(vi:) The limit k0/¢,—> 00. 


The expression for internal energy in this limit may be derived from 
(3.1), (3.2) and (3.3), giving 


E/N=kT—1h022, 
so that 7 Beta Ra ERT 
C/Nk=1—}dejd(T/9), | 


since here 6’>6. Relations (3.15) for rectangular bands may be compared 
with the corresponding limiting relations for parabolic bands (Stoner 
1939, relations (3.1), (3.2)) 


E/N= 8kT—3k002, 
C/Nk= 2 —1d@/d(T/é). 
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The second term in the expression for C/Nk, the classical “ magnetic ”’ 
specific heat, is the same for both types of band form. The first term, 
however, the electronic specific heat, is different in the two cases. This 
dependence on band form also appears in the limiting values of the various 
energy and specific heat characteristics derived above. For rectangular 
bands, in the limit £6/ Ey > ©, 


(E,—E,)/NkO=1-5, (C,/Nk)o_o=(C/N&)o.9=(C/Nk), =1, 
(C./Nk)o_9=9, (Cy/Nk)p_o=1'5, (C/Nk)p_9=2°5, (AC/Nk)=15; 
the corresponding limiting values for parabolic bands are 2, 1-5, 0, 1-5, 
3, 1-5. A similar dependence on band form is immediately apparent in the 


respective expressions for the electronic specific heat at low temperatures 
(cf. (3.10), (3.14) and Stoner 1939, relation (4.15); also Appendix). 


Fig. 2. 
30 
a © Grew (smoothed values) 


0-0 0:2 0-4 06 TH 08 [0 2 


Temperature variation of specific heat. 


C/Nk, reduced specific heat, T/#, reduced temperature. Calculated values are 
for k6’/eg=1-10. Experimental values are for nickel, due to Grew (1934; 
of. also Stoner 1936, and § 4). 


Computations. The internal energy and specific heat expressions 

_ derived above were evaluated numerically for 'T/@=0-0 (0-1) 0-9 (0-05) 1-0 

and for k6’/ej=1-10, using the corresponding values of ¢ given in Table J 

The particular value of &0’/e, chosen is indicated by the susceptibility 

results to be close to the actual value for nickel (cf. §4). The calculated 

variation of C/Nk with T/@ is shown in fig. 2, together with experimental 
values for nickel discussed in § 4. 


2D2 
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The curve may be compared with Stoner’s (1939, fig. 2), from which it 
differs only slightly in shape. Watanabe (1948) also gives a calculated 
specific heat curve for nickel, using a value of k6"/<) larger than the one of 
fig. 2 owing to an incorrect analysis of the susceptibility results (cf. §4). 
He gives, moreover, no indication in his paper that values of Cas dependent 
on T/0, which are clearly necessary, have been calculated. Watanabe's 
final curve, however obtained, agrees quite closely in shape with Stoner’s 
and the present curve. 


Fig. 3. 


0 
0-0 0:5 KOK, +0 5 
Energy and specific heat characteristics. 

1, (E>—E,)/Nk@, reduced internal energy change between absolute zero and the 
Curie point, 0; 2, {(C/Nk)/(T/0)}pjo,9, ratio of reduced specific heat to 
reduced temperature in the limit T/@+0; 3, (C/Nk)g», reduced specific 
heat just below the Curie point; 4, (C/Nk)g:9, reduced specific heat just 
above the Curie point ; 5, 4C/Nk, reduced specific heat discontinuity at 
the Curie point. For limiting values as k0/e, > 00, see text. 


Computations have also been made of the various energy and specific 
heat characteristics given by the relations (3.9), (3.10), (3.11), (3.12) which 
were evaluated numerically for £0/<e97=0-0 (0-05) 0-5 (0-1) 1-5. The results 
are shown in fig. 3. 

The curves may be compared with those derived by Stoner (1939, figs. 
3, 4 and 5) for parabolic bands. The most marked differences appear in 
the limiting values as k6/e9—> 00, already noted, and in the low temperature 
specific heat curve (l.c. fig. 5, curve 2). 
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§4. Discussion OF SOME EXPERIMENTAL RESULTS. 


In previous papers (for summaries see Stoner 1948, 1951) results of 
experimental investigations into the magnetic and thermal behaviour of 
ferromagnetic metals and alloys were fully discussed on the basis of the 
collective electron treatment with the assumption of a particular, parabolic, 
band form. In this section it will suffice, therefore, to consider very briefly 
certain of the experimental data in the light of the theoretical results 
obtained above with the assumption of a rectangular band form. As 
before, it is mainly nickel and its alloys which are to be considered. 


(i.) Susceptibility above the Curie point. Tt has been shown (Wohlfarth 
1949 a) that for nickel the experimental reduced susceptibility curve lies 
very close to the curve of fig. 1 for a parabolic band with k6’/e,=2-18 
if this is corrected for band overlap effects. Inspection of fig. 1 thus 
indicates that with rectangular bands the appropriate value of k6’/e, for 
nickel is about 1-10, and hence from (2.8) k6/«)==0-42. For iron and cobalt 
the analysis is, for several reasons, less conclusive, but Stoner’s discussion 
(1938 b) shows that the reduced susceptibility curves for these metals lie 
very close to that for nickel. Although it is questionable whether for 
these two metals the assumption of either a parabolic or a rectangular 
band form over the whole of the relevant energy range is justified, an 
estimate for k6/e) of about 0-4 may be accepted provisionally for all three 
metals. It may be noted, incidentally, that the value of k0/<) corre- 
sponding to k@’/e,=2-1/? for parabolic bands is also close to 0-4. Watanabe 
(1948), on the other hand, gives the following as estimates for iron, cobalt 
and. nickel respectively: k6/<,=0-45, 0-37, 0-50, basing his estimates on a 
Curie-Weiss representation of the susceptibility results. With 0/<«,=0-4 
the values of «,/k are 2-6, 3-5 and 1-6 10° °K. for iron, cobalt and nickel 
respectively, the nickel value being close to that estimated previously. 

The analysis of the susceptibility results given previously for nickel 
alloys shows that in general alloying is accompanied by a shift of the 
reduced 1/X, T curves towards the T axis as the alloying content increases. 
For both types of band form (cf. fig. 1) this decrease indicates primarily a 
decrease of k6/e, on alloying, not, as Watanabe believes, an increase in 
importance of band overlap effects (cf. also Wohlfarth 1949 a). 

It may be concluded that the analysis of the susceptibility results is 
essentially similar for both types of band form. 


(ii.) Spontaneous magnetization. It has already been shown (Wohlfarth 
1949 b) that the experimental curve relating o/c) and T/@ for nickel lies 
considerably above that calculated for parabolic bands (with kO'/eg=2-*", 
as given by the susceptibility results). The observed curve lies in fact, 
very close to the classical curve (k0’/eg >). Inspection of fig. 1 shows 
a similar discrepancy although slightly less serious, between the observed 
curve (close to the classical curve) and that calculated for a rectangular 
band (with £6’/ej=1-10). The main reason for the discrepancy is believed 
to be the temperature variation of the effective 0’ below the Curie point, 
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owing to the fundamental relation (2.2) being an oversimplified represen- 
tation of the interchange interaction energy. ‘This effect, called for brevity 
the “ 6’ effect ” below, is apparently of greater importance in determining 
the shape of magnetization, temperature curves than is the band form. 

Only two complete sets of curves for ferromagnetic alloys are available, 
for Ni-Cu and Ni-Co. These have been analyzed by Wohlfarth (1949 c). 
For Ni-Co the reduced curves were shown to coincide over the whole 
concentration range (/.c. fig. 6), in agreement with the coincidence of 
the reduced susceptibility curves (/.c. fig. 3). For Ni-Cu the reduced 
experimental magnetization curves shift towards the temperature axis as 
the copper concentration increases (/.c. fig. 7). For parabolic bands, 
as already noted, the reduced theoretical curves for different Cy values 
coincide very nearly if k0'/e,<2-1/. The interpretation of the observed 
shift for Ni-Cu was here based on the decrease of ¢) with increasing 
copper concentration and the consequent increasing accuracy of (2.2) in 
representing the interchange interaction energy. For rectangular bands 
the theoretical curves themselves shift as k0’/e) decreases (cf. fig. 1), but 
¢) remains equal to unity. The interpretation of the experimental shift 
is thus more straightforward, but it is perhaps easy to see that 
fundamentally it is similar to that already put forward for parabolic bands. 

(iiti.) Energy and specific heat characteristics. The values of the five 
characteristics given in fig. 3, with k0/e)=0-42, the estimated value for 
nickel, Stoner’s values for parabolic bands (k6’/e9=271/%) and experimental 
values for nickel due to Grew (1934) and Keesom and Clark (1935), 
analysed by Stoner (1936), are shown in Table IT. 


TaBLE II. Calculated and observed energy and specific heat 
characteristics for nickel. For symbolism see caption to fig. 3. 


E,—E, eS) ( i G AC 
Nké T/@ T/0+0 Nk 6—0 Nk 6+0 Nk 


Rectangular band 


(0 /e>=0-42) 0-67 0-69 1-52 0:87 0-65 
Parabolic band 

(k0'/eg= 2-148) 0-86 1-26 1-76 1-19 0-57 
Experimental (nickel) 0-8 0-91 2:5 1-0 1-5 


It is seen (a) that the theoretical values, for the two types of band form, 
for (E,—Ko)/Nk@ and (C/Nk),,) agree quite closely with each other and 
bracket the experimental value in either case, (b) that the theoretical values 
for {(C/Nk)/(T/@) }njo+9 differ by nearly 100 per cent but that they again 
bracket the experimental value, (c) that the theoretical values for (C/Nk)y_» 
and AC/Nk are quite close to each other, but are considerably smaller than 
the observed values. The discrepancies in (c) are doubtless due to the 
“ 0’ effect ” mentioned above, which brings about discrepancies between 
observed and calculated magnetization, temperature curves. The effect 
of assuming a particular band form is shown in Table II. to beimportant only 
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in the case of low temperature specific heat. Here, as in (a), the experi- 
mental results for nickel seem to show that close agreement between 
theory and observation should be obtainable if the assumed band form is 
intermediate between parabolic and rectangular form, bearing out the 
conclusion of Fletcher and Wohlfarth’s calculation for the metal (1951) 
It may be noted, on the other hand, that the values of {(C/Nk)/(T/) }pi6->0 
for cobalt and iron are much smaller than that for nickel, 0-49 and 0-28 
respectively, so that these metals must have a much more complicated 
band form than that considered satisfactory for nickel. 


(iv.) Temperature variation of specific heat. Grew’s results for nickel 
(1934; ef. Stoner 1936) are shown in fig. 2 together with calculated 
values for a rectangular band form (k6’/ej=1-10). It is seen that the 
agreement between the two curves is very close up to T/9=0-8, but that 
for 0-8<T/@<1 the observed values lie markedly above the theoretical 
curve. The disagreement close to the Curie point is, as already stated, 
believed to be due to the “ 6’ effect”. Although this has been shown to 
introduce discrepancies between observed and calculated magnetization, 
temperature curves over the whole T/@ range, the magnetization and its 
temperature derivative enter in the expressions for specific heat in such a 
complicated way (cf. relations (3.6)) that here the discrepancies seem to 
cancel over most of the lower temperature range. 

(v.) Magneto-caloric effect. As stated by Stoner (1951) calculations on 
the magneto-caloric effect on the basis of the collective electron treatment 
are now in progress at Leeds. It will suffice, therefore, to indicate by 
an example how such calculations could be carried out in a particularly 
simple way for a rectangular energy band. The example chosen concerns 
the magneto-caloric effect in the temperature range below the Curie point 
where the applied fie!d is negligible compared with the “molecular field”’. 
The relevant experimental results for nickel have been considered by Stoner 
(1936). 

If classical statistics were applicable the effective molecular field 
coefficient, v, could be derived immediately from the observed adiabatic 
change of temperature with magnetization, since in negligible fields (H <vI) 


»=—3(F), =z), . ° ° . . e ( sL) 


where C; is the specific heat at constant magnetization I and p the density. 
The classical relation for internal energy, E, is just that derived above 
(relation (3.15)), giving the well known relation between v and the Curie 
temperature 6, 


Nké 
= 4.2 
where (=I/I,. With Fermi—Dirac statistics applicable, however, the 
correct relation for internal energy is E=E,+Ey, given by (3.1) and (3.2). 
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Hence the apparent molecular field coefficient, derivable from the magneto- 


caloric data, becomes 
, Nk’ 11 /eK, 
aan eae T2 Fy ee 


The second term is closely related to 5s eae heat term C, given in 
(3.6) and it is readily found that 


7 
dad <r $(T/0, k6" ey), | 
+ (4.3) 
sinh a ; 
where ¢=1—{f cosh (n—1)¢/7—sinh (n —blt} ar ainh Ge’ | 


and where n=k6'/e), t=KT/€9=(kO/e9)(T/@) and ¢ is a function of k6’/e9 
and T/6 (ef. relation (2.5)). It follows from (4.2) and (4.3) that the ratio r 
of the apparent molecular field coefficient, v’, to that derived from classical 
statistics, v, is given by 


r=v'jv=(6'/0)$(T/0, bO"/e), - %.  weten (44) 


a function of temperature. If the approximate value of /6’/¢, for nickel, 
1-10, is inserted into (4.3) and (4.4) it is found that 7 increases very nearly 
linearly with T/@ from 


r= (8'/0)4(0, 16’/€y)=(€9/KkO)(K0’/egp—1) =0-24 to 
74=(6'/0)h(1, £6") <>) = (8/8) (€o/k0)(k0'/_g—1)=0-63 


as the temperature increases from 0°K to 6. Stoner’s analysis (1936) of 
the results for nickel indicates a temperature increase of the apparent 
molecular field coefficient below the Curie point of this order of magnitude, 
although close to the Curie point v’ is strongly field dependent (/.¢. p. 180). 

(vi.) High temperature specific heat of paramagnetics. Relation (3.13) 
shows that for rectangular bands the electronic specific heat of paramag- 
netics tends to the value Nk at high temperatures (kT/e,>1). The corre- 
sponding limiting value for parabolic bands is 1-5Nk. This dependence 
on band form has a bearing on the experimental finding that the high 
temperature electronic specific heat of palladium is considerably smaller 
than that of nickel above the Curie point. The observed values are 
about I and 2 cal. mol.~! deg.-! respectively, whereas the number of 
holes per atom, proportional to N, is roughly the same for both metals, 
about 0-6. This anomaly has been variously mentioned (e.g. by Mott 
1949), but does not seem to have been interpreted by an explicit con- 
sideration of differences in the band form for the two metals. It must be 
pointed out, however, that other effects, such as the anharmonicity of the 
lattice vibrations and band overlap effects. make a quantitative dis- 
cussion of high temperature electronic specific heats of much less direct 
significance than discussions of thermal and magnetic properties at lower 
temperatures. 
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Results for a more general band form. 


The discussion of the preceding sections indicates that some of the results 
of the treatment do not depend very markedly on the band form, parabolic 
or rectangular, but that others are, in fact, quite sensitive to it. In this 
Appendix a very brief discussion is given of such results as are easily 
obtainable for energy bands of more general form. 

If the energy dependence of the density of states, v(e), were only 
available in numerical form, all but a very few results would have to be 
derived by a long process of numerical integration. It is easily seen, 
however, that in the limit 46/<c)—>< the magnetic properties, spontaneous 
magnetization and susceptibility, are independent of the band form and 
are given by the relations based on classical statistics. The thermal 
properties, however, depend on the band form even in the classical limit, 
as shown in § 8, and it does not seem possible to derive any general 
results in this case. For a band for which the density of states is given 
by the simple relation 


v(e)=ae", Pee ere a «ers (A 1) 


the expressions for internal energy and specific heat in the limit 
k@/e >> 00, corresponding to (3.15), are 


E/N=(m+1)kT—4k02, | | is 


C/Nb=(m+1)—}d22/d(T/6). 


In the low temperature limit the internal energy and specific heat are 
easily calculable for an energy band of general form, such as that which 
has been obtained for nickel. If ¢), and «9, denote the energy widths of 
the unoccupied states for the two spin directions at 0° K. and «9 the 
corresponding width for zero spontaneous magnetization, then 


E 


Pieide. (A3 
5 in F v(e) de (A.3) 


0 
giving €y;, €92 a8 functions of Co, the relative magnetization. The internal 
energy is therefore an implicit function of C5, 


Eo 


Eanes | Paneer [oe TEEN TE 2, coe (A 4) 


0 


and the condition for E to be a minimum gives the relation between Cy 
and k6’/e, as well as the values of €; and éq9. The low temperature 
specific heat is then given by 


C= hr? 7T fy(€9;)+(eo2)}- = 2 2 ees (A.5) 
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If (ec) is given by (A.1) it is easily found that 


Ne = 3 (gers) {0+ Lorin. famines — 5. (6) 
0 
and 

aa {+6 prt tL (10) ae en 


The relative magnetization at absolute zero is found from (A.6) and is 
given by 
ko’ _ 1 1/m+1} - : 
— = s5 {(14 0) mt (1— Eg) m4}; . . . (A.B) 
€0 200 
hence the condition for ferromagnetism is k6’/e9>1/(m+1) and the 
range of incomplete saturation at absolute zero (0 <<, <1) is given by 


1 kd’ 
rea —ml (m+), af gl ee eee 9 
Sag (A.9) 


The relations given above for rctangula and parabolic bands may be 
derived at once by putting m—0, $. 
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SUMMARY. 


A derivation is given of a form of the virial theorem in quantum 
mechanics suitable for the discussion of a system of particles in a box. 


Ir T and V are the kinetic and potential energies of a system, and r; is 
the position vector of the ith particle, the virial theorem in classical 
mechanics states that 


MoTs —2; . WEN . . . . ° . ° (1) 


where the bars denote time averages, and the R.H.S. is known as the 
virial. This result has been used in the discussion of equations of state 
(cf. Jeans 1925). In the discussion of the thermal equation of state 
derived from quantum dynamical principles, de Boer (1949) has used the 
virial theorem, a procedure which is only justifiable if it can be shown 
that it holds in quantum mechanics. In view of the recent controversy 
concerning quantum mechanical equations of state (cf. Green 1949, 
Uhlenbeck and Riddell 1950) it seemed particularly desirable to examin 
the validity of this assumption. \ 

de Boer justifies the use of the theorem in quantum mechanics by 
referring to Slater (1933), who concluded that it is valid ifthe time averages 
of classical theory are replaced by quantum mechanical averages. A 
similar conclusion had been reached by Finkelstein (1928). It appears, 
however, that these proofs apply only to systems where there are no 
infinite discontinuities of potential, and are thus not applicable to a 
system in an enclosure where the potential is infinite at the walls. As it 
is under just these conditions that the theorem is applied in the discussion 
of equations of state, a further proof is required. The subject has been 
touched upon again in more recent publications, though only incidentally, 
and these contributions will be discussed in relation to our treatment. 

After some manipulation of the Schrodinger equation, Slater obtains 
the result 

2 OM Ors yp - ( ag a 
2 ay, 2 (Ogata, ~ Be a) tH" (28g, )¥=% =e) 
ee > oh ee SA eo Lege SA Te 
+ Communicated by James Taylor. 
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where the symbols have their obvious significance. To integrate this, 
use is made of the identity 


eee 


ey * OMY Os OPyp* — 9ify* ony ase yor 2 5) On; 
*3 (v Ox?dx, Ox; Ox? io Ox? i Ou, . b* 


(3) 


The second term on the R.H.S. is set equal to zero after integration with 
respect to x,, because s*2 is zero at the limits. It may be seen on 
performing the differentiation inside the brackets, however, that the term 


[ 25,52 FE |. a 


1 as 00 


where /, to /, is the range of integration, remains, and that this is zero 
only if os*/dv; vanishes. In free space where the range of integration 
is —oo to +o, the derivative of the wave function vanishes at the 
limits, but this does not happen when we are concerned with a system 
inside a box. Here % is typically a sine function with a non-vanishing 
derivative at the limits. This term must be concerned with the virial 
of the external forces acting on the system by virtue of the walls, and 
we shall now show that, if @ is a parameter defining the scale of the box, 
the neglected term is equal to a(dE/da), giving the final result 
4 ay) 


oF — aa, tar VN, 1 (8) 


which is the form in which the theorem has, in fact, generally been used. 


/ : 9 Ty 
If HaT+V— 3 Ve re 
and the wave equation Hj—Ey and its complex conjugate hold, it is 
easy to show, by differentiating with respect to a, keeping r,; constant, that 


ett: ae J(#s eV), Be ae etiy 3| 


re 


In the above, as in all following equations, the volume integral is over 
the volume of the box for each particle. This volume integral may be 
transformed into a surface integral using Green’s theorem, the surface 
being that corresponding to the volume through which the r, integration 
is carried. This allows one of the terms to be dropped because b* is 
zero when any position-vector r; terminates on the surface S,;.. We 


obtain 
OK ~ _ fh? Ou 
da om 2m, | Le Wah"). d 8; | dz) en a een) 
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where dr has been replaced by dr; . dz; in the usual convention. 
Similarly, starting with the wave equation and proceeding in the 
same manner as Slater, we obtain 


o ade aie Re. 
br, V;Vap*r, . V; (z- a vib) —(2- om viv") r, . Vib 
(9) 


which is equation (2) expressed in vector notation. It can easily be 
shown by expansion that 


rae h? h? 
ar, . V;{ 2— Val |= 2 — PRA ero se Ve 
oN i(2 a =e Tah (ss,.¥) 2g Badin nts 

(10) 
On substituting in (9) we obtain, after integration throughout the 
available configuration space, 


—2T +51, .V,V 


hi2 
=—2— 5 [lWtVien;. Vil) (Bey VV] > 
hi2 
=E— 5 | UlEn; . Vil) Vab* 48] dr ar: sine deatti() |) 


If a defines the scale of the box then it may be shown (see Appendix) 
that 2; . V j= —a(0}/da). Comparing this result with equation (8) we 


Be aciude that 
— ———— db 

bee : Nee ee es (12) 
which is the form in which the virial theorem is usually used. 

Since this work was carried out, it has come to our notice that Price 
(1950) has derived an expression for the stress tensor at the boundary in 
a fluid using a similar transformation to a surface integral.. By applying 
_ Fréhlich’s method of boundary perturbations, he obtained a result which 

can, by combining his equations (25) and (39), be put into the usual 
form of the virial theorem. 

It is true that if the virial in Slater’s discussion is regarded as the sum 
of the virial of the external and internal forces, the result is formally 
correct, but because of the infinite discontinuity at the walls, the integration 
of equation (2) throughout infinite space cannot properly be carried out, 
as has been pointed out by Uhlenbeck and Riddell (1950). These authors 
attempted to get round this difficulty by considering a smooth though 
nearly step-like potential at the wall of the box. This does not define 
the box properly, as pointed out by Green (1950), but it is possible to use 
Uhlenbeck and Riddell’s analytical method together with an infinite 
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potential wall, to give the correct result. Their method, which has also 
been applied by Zwanzig (1950) in a similar connection, used a device 
due to Green (1949) to differentiate the Hamiltonian, but as there seems 
to be some confusion about the Hermitian character of one of the terms 
involved, and also because this seems to be a more elegant, if less super- 
ficially convincing method, we apply it here. 

In Cartesian coordinates the Hamiltonian is 


He Tey eta (a 1V(x z,) (13) 
aver? (24242) 7 


We introduce new variables €;,7;, ¢; defined by x,=a€;, y;=an;, 2;=al % 
and differentiate the Hamiltonian with respect to a, keeping €,n, ¢ constant. 
Thus 


dH Re (/@ 8 B\ av 
Ba 2 2 Sina (368 + Ore aa) ta 8 

and so 
a5, = —2T+Er,. WY, Pee ee 7s) 


a fy Jq bp r=—2T + Ee, AR ho. le 


By differentiating the Schrédinger equation, and integrating the result, 
we have 


a a= [eS Hp dr+ [ye B) Har, any? aes 


Since (0//a). ,,-=0 at the wall (see Appendix), the second term on the 
R.H.S. is Hermitian, and so disappears, giving a result which may be 
combined with (16) to give the virial theorem (12). Green (1950) has 
stated that when the potential at the wall is infinite, this term is not 
Hermitian because 0y/da is not zero at the surface. This seems to be 
due to a confusion about what is kept constant in the partial differ- 
entiations. In the present treatment we are concerned with (0%5/0a); ,, , 


which is zero at the wall. (d/da),,,,, on the other hand, which we 
used in the previous treatment, is finite at the wall. 


A DabutiNGDPeX 
Let the equation of the box in rectangular Cartesians be 


SiS, uGleOian ee ats Ge meres) 


where =z/a, etc. We can also write 


p(x;, Yn Cae a) — ol Nis ie a), Say. . e . e e e e (19) 
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Then provided é,;,7;, ¢; satisfy (18) for any single 7, 7. e. provided the ith 
particle is on the surface, Y is identically zero. 


2 "Vy 


ase = =2(e i +u5t +2 5) ass Zr, . Vib+a (=) 


II 


0 when SH Siti eG ams a ce (20) 


X,Y, 2 


(21) 
The R.H.S. of (21) is, therefore, also zero. 
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XLV. The Processing of Thick Photographic Emulsions. 
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[Plate XVIT.] 


SUMMARY. 


An attempt has been made to determine satisfactory processing 
conditions for Ilford G5 electron-sensitive emulsions up to 600m in 
thickness. The penetration and temperature characteristics of metol- 
hydroquinone, para-aminophenol, and p-diaminophenol developers have 
been investigated. The observations show that a p-diaminophenol 
developer gives satisfactory results and a suitable procedure for processing 
is described. In addition the problems of fixing, washing and drying, 
have been studied. 


INTRODUCTION. ; 

PHOTOGRAPHIC emulsions have been used as detectors of charged particles 
ever since Reinganum, in 1911, discovered that an alpha-particle was 
capable of producing tracks which were observable under a microscope. 
The early emulsions did not possess great sensitivity; consequently 
accurate measurements of the properties of nuclear particles could not be 
made. In 1946, however, Ilford Ltd. in Great Britain produced ** Nuclear 
Research ’* emulsions containing about eight times the normal amount 
of silver bromide per unit volume (Powell et al. 1946). These emulsions 
were capable of recording the tracks of particles with a specific ionization 
of five or more times the minimum value. Their thickness was usually 
50-100 and no particular difficulties were encountered during processing. 
Owing to the high concentration of silver bromide, however, it was found 
advantageous to agitate the fixing solution and to change it fairly 
frequently (Powell and Occhialini 1947). 

When emulsions of 150, thickness or greater were produced, the 
application of the processing methods previously employed resulted in 
a large variation of development with depth. This was found to be due 
to the fact that the developer took a considerable time to penetrate the 
emulsion. This difficulty was overcome by Dilworth, Occhialini and 
Payne (1948) by the method of * temperature development ” in which 
the developer is allowed to penetrate the emulsion at a low temperature, 
and is subsequently warmed to allow the main development to occur, 
Details have been given by Wilson and Vaneslow (1949), while the results 
of an extensive investigation have been published by Dilworth, Occhialini 
and Vermaesen (1950). 

An alternative approach has been suggested by Blau and De Felice 
(1948). They first soak the emulsion in a developing solution made up 


* Communicated by Professor ©. F. Powell, F.R.S. 


On the Processing of Thick Photographic Emulsions 397 


without alkali. When permeation is complete the plate is transferred to 
developer containing excess alkali. This is known as the “ two bath ” 
method. For the successful application of this method to thick emulsions, 
it is necessary that the speed of travel of a pH change through gelatin 
from the surface should be much greater than the rate of diffusion of the 
developer as a whole. This condition is not satisfied (Waller 1951). The 
“temperature development ’’ method has therefore come into general 
use, and forms the subject of the work to be described in this paper. 

In 1948. Kodak Ltd. produced an emulsion, 200 » in thickness, capable 
of recording the tracks of particles at minimum ionization (Berriman 
1948). Similar emulsions are now manufactured by Ilford Ltd. The 
production of such emulsions represented an important technical advance 
and has made possible a considerable extension of the field of research 
in nuclear physics. The analysis of the tracks of charged particles by 
a study of their multiple scattering and grain density requires tracks of 
great length in the emulsion. With the early thin emulsions, the number 
of tracks of sufficient length was, for simple reasons of geometry, relatively 
small, and it became desirable to employ emulsions of greater thickness. 
The first successful steps in this direction were taken by Dilworth et al. 
(1948). 

For the successful application of the scattering method, distortion of 
the emulsion must be reduced to a minimum, for at the present time 
errors due to this cause commonly determine the maximum energy of 
a particle of given track length which can be measured, and contribute 
to the errors of any given measurement. These considerations emphasize 
the importance of processing emulsions of relatively great thickness with 
a minimum of distortion. 

In this paper we give an account of the processing methods which have 
been found satisfactory for Ilford G5 emulsions up to 600 « in thickness. 
These methods are not necessarily the best or the only ones available but 
they may be employed in dealing with large numbers of plates. It 
must be emphasized that each batch of plates is slightly different from 
any other, and the optimum processing conditions must be found by 


- trial and error. 


$1. THH DEVELOPING AGENT. 
General Considerations. 


Following Dilworth et al. (1950) we may specify several conditions 
which an ideal developer should satisfy. For the “ temperature 
- development ’”’ method a satisfactory developer should :— 

(i.) be capable of penetrating the emulsion reasonably quickly. 

(ii.) possess a moderate or high-temperature coefficient. 
(iii.) produce negligible background fog and stain in the emulsion. 
(iv 
(Vv 
SER. 7, VOL. 42, NO. 327.—APRIL 1951 


-) 

) 

.) be stable at temperatures up to 30° C. 
:) 


have a pH as close as possible to that of the emulsion. 
25 


398 A. D. Dainton, A. R. Gattiker and W. O. Lock on the 


Each of these factors will be considered in turn. The developers 
chosen for investigation were: D19b (metol-hydroquinone), Azol (para- 
aminophenol), Amidol (p-diaminophenol) and Amidol-bisulphite. Details 
of the constitution of the solutions are given in Table I. 


The Characteristics of the Developers. 


(i.) Penetration characteristics. 

There has been little published work on the penetration times of various 
developers in “ electron-sensitive ’ emulsions. Payne (1949) has shown 
that for thin emulsions the rate of diffusion of any developing agent 


29 


TABLE I. 


The Developing Solutions Investigated 


Name D19b Azol Normal amidol Amidol bisulphite 


Developing Metol-hydro- Para-amino- p-diamino- p-diamino- — 
agent quinone phenol phenol phenol 
Formula | Metol 2:2 g. Johnson’s Azol solu- | Amidol 3-0 g. Amidol 3-0 g. 
used Hydroquinone 8:8 g. tion 16 mls. Anhydrous sodium | Anhydrous sodium 
Anhydrous sodium | 1 per cent potassium sulphite 12-0 g. sulphite 6-7 g. 
sulphite 72-0 g. bromide 88 mls. Distilled water (add) | Sodium _bisulphite 
Sodium carbonate | Distilled water (add) 1000 mls. liquor (8.G. 1-:34)* 
48-0 g. 384 mls. 1-4 mls. 
Potassium bromide Distilled water (add) 
40g. 930 mls. 
Distilled water to 
2000 mls. 


The above stock 
solution diluted 
1: 1 with distilled 
water 


pH 10-0 


— 
= 
or 

1 

bo 

fon) 

a | 


* The sodium bisulphite liquor is a commercial preparation supplied by British Drug Houses Ltd. 


increases with decreasing pH of the solution ; and that pre-soaking the 
emulsion in water reduces the developer penetration time by as much as 
40 per cent ; see also Dilworth et al. (1948). 

The method adopted for the determination of the penetration times of 
various developers in Ilford G5 emulsions was the following. Plates 
coated with emulsion of different thicknesses up to 1000 were cut 
into pieces of dimensions } in. x #in. The back of each plate was exposed 
to the image of a characteristic pattern, produced by a photographic 
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enlarger. As the light impinged on the emulsion from the back of the 
plate, a latent image was produced in the emulsion at the glass-emulsion 
interface. The exposure time was such that a clearly defined image was 
produced at the bottom surface of an emulsion 5 « thick, when developed 
in the D19b developer at 20° C. for one to two seconds. Even at 5° C. 
the corresponding period is only about 70 seconds. The same exposure 
was given to all the emulsions, of various thickness, by means of an 
accurate photographic timer, the voltage supply to the light source and 
“timer” being accurately controlled. 

To ensure that the photographic latent image was confined to the 
glass-emulsion interface, the light from the enlarger was passed through 
a Wratten No. 35 filter before striking the plate. Such a filter transmits 
red and blue light but not the intermediate wavelengths. The silver 
halide is sensitive only to the blue light which is strongly absorbed, and 
the latent image formed is confined to a layer approximately 3 thick 
at the glass-emulsion interface. Thus, no trace of the geometrical 
pattern appeared on the surface of the emulsion 5 u thick exposed in this 
way. The first appearance of the pattern at the bottom of an emulsion 
is therefore a good indication of the time the developer has taken to 
penetrate the full thickness of the emulsion in sufficient strength to 
develop the exposed grains of silver halide. For the purpose of the present 
experiment, we may define the penetration time as the period the developer 
takes to penetrate to the bottom of the emulsion and develop the image 
of the geometric pattern. It may be emphasized that the penetration 
time so defined is not appreciably affected by the variation with 
temperature of the time taken to develop the image once the developer 
is present at a given concentration ; for all temperatures this time is 
relatively very short. 

Details of the emulsion thicknesses used are given in Table IT. 


TABLE II. 


Emulsion Thicknesses Used. 
Nominal Thickness Actual Thickness 


50 w 45 pw 
100 105 pe 
200 pu 230 pu 
300 260 je 
400 400 pw 
600 675 pu 
800 ju 720 pu 

1000 1050 pu 


The measurements were first carried out at 20° C. with the different 
developing solutions. One complete set of plates was first pre-soaked in 
distilled water for three. hours, whilst a second set was not pre-soaked. 
The results obtained are given in Table III. and in fig. 1. 


4Z4h2Z 
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TABLE IH. 


Penetration Times of the Different Developing Solutions at 20° C., in 
Pre-soaked and non Pre-soaked Emulsions. 


Developer Soaking Penetration Time in Minutes for Different 
Thicknesses 

, 105 p 2302 260 400 675. 720p 1050 p 
Azol Pre-soaked 44 17 22 38 96 96 210 
Non Pre-soaked 6 21 27 00) ISS sel 330 
D19b Pre-soaked 34 94 144 21 50 64 140 
Non Pre-soaked 5 12 20 37 110" 125 270 
Normal Pre-soaked 14 6 84 12 26 26 58 
amidol Non Pre-soaked 34 9 114 22 49 52 120 
Amidol Pre-soaked 2 64 8 12 25 25 55 
bisulphite Non Pre-soaked 2 84 10 204 = 53 46 110 


All pre-soaked emulsions soaked in distilled water for three hours irrespective of 
thickness. 


Pints gl 


in MINUTES 


PENETRATION TIME 


| 
2 ry Cars 
100 1000 


EMULSION THICKNESS in MICRONS 


Penetration time of the developer as a function of emulsion thickness for normal 
amidol at 20° C.; curve A for emulsion without pre-soaking; curve B 
emulsion pre-soaked in distilled water at 20° C. for three hours. 


These show that pre-soaking the emulsion in water considerably reduces 
the penetration time of the developer. The ratio of average penetration — 
times, without and with pre-soaking, for the various developers is :— j 
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Azol 1:34-L0:05 
D19b 1-70 +0:°13 
Normal amidol 1-86 +0-12 


Amidol bisulphite 1:65+0:14 ' 


Tt will be seen that the amidol and amidol bisulphite developers are 
equivalent in penetrating power, and they penetrate the emulsion more 
rapidly than D19b or Azol. 


200 


in MINUTES 


PENETRATION TIME 


2 4 Gi 1006 


© EMULSION THICKNESS in MICRONS 
Penetration time of the developer as a function of emulsion thickness for normal 
amidol at 5° C., 10° C., 15° C., and 20° C. All emulsions pre-soaked in 
distilled water for three hours at the appropriate temperature. 
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In view of these results, the penetration measurements at 15° C., 
10° C., and 5° C. were carried out, using only the first three developing 
solutions specified in Table I., and always pre-soaking the emulsions in 
distilled water for three hours at the appropriate temperature, whatever 
the thickness of the emulsion. The results obtained are displayed in 
Table IV. and plotted graphically in figs. 2 and 3. The penetration 
time, T, may be expressed in terms of the thickness of the emulsion, t, by 
a relation of the form 

Ta-keite, 
where x is of the order of 1-4, and k is a constant for a given developer 


at a given temperature. 
Fig. 3. 


PENETRATION TIME in MINUTES 
o @ 
T T 
| oy 
| > i 
N 
te) 
FS 
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| 2 4 6 68 
100 1000 
EMULSION THICKNESS in MICRONS 


Penetration time of the developer as a function of emulsion thickness for normal 
amidol, D19b and Azol developing solutions at 15° C. All emulsions pre- 
soaked in distilled water for three hours at 15° C. 


It must be emphasized that T is the penetration time as defined above 
and not necessarily the time required for.satisfactory processing. The 
latter periods are given in Table VIII. 

(ii.) Temperature characteristics. 


It was also desirable to determine the temperature characteristics 
of the different developers in the G5 emulsion. Strips of dimensions 


et 
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g In. xX? in. were cut from 2 in. 2 in. plates with emulsion 5 y thick. 


These were exposed to light passed through a step wedge (15 steps of 
density range 0-2) held in a photographic enlarger, the image of this 
step wedge being formed on the glass-emulsion interface of the plate. 
The exposure was accurately controlled as in the previous experiment 
and was such that the image of the middle portion of the wedge was 
satisfactorily developed using the D19b developer at 20° C. for 10 minutes. 

For any one developer at any one temperature, five emulsion strips were 
used. These were immersed in the developer in question for various 
times, depending on the normal time used for that developer. For 
example, the D19b developer is usually used for 20 minutes at 20° C. ; 
thus strip no. 1 was removed from the developer and put to fix after 1} 
minutes, no. 2 after 3 minutes, no. 3 after 7 minutes, no. 4after 15 minutes, 
and no. 5 after 30 minutes. At a lower temperature, these times were 
increased, e.g. at 5° C. using the same developer the times were twice those 
used at 20° C. 


Tass LV. 


Penetration Times, in Pre-soaked Emulsions, of the Various Developing 
Agents at 5° C., 10° C., 15° C., and 20° C. 


Developer Penetration Time in Minutes 
105 uw 230 260 400 675 720 1050 
Azol De Oe it 45 50 100 270 #300 300 
hoes 84 9-27 34 70.) L85e 208 210 
15°C: 8$ 25 30 63 145 150 210 
DO- C. 44 17 22 38 96 96 210 
D19b uae 6 21 27 58 165 172 270 
107 C: 5 164 = 19 S712 s1085 S114 210 
Loa: 4 [ke ao ae ae) 27 72 68 220 
20° GC. 34 9$ 144 21 50 64 140 
Normalamidol 5° C. iis 16 20 a 80 84. 190 
10° G: a+ oi 134 22 51 534 ~=—-110 
157. C: 3 8 114 18 374 45 95 
20° C. 14 6 84 = 12 26 26 58 


The figures for 20° C. are given again for comparison purposes. 


The above procedure was carried out at 5° C., 10° C., and 20° C. for 
each of the four developers under consideration. When all the plates had 
been processed and dried the density of each step in each wedge was 
determined using an Ilford densitometer and recorder. 

The usual method of determining temperature coefficients was found to 
be inapplicable in this experiment owing to the extreme contrast of the 
plates. Instead, the times to produce equal densities on the same step 
(no. 6 in this case) of each wedge were used as a basis of comparison between 
developers. Equivalent readings were not available in the case of Azol, 
which was found to be incapable of giving images of sufficiently high 
density. The results are given in Tables V. and VI. and are presented 
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TABLE VI. 
Development Time for Unit Density as a Function of Temperature. 
Developer Time for Density Figures Normalized 
1 in Minutes to Density 34 at 20° C. 
D19b jr — — 
10° C. 30 30 
150; 8 8 
20° C. 34 34 
Normal amidol Oe.U: 45 45 
10° C. 18 18 
15°C. = —- 
20° C. 34 34 


Amidol bisulphite oC, — aa 


Azol 


DENSITY 


DEVELOPMENT TIME in MINUTES 


Density of step 6 of the step wedge as a function of development time for amidol 
bisulphite developer at various temperatures. 


graphically in figs. 4 and 5. It may be seen that the rate of action of a 
the developers changes rapidly with temperature. The i ag i 
suggest that D19b is more sensitive than the amidol developers. zol, 
although very sensitive to temperature, has very low activity. 


(iii.) Background formation and staining properties of developers. 

Background grains not associated with recognizable tracks are ie in 
part, to particles traversing the emulsion at a large angle of gies ° fe 
grains arising from (a) chemical action during the preparation 0 
emulsion, and (b) the action of the developing agent. 
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A measure of the number of background grains per unit area was 
obtained by counting the number of such grains in focus under specified 
optice! conditions. The microscope used was a Cooke M4000 with 
x95 objective and 15 eyepieces. Some results for emulsions exposed 
and processed under different conditions are given in Table VII. 
A background of as much as 80 grains per 10,000 u? can be tolerated 
in emulsions in which the grain density in the track of a particle with 
minimum specific ionization is of the order of 15 grains per 50 p. 


in MINUTES 


DEVELOPMENT TIME 


S 10 IS 208C 
TEMPERATURE 


Time for unit density of step 6 against temperature, for amidol type and 
D19 b developing solutions. The figures are normalized at 20° C. 


TaBie VII. 
Background Percentage of 
a fies. wey eee Grains/10,000 p22 Light 
Exposure Developer Stop Bath Su from 5, trom ‘Transmitted 
surface glass 

Ll. G5 400 ph Azol 5 per cent 35+2 Too misty 27 
1} hrs, at 68,000 ft. sodium 

bisulphite 
2. G5 400, Normal pH3 HCl 95+5 7443 a 
53 hrs. at 64,000 ft. amidol 
3. Gd 400. Amidol 2 per cent 34+2 33+2 40) 
55 hrs. at 64,000 ft. bisulphite acetic acid 
4. G5 500 pu Amidol 4 percent 59+3 6143 25 


Glass 2 mm. bisulphite acetic acid 
6 hrs. at 75,000 ft. 
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The degree of stain in an emulsion is difficult to define quantitatively. 
We may estimate the combined effect of background and stain by means 
of an Ilford densitometer which measures the percentage of light trans- 
mitted by the photographic plate. The figures given in Table VII. show 
that the stain is not necessarily related to the background grain density. 
It may also be seen that the use of amidol bisulphite developer results in 
plates of good optical quality. 


(iv.) Stability of the developer as a function of temperature. 


Following Dilworth et al. (1950), who drew attention to the possibility 
of employing developing solutions at’ bigher temperatures than those 
commonly used, we have found thet very satisfactory results are obtained 
when a temperature of the order of 27° C. is used in the hot stage of the 
temperature cycle. For such work it is necessary that the developer 
should remain stable at these temperatures, and all the developers given 
in Table I. are satisfactory in this respect. 


(v.) pH considerations. 


It has been shown that when Ilford Nuclear Reseerch emulsions are 
immersed in a solution the swelling of the emulsion decreases as the pH 
of the solution is reduced. In fact, the emulsion will even shrink when 
placed in dilute acid (Waller 1951). These considerations make it 
desirable that the developer employed should be slightly acid, for the less 
the swelling of the emulsion, the smaller the distortion introduced at this 
stage of the processing. 


CONCLUSIONS. 


From a consideration of the factors discussed above, it appeared to us 
that amidol promised to give the most satisfactory results. It was 
therefore adopted as a standard developer for emulsions up to 300 p in 
thickness. For thicker emulsions it was found that normal amidol 
produced a slight variation of development with depth. We have 
succeeded in reducing this variation by adding sodium bisulphite to the 
_ developing solution (Balagny 1912). Balagny employed normal photo- 
graphic emulsions and his recommended value for the pH of the developer 
was not found to be suitable for thick electron sensitive emulsions. 
Using the amidol bisulphite developing solution given earlier, with a 
pH of about 4, we have found it possible to develop the lower regions of an 
emulsion without developing the surface layers. The pH of the solution 
which gives satisfactory development throughout the entire thickness of 
the emulsion is 6°7. 

Experiments were also carried out in an attempt to reduce the formation 
of fog in the emulsion by reducing the concentration of amidol, but tracks of 
particles of minimum ionization were not recorded for concentrations less 
than 3 gm./litre. 
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§ 2. THE DEVELOPMENT PROCESS. 


We may now consider the details of the development procedure. 
The results given in the previous section suggest that : (a) the emulsion 
should be pre-soaked in water before being placed in the developing 
solution ; and (6) that the developer used should be amidol in some form. 

In our procedure the plate is therefore first immersed in distilled water 
at room temperature. The water is then cooled to 5° C. and maintained at 
this temperature for a time which depends on the thickness of the emulsion. 
It is then transferred to the developing solution, which is also at 5° C., and 
remains in this solution until the developer has completely permeated the 
emulsion. During this time little development will have taken place. 

For the main development process two alternative methods may be 
employed. The developer and plate together may be warmed 
to about 18° C., and the temperature maintained at this value 
for the time required to give satisfactory development. Alternatively, 
following Dilworth et al. (1950), the plate may be separated from the 
developing solution, any excess solution being removed by means of filter 
paper, and the plate then placed glass downwards on a dry surface. This 
surface is then brought to a temperature between 25 and 30°C., and 
maintained at a constant value for the requisite time. In this method 
one thus relies upon the developer already present in the body of the 
emulsion and avoids over-development at the top of the emulsion due to 
the entry of fresh developer. This is known as the “ hot-plate ’ method. 

In our experience the second method gives the more satisfactory results. 
‘To obtain tracks of particles with minimum ionization, the actual value of the 
temperature of the ‘* hot plate’ does not appear to be critical. If, however, 
it is important to preserve uniformity of development among plates from 
different batches it is necessary to introduce an accurate system of tempera- 
ture control. Care must also be taken to ensure good thermal contact 
between the plates and the “‘ hot plate ”’. 

If amidol is used as the developing solution, the surface of the emulsion, 
after processing, is usually covered with a deposit of silver. This may 
sometimes, although not always, be swabbed off with cotton wool soaked 
in methylated spirits, when the emulsion is thoroughly dry. For some 
experiments, such as those with artificially accelerated particles, it is 
important to know the exact point where a particle entered the surface. 
Any surface silver deposit is therefore to be avoided if possible. It has 
been suggested by Dilworth et al. (1950) that developing the plates in 
an inert atmosphere may lead to a diminution in the amount of the surface 
deposit. Alternatively, for emulsions less than 400 « in thickness, azol 
may be employed, using the first of the development methods described 
above, since it commonly leaves little or no surface deposit. 

The next problem is to stop the development at the correct moment. 
The general procedure is to transfer the plates to a ‘‘ stop bath ”’ which is 
at a temperature near to that of the ‘‘ hot plate ’’; the temperature of the 
solution is then lowered to about 5° C. and kept at that value until all 
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development has ceased. The solution and plates may then be warmed to 
the temperature of the fixing solution, and after a short wash in water, 
transferred to it. Thus development is stopped both by the application of 
a stopping agent and by the use of low temperature. We have tried 
several stopping agents including distilled water, a 1 per cent potassium 
bromide solution, a5 per cent sodium bisulphite solution, } percent solution 
of hydrochloric acid and solutions of acetic acid. The best results were 
obtained with the acetic acid. The concentration did not appear to be 
critical, and we find that a 4 per cent solution is satisfactory. 

In all the processes mentioned above, the plates are kept horizontal. 
This is especially important when dealing with very thick emulsions. 


§ 3. FIXATIon. 

The critical importance of avoiding distortion has already been stressed. 
The main causes of distortion are the stresses introduced when the emulsion 
suffers volume changes including those during manufacture. These 
changes may be caused by mechanical agitation, temperature variations, 
pH changes, and the removal of the unsensitized silver halide from the 
emulsion. Since most of the silver halide, which constitutes a large 
fraction of the unprocessed emulsion, is removed during processing, it is 
probable that most of the distortion is introduced during fixation, washing 
and drying. After the fixing solution has been washed out and the plate 
allowed to dry, the emulsion has less than half its original thickness. 

The plate is fixed with sodium thiosulphate— hypo”. Camerini 
and Franzinetti (1947), working with thin emulsions, found that a 
concentration of 40 grams of hypo to 100 mls. of water gives the most 
satisfactory results. We have confirmed this result for emulsions 
400 w thick. 

Owing to the large amount of silver halide in electron-sensitive emulsions 
the fixing process takes a relatively long time. It is well known that 
ammonium thiosulphate acts rapidly as a fixing agent. The addition of a 
small amount of ammonium chloride to the hypo solution also shortens the 
total fixing time. We have found, however, that a fixing solution of this 
nature tends to remove some of the developed grains near the surface of 
the emulsion, especially in the case of thick emulsions where the total 
fixing time may be of the order of two or three days. This has also been 
observed by Wilson and Vaneslow (1949). Because of this effect we do not 
favour the use of ammonium chloride. 

On the other hand, we have found it advantageous to add sodium 
bisulphite to the hypo solution ; this leads to a reduction of the stain in the 
finally processed emulsion (Mees 1942). The fixing solution which has 
been adopted as standard in this laboratory is :-— 


Sodium thiosulphate 400 grams. 
Sodium bisulphite 30 grams. 
Water to add 1,000 mls, 


The pH of this solution is 5-0, 
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Agitation. 

It has been found by several workers that gentle agitation of the 
fixing solution leads to a reduction in the total fixing time (Powell and 
Occhialini 1947; Wilson and Vaneslow 1949). Itis essential that the plates 
are kept horizontal and that the agitation, if any, is very gentle. This is 
especially important for emulsions thicker than 400. If, for example, a 
500 4 emulsion is held vertically during the fixing process, the emulsion 
tends to ‘‘ creep’, and it is frequently observed that tracks near to the 
bottom are discontinuous. An example of this nature which we call 
‘chopping’ is shown in Pl. XVII.a. Similar effects sometimes occur 
in emulsions processed horizontally (Herz 1950) and have been seen in 
emulsions 200 yz thick. 

We have found that suitable agitation is provided by a slow steady 
flow of the fixing solution over the surface of the plates. Satisfactory 
results may also be obtained by placing the plates in a static fixing bath, 
provided the bath is changed fairly frequently. A fixing bath is 
considered to be exhausted when the concentration of silver in it rises 
above between five and seven grams per litre. 


Fixation time. 

The period for complete fixation is usually taken as one and a half times 
that for the emulsion to clear, that is, when the last trace of opaque silver 
halide has disappeared from the bottom of the emulsion. For a 400 u 
emulsion, the total fixation time is of the order of 24 hours, while for a 
1,000 . emulsion it may be as long as 100 hours. 


$4. WasHine anp Drytne. 
Washing. 

When all the surplus silver halide has been removed from the emulsion, 
the fixing solution must be washed out. Thin emulsions may be washed in 
a stream of tap water. When this method is employed with emulsions 
600 1 thick, however, the emulsion frequently distorts, and “‘ bubbles ”’ 
may be formed on the surface. 

We have therefore found it necessary to commence washing with water 
at a low temperature (5-10° C.). The plates are placed in a large dish of 
fixing solution and the cold water is allowed to trickle in slowly at: first. 
As the concentration of the fixing solution is reduced, the rate of flow of the 
water and its temperature may be increased. The washing process may 
take as long as two or three days, and it is only considered complete when 
the washing solution gives no indication of the presence of hypo as judged 
by the potassium permanganate test. 


Drying. 
Following a suggestion of Mr. C. Waller, we next soak the plates in a 
solution of glycerine for approximately an hour. This prevents the 


emulsion from stripping off the glass after it has dried. A 1 per cent 
solution has been found to be satisfactory for this purpose. We have 
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found some indication that the use of a 5 per cent solution tends to reduce 
the distortion, although the processed emulsion is then rather soft. 

After removal from the glycerine solution, the plates are placed side by 
side in the form of a square or a rectangle in an enclosed box. If they are 
of a small size it is of some advantage to place a “ guard ring” of wet 
emulsion around the plates to be dried (Dilworth 1950). It is more 
convenient, however, to use large plates, say 8 in. x 8 in. and then cut them 
to 6 in. x 6 in. after they have dried. The distortion effects at the edges 
are thus almost completely eliminated. 

A preliminary investigation has shown that distortion is reduced if the 
drying of the plates is not unduly prolonged. A batch of plates was 
allowed to dry for seven days at a relative humidity of 90 per cent, and 
then for a further seven days during which the relative humidity was 
gradually lowered to about 55 per cent. The distortion in these plates was 
a little greater than that in emulsions of equal thickness which had been 
dried over a period of three days. Drying by gently blowing air over the 
plates is undesirable and introduces severe C-shaped distortion. Warming 
the plates by means of a 25 watt electric light bulb at a distance of a few 
feet appears to give satisfactory results. 


§ 5. CONCLUSIONS. 
Distortion. 

The effect of distortion on the measurement of the scattering of a track 
has been briefly discussed by Fowler (1950), while Cosynsand Vanderhaeghe 
(1950) have described a method of estimating its magnitude. It is, of 
course, most severe at the edges of a plate ; examples of both C- and 
S-shaped distortion at the edge of a 400 » plate are shown in PI. I. 6 and c. 
The elimination of these edge effects has been discussed in the previous 
section. 

The distortion in the centre of a large emulsion is usually of the C-type 
and its effect on the scattering of a track may then be estimated. On the 
other hand, it may occasionally take the form of irregular and local 
‘‘ chopping ”’ for which no correction can be made. 


General considerations. 

The results given above have been summarised in Tables VIIT. and IX. 
which detail satisfactory processing methods for G5 emulsions of various 
thicknesses. 

The methods given in Table VIII. have been used in the processing 0 
Ilford C2, D1, and El emulsions 200 thick (Danysz 1950; Muirhead 
1950), Ilford C2 emulsions 600 4 thick (Smith 1950) and Kodak NT4 
emulsion 200. thick (King 1950).° Satisfactory results were obtained in 
all cases. 

In collaboration with Mr. M. G. K. Menon, and Mr. G. T. Zorn, the 
investigations described are being extended to still greater emulsion 
thicknesses, and to the diluted emulsions recently produced by Ilford 
Ltd, (Dodd and Waller 1950). 


. ae 
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TaBLeE VIII. 
Amidol Type Developer. 


Emulsion Thickness. 


100 200 400 pe 600 1006 pe 
Soak in distilled water 15 mins. 25 mins. 90 mins. 120 mins. 240 mins. 
Developer at 5° C. ED oeee WI Fe 90 _,, IWAU) as 300s 
‘> Hotplateta2keO rr 2b ae: DASH 9 ey 30 eee 30 me Jom 
Stop bath LOE ug, a0 cd OD eds 120, 180i 
Fixing 2hours 5hours 24 hours 70 hours 100 hours 
Washing LS oe ee OU) ee 12 a eee 2) es 
Drying Game. Lg ee 60° 55 LOO 20 eee 

TABLE IX, 


Azol Developer. 
Emulsion Thickness 


100 pw 200 pu 400 pu 
Soak in distilled water 15 mins. 25 mins. 90 mins. 
Developer at 5° C. 30 90. 120 a8 
Warm stage 18° C. LOD* 5 1055 105 gee 
Stop bath aN 33 AD. ., OSs 
Fixing 2 hours 5 hours 24 hours 
Washing Ase 122 GOmmee 
Drying Gaur 2s OO ae 


The stop bath may be 5 per cent sodium bisulphite or } per cent acetic acid. 
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(@) (0) (c) 
Examples of (a) “chopping”, (6) C-shaped distortion, and (c) S-shaped 
distortion in 400 » G5 emulsions. 
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SUMMARY. 


A general relationship between stress and plastic strain in a 
polycrystalline aggregate is derived for any metal in which individual 
crystals deform by slipping over preferred planes under a critical shear 
stress. Full account is taken of the non-uniform distortion due to mutual 
constraints between the grains of an aggregate. It is shown that a 
plastic potential exists which is identical with the yield function. Upper 
and lower bounds are obtained for an approximate calculation of this 
function for any applied system of combined stresses. 


§1. INTRODUCTION. 


THE macroscopic theory of plastic deformation in a polycrystalline metal 
is based on observations of the behaviour of the metal in bulk. It rests 
also on the simplifying hypothesis that the material is _ locally 
homogeneous (though not necessarily isotropic). The theory so 
constructed is found to be adequate as a first approximation when 
applied to many problems in engineering science and metal technology. 
There remains, however, the task of relating the macroscopic observations 
to ones more fundamental, for example, the mechanisms by which a 
single crystal deforms under load. Such an investigation might hope to 
attain two main objectives: first, to show which of the several 
microscopic modes of distortion are chiefly responsible for the behaviour 
of the aggregate as a whole ; and second, to indicate, more precisely than 
experiment at present can, what refinements should be added to the 
macroscopic theory and what these might be for any particular metal. 
Several previous writers have dealt with related problems, and the most 
important papers will be discussed briefly. 

In single crystals of many metals it is known that the main mechanism 
of plastic distortion on a microscopic (though not atomic) scale is simple 
shear parallel to preferred planes and directions. At ordinary temperatures 
both the planes and directions are generally those of closest atomic 
packing. It is observed that pronounced slip (or glide) is initiated along 
the particular plane and direction for which the corresponding component 
of shearing stress first reaches a certain critical value under increasing 
external load. This critical value depends on the temperature and rate 
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of straining, but not on the crystal orientation or the type of load 
(provided the loading-path is not such as to induce hysteresis effects). 
This has been established for metals such as copper, aluminium, nickel, 
magnesium, zinc and cadmium by the experiments of Mark, Polanyi 
and Schmid (1922), Taylor and Elam (1923), Polanyi and Schmid (1923) 
and Schmid (1924), among others. Double slipping is observed if the 
critical shear stress is reached simultaneously on two of the possible 
planes and directions. 

The problem of predicting the tensile yield stress of a polycrystalline 
face-centred cubic metal was considered by Sachs (1928) and by Cox and 
Sopwith (1937) on the assumptions that (i) each grain is subjected to a 
uniaxial stress parallel to the specimen axis sufficient to initiate slipping 
in the most highly stressed direction ; (ii) all orientations are equally 
likely ; (iii) each grain is in the same state of work hardening. The 
resultant axial load on the specimen was calculated as the sum of the 
individual loads in the constituent grains. The value obtained for the 
macroscopic uniaxial yield stress was 2-27 (approximately), where 7 is 
the shear yield stress of a single crystal. This treatment is open to the 
objections that the grains could not form a coherent whole if only a single 
homogeneous glide occurred in each, and that the necessary conditions 
of equilibrium could not be satisfied across grain boundaries if the stress 
in each were a simple tension of varying amount. 

Kochendérfer (1941) added the further hypothesis that the grains 
extend by the same amounts in the direction of the specimen axis. This 
enabled bim to calculate the stress-strain curve of the aggregate. For in 
each grain 

o/7=dy/de=m 
where o is the axial stress in a grain, de is the common axial strain- 
increment, dy is the increment of shear strain in the slip direction, and 
m is a dimensionless factor depending only on orientation. (In passing, 
it may be noted that the virtual work equation ode=rdy is naturally 
satisfied). We now interpret the various symbols as mean values over 
all orientations, so that m=2-2 from Sachs’ calculations. It follows that 
the stress-strain curve of the aggregate is 


o=mr(y)=m7(me) 


where 7(y) is the shear-hardening curve (the transition from dy=mde 
to y=me is legitimate since m is independent of the amount of strain 
on the assumption of continuing isotropy). However, the best fit of 
this o(e) relation to Taylor’s experimental curve for aluminium gives a 
value of about 3-1 for m. 

In a recent paper (1950) Calnan and Clews considered the particular 
directions in which a pure tension would need to be applied to operate 
simultaneously either 4, 6, or 8 shears in a face-centred cubic lattice. 
They assumed the stress state in each grain to be a uniaxial tension of 
this kind and of the requisite amount. Each grain is assumed to undergo 
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the specimen extension (and not the strain as a whole). Thus, the 
objections raised against previous work with regard to the independent 
action of the grains apply here, also. Despite this Calnan and Clews’ 
computations led to a stress-strain curve in fair agreement with that 
measured for polycrystalline aluminium. 

A more realistic calculation of the tensile yield stress of a face-centred 
cubic aggregate is due to Taylor (1938), who assumed that each grain 
undergoes the same uniform strain (so maintaining cohesion). This 
requires, in general, the operation of at least 5 independent shears out of 
the 12 operable at ordinary temperatures. Hence, if only geometry is 
considered, many choices are possible for a given strain. Taylor 
introduced the hypothesis that the actual active set of shears is that for 
which the sum of their separate magnitudes is least. The hypothesis 
has no obvious a priori justification, but Taylor based it on observations 
of single crystals under uniaxial stress and on a postulated analogy with 
the dynamics of a non-conservative mechanical system. By adding the 
further assumption that the shear hardening of a crystal during 
simultaneous glide depends only on the sum of the absolute shears, 
Taylor was able to compute a stress-strain curve for an aluminium 
ageregate in fairly close agreement with the measured curve. There 
are two main criticisms of this approach: first, it is not proved that 
a combined stress could always be found to operate any geometrically 
possible set of shears (without exceeding the critical shear stress 
in non-active directions); and, second, no account is taken of 
stress-continuity conditions across grain boundaries. 

All these previous investigations have been confined to an aggregate 
under simple tension or compression. In the present paper no such 
restriction is placed on the applied load. 


§2. THE Macroscopic THEORY OF PLASTICITY. 


In order to appreciate the aim of the present work it is necessary to 
be familiar with the premises of the macroscopic theory of plasticity. 
These will be summarized here; a full discussion has been given by 
Hill (1950). 

The theory assumes a locally homogeneous material (possibly 
anisotropic) having a sharp yield point after work-hardening. The 
yield criterion for a metal in a given state of hardening is assumed not 
to depend on the hydrostatic component o of the applied stress o,,, and 
is usually taken to have the form 


f(oj)=¢3 oj=o,—08 4. 5 sey ite. ~~! CY) 


This equation ih be regarded as defining a surface in stress nabs 
the surface is ‘‘ cylindrical ’’, with generators parallel to the direction 85 
An assumption adequate for many applications, though not necessary in 
principle, is that the function f does not involve parameters depending 
on the previous strain-history. The dependence of the shape and size 
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of the yield surface on strain-history is thus taken to be expressible 
through the single parameter c. Infinitesimal changes of stress for 
which df=0 maintain the element on the point of yielding, and cause. 
only elastic changes of strain in a hardening material. The sign of f is 
conventionally chosen so that df>0 for an increment of stress producing 
further plastic strain, while if df<0 the element unloads elastically. 

The elastic part of the strain is neglected henceforward. The relations 
between stress and strain are taken to be 


4 9 ae. 
dey=ha afi df>0. . . 2 ss Q) 


g (the plastic potential) and h are scalar functions of the reduced stress oj; ; 
parameters in g depending on strain-history are held constant in the 
differentiation. When the element is isotropic f, g and h must be functions. 
only of the two independent invariants of the tensor oj; Since g does 
not involve the hydrostatic part of the stress, (2) implies that there is no 
plastic volume change. 

It is normally assumed, in the absence of definite experimental evidence, 
. that the yield function f and the plastic potential g are the same*. With 
this simplification equation (2) becomes ~ 


af 
dey—hae- af; df>0. . . (3) 


Certain extremum principles (and an associated uniqueness theorem) may 
then be proved (Hill 1950). The proof of one of them will be repeated 
here for the sake of a later analogy. 


(i) Maximum Work Principle. 

Equation (2) may be interpreted as asserting that the strain increment 
(represented in stress hyperspace) is parallel to the corresponding outward 
normal to the yield surface. In order that a unique state of reduced 
plastic stress should correspond to a prescribed increment of strain it is 
evidently necessary for this surface to enclose, and be concave at all 
points to, the origin. 

There will be two generators on the surface where the normals are 
parallel to a given direction ; uniqueness is secured by the condition 
that the plastic work is positive. 

Suppose, now, that an element is in a given state of hardening, defined 
by a certain yield-surface. Consider a plastic stress state o,,; and 
corresponding strain-increment de,;,. If oj; is a stress state lying within 
or upon the yield surface, the vector o,;—0;; is inclined at an acute angle 
to the outward normal at o;; (the yield surface being cylindrical and 
concave to the origin). Hence (using the summation convention) 


(o,,;—o;,)de;;>9, de,; £0, eee aaa hte (ay 


* More generally, g is taken equal to any convenient multiple or power of f. 
The resulting stress-strain relations are always the same if h is suitably altered. 
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the equality holding only when o,; and o;, differ by a hydrostatic 
component. This has the interpretation that the work done by the 
actual stress in a prescribed strain-increment is greater than that done 
by any other stress not violating the yield condition. The result does 
not involve any assumption about the dependence of the yield surface 
on the strain-history. 

Consider, next, a finite mass undergoing plastic deformation at 
all points. Let o;; be the actual stress, satisfying the yield condition and 
the equilibrium equations 


oy 0... en 


at all points. The associated strain-increment must be derived from a 
continuous displacement increment du,, such that 


lfa 7) 
If oj, is any stress state satisfying (5) and lying within or on the yield 
surface for each point of the mass, we have the virtual work equation 


S{(o3—o},)de,;} dV = f{(F ,—F; )du,} ds, 


where the first integral is taken through the volume, and the second 
over the surface of the mass, F,; and F; being the ith components of the 
respective outward forces acting on unit surface area. Hence, from (4) 


f{(F HF) dujfdSS0, se eee) 


the equality holding only when the stress states differ by a uniform 
hydrostatic stress. Equation (6) is true whether or not the states of 
hardening and anisotropy are non-uniform, and no assumption is made 
as to how they depend on strain-history. 


(ii) Converse of Maximum Work Principle. 


Suppose it to be given that, for any prescribed strain-increment, the 
corresponding state of plastic stress in an element is such that the work 
done is stationary with respect tc infinitesimally near stress states, lying 
on the yield surface. Then it is evident that the normal to the yield 
surface at the corresponding point must be parallel to the strain-increment 
vector. A plastic potential therefore exists and it is identical with the 
yield function. The latter, moreover, must be cylindrical with generators 
parallel to the direction 5,;,, since any possible strain-increment has zero 
hydrostatic part (there being no plastic change of volume). If the work 
done is an absolute maximum, with respect to all stress states within or 
on the yield surface, the surface must be concave to the origin at every 
point. 

This is the basis of the method adopted later to derive equation (3) 
from the properties of single crystals. 
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(ili) Complementary Minimum Principle. 

Consider a mass in an equilibrium plastic state (o;;, de,;) and let de;, 
be any strain-increment; de;; and de;,; are derived from continuous 
displacement increments ae oad du; taking the same values on the 


surface and having zero eaten everywhere. It may be proved 
(Hill, op. cit., p. 67) that 


S (| ol] des]) dV<J]( (o;;|| de;,| ) d say Dlsp a rgeane (el) 
where |o:;[ =~+/(e;,04;), ete., 


y-y 


provided the material is isotropic and the i Rae is a circular 
cylinder in principal-stress space ; i.e. f=oj; o/, in (3). If the state of 
hardening is uniform, (7) reduces to Markov’s att 


S]de,,|dV< f|dej;| UN Sie Weed Mele ae ick ial 3) 
The analogue of (7) for any ee Spa is 


where (4) is used to derive the See eh of foe the ea stress-state 
corresponding to de;,. 


§ 3. PLAsTIcITY OF A SINGLE CRYSTAL. 


We assume that the only mechanism of plastic distortion in a single 
crystal is by glide parallel to preferred planes and directions, and we 
examine the consequences of such an assumption. Let there be n possible 
slip directions a,,...a,,. (on associated glide planes), for a certain crystal 
lattice. The displacement of any point of the crystal, due to 
simultaneous infinitesimal shears dy,, . . . dy, (positive or negative), 
in the respective directions may be found by simple geometry. The 
corresponding tensor representing the homogeneous strain has five 
independent components (its hydrostatic part being necessarily zero), 


each of which is found as a certain linear combination of dy,, . . . dy. 
For given dy, (k=1, . .. n) the strain tensor is thereby uniquely 
determined. 


' On the other hand, when the components of the strain tensor are 
given (the hydrostatic part being zero), we have » unknown dy, and 
5 equations between them. If 7 is less than 5 (as in a hexagonal metal) 
a combination of shears cannot be found to produce an arbitrary strain. 
If n is equal to 5 there is a unique set of the dy, if the determinant of 
the coefficients is non-zero. If n exceeds 5, a set of 5 slip directions can 
be selected in any one of "C, ways, but the corresponding set of shears 
can only be found when the associated’ determinant is non-zero, 
i.e. when the slip directions form an independent set. It may also be 
possible to find combinations of 6 or more shears to produce the 
prescribed strain. 
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In a fully annealed crystal the critical shear stress appears not to 
depend on which slip-direction is activated, nor on the sense of the shear. 
During continued deformation active and non-active directions harden 
equally, to a first approximation (Taylor and Elam 1925). Small differences 
are observed, however; for example, double slipping does not always 
begin at the stress which would be expected if the previously non-active 
direction hardened at the same rate as the active one. Differences 
between the critical shear stresses may also be due to microscopic internal 
stresses set up during cold-working. These stresses can be expected to 
favour slip in certain directions and senses, and they will therefore be ~ 
effective whenever the load is varied in orientation or sense. At present 
there is no clear evidence as to the importance of such effects in crystals 
of pure metals, but a Bauschinger effect has been observed in cold-worked 
brass, from which it could be largely removed by a mild annealing not 
reducing the average hardening (Sachs and Shoji 1927). 

In the subsequent analysis we allow for these possibilities by 
associating two critical shear stresses with each of the preferred 
directions : 7, for slip in the sense a,, and 7; for slip in the sense —a,. 
It is assumed that the hydrostatic part of the applied stress is without 
influence on these critical values ; this has been found true in many pure 
metals (as distinct from alloys) provided the hydrostatic stress is not more 
than a few times the shear yield stress (e.g. Polanyi and Schmid 1923, 
Taylor 1927). 


(i) Maximum Work Principle for a Single Crystal. 
Let de,;; be a prescribed increment of strain in a single crystal, and 
let o,;; be a stress (not necessarily unique) which will produce this strain. 


Let o;, be any other stress which does not violate the yield conditions. 
Then 


dW —dW*=(0;;—0;;) de,,=2(7t—1*) dy, 
where the summation is over all shears of a set equivalent to de;; and 
activated by o,;;. We take each dy to be the shear (positive or negative) 
in the sense a of the corresponding slip-direction, so that 7 and r* are 


then the shear stresses in the same sense. Then, when dy>0 in the kth 
direction, 


t—T*=7,—T*>0 
and when dy <0 in the kth direction, 
T—T*=—7T1—T*<0 
since, by hypothesis Trt <7;: 
All products (r—r*)dy are therefore positive or zero, and so 
(G50) Aes ee0t amity cuaot oot swe ato) 


The equality holds only when 7 and +* are equal in all active directions, 
and, if there are at least 5 of these, o,; and o}, will generally differ only 
by a hydrostatic stress. Particular strains may be obtainable by applying 
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any one of a number of stresses (not differing merely by a hydrostatic 
component), but these stresses all do the same work. For if o;, is also a 
physically possible stress, 


in addition to (9), and therefore oj,de;;=0;; de,;. 

It is also to be remarked that many sets of shears (even associated 
with the same stress) may produce, in certain cases, the same strain. 
That is, more than one set of shears may be physically, as well as 
geometrically, possible. However, since there are only 5 disposable 
stress components it is clear that only for certain sets of 5 slip-directions 
is it possible to find a stress to activate them in prescribed senses, since 
the critical shear stress must also not be exceeded in the non-active 
directions. 


(ii) Minimum Shear Principle for a Single Crystal. 

Let de;; be a prescribed increment of strain and let o,, be a stress which 
will produce this strain by activating a set of shears dy. Let dy* be any 
set of shears which are geometrically equivalent to the prescribed 
strain, but which are not necessarily operable by any stress satisfying 
the yield conditions. The virtual work equation can be written in 


two ways : 
04; dé y= Lit dy=2't dy*. 


If +,=7; for each direction (no Bauschinger effect), 
a7 dy=27,|[dy|, and 27dy*<2r,|dy* |. 
Hence 27, | dy |< 27; | dy* |. oO Bee) 
In particular, when all directions are equally hardened, 
ideale WA | jhe Mgmt ogy sdie stein dbl) 


This states that the sum of the absolute values of a physically possible 
set of shears producing a given strain is less than that for a set which is 
only geometrically possible. If more than one set is physically possible 
the sums of their absolute values are equal. 

We have thereby proved the hypothesis (11) suggested by Taylor (1938), 
and extended it (10) to the case when the hardening is unequal. 


§ 4, POLYCRYSTALLINE AGGREGATE. 


The experimental laws of plastic deformation in an aggregate (as a 
whole) express relations between macroscopic measures of stress and 
strain. An ideal experiment is designed so that the measuring device 
gives an average value of stress or strain over a large number of crystal 
grains. In the interpretation of such measurements the following two 


assumptions are involved. 
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(a) The measuring device extends over a volume of such size that the 
distribution of orientations and hardnesses of crystals (or perhaps 
fragmented “crystallites ’) does not vary significantly from one part 
to another. That is, the specimen is regarded as homogeneous in the 
macroscopic sense. This does not preclude a macroscopic state of 
anisotropy, since the actual orientations are not necessarily randomly 
distributed over all possible ones. 

It is convenient to refer to the smallest cubical volume possessing this 
property as a “‘ unit ” cube. 

(b) No correlation exists between microscopic stress and position 
over any plane section of “unit” area. This is the necessary and 
sufficient condition for the stress resultant over such a unit section to be 
a single force through the centroid of the section. The Cartesian 
components of the resultant force per unit area are, of course, the 
conventional measures of the corresponding components of macroscopic 
stress. ‘These form a tensor (as a direct consequence of the conditions 
for linear equilibrium) and the tensor is symmetric if the stress resultant 
on a plane is a force, and not a force and irreducible couple. If a 
correlation did exist between microscopic stress and position (and this 
appears possible in principle), the macroscopic stress tensor would not 
necessarily be symmetric, and the conditions for angular equilibrium 
would take another form. These remarks are equally pertinent, of course, 
to the concept of microscopic stress in a single crystal and its dependence 
on the distribution of interatomic forces. 

The definition of a normal component of an infinitesimal homogeneous 
strain is the increment of distance between two points initially unit 
distance apart in the direction under consideration. The definition 
of a component of shearing strain is based on the total relative tangential 
displacements of two pairs of perpendicular planes at unit distance 
apart. The analogous definition of the macroscopic strain-increment 
components for a “ unit cube ” of an aggregate is 


dB j;=4 J (l;dujtldu,d8, . . . . . . (12) 


where the integral is taken over the surface of the unit cube, 1; being 
the unit outward normal, and du; the incremental microscopic 
displacement (the coordinate axes are taken parallel to the cube edges). 
When the microscopic displacement is a continuous function of position 


(no opening of cavities or sliding of one grain over another), Green’s 
theorem gives 


irfa a 
dE y= 5 | {an (du;)-+ = (an) dV, 


“J 
3 | (de,;) dV, 


taken through the volume of the unit cube, where de,, is the microscopic 
strain-increment. Since the latter is a symmetric tensor with zero 
hydrostatic part so also is dE;;. 
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Consider the work of deformation within a unit cube. It is 


provided the displacement function is continuous and the equilibrium 
equations are satisfied. At grain boundaries, or elsewhere, the 
microscopic stress is not necessarily continuous in all its components, 
and only the stress resultants acting on opposite sides of such a TS 
must be the same for equilibrium. 

We now postulate that the grains are distributed in the unit cube in 
such a way that there is no correlation between any component of. the 
microscopic stress and any component of the displacement over any 
plane section of unit area. This is so provided that 


foydAX fdu,dA=f(o,;du,)dA(t, j, k=1, 2,3), . . (14) 


where the integration extends over the unit section. In particular, this 
is satisfied when the stress is uniform or when the strain is uniform (the 
stress being distributed according to (b) above). Applying (14) to each 
face cf the unit cube, we find from (13) that 


dW= (oe OV=S,0H,,, . . > » ~ (15) 
where S,,; is the macroscopic stress tensor defined in (6). 


(i) Maximum Work Principle for an Aggregate. 

Let S;; be any sia croscopic stress corresponding to an equilibrium 
microscopic distribution o;, not violating the yield conditions at any 
point of a unit cube. Then, from (15), 


(S83) dye tlc; o%) deydV20, a. ¢! (18) 


since (9) holds at every point. We have thereby derived (4) from assumed 
properties of a single crystal. It is worth noting that (16) is still true 
even if du, is tangentially discontinuous across certain surfaces provided 
the corresponding frictional stress vanishes. This is approximately true 
in high-temperature creep where relative slip occurs between grains. 


(ii) Minimum Shear Principle for an Aggregate. 

Let du; and du* be two continuous displacement distributions, with 
ZeTO divergence, taking the same values on the surface of a unit cube. 
du, is associated with an equilibrium distribution of stress o;; satisfying 
the yield conditions. Then 


f (04; de;;) V=J (o,,de%) dV or 8, dEy=8; ,dk;;, 
and so f(2r a )dV= (27 dy*) dv. 
If 7,7; for each direction at any point (no Bauschinger effect), 
S17 dy=Zr1,|dy|, 2X7 dy*<27;,|dy* |. 
Hence 8, GE y= J (27, | dy |) dV <J (27% | dy*)) ain. eens 1) 
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When the critical shear stress is uniform throughout the aggregate : 
f(2[dy|)dV<f(Z[ay* dV.  . . ~~ « (18) 


This is true whether or not the aggregate has the properties ascribed 

to the unit cube. There is a close formal analogy between (17) and (7), 

and between (18) and (8). . 
On combining (16) and (17) we have for a unit cube, when 7,=7% 


Another inequality that will be used later is 
Si; dHi;<J (27, | dy* |) dV, . . . . . . (20) 


where du; is any continuous distribution. 

It is to be noticed carefully that in establishing these extremum 
principles for an aggregate there is an implicit assumption, not merely 
that slip is the only mechanism of distortion, but that it is also a sufficient 
mechanism. In other words, it is assumed that the equations of 
equilibrium (5) can be satisfied throughout an aggregate by a stress 
distribution which, at the same time, will everywhere operate sufficient 
shears to produce a continuous displacement. For this to be 
mathematically possible, it may well be that the microscopic stress and 
strain must be allowed to vary continuously, and not be restricted to take 
constant values within each of a finite number of regions (which may be 
grains or parts of grains). If a continuous variation is envisaged, the 
relations between microscopic stress and strain, established for a crystal 
of a certain finite size, must be assumed to apply at a point. This, 
of course, is a procedure that is adopted in every branch of mechanics. 


§5. CALCULATION OF THE YIELD SURFACE. 

In § 4 (i) a maximum work principle for an aggregate has been shown 
to hold when slipping under a critical shear stress is the only microscopic 
mechanism of distortion. It follows from § 2 (ii) that a plastic potential 
governs the relation between macroscopic stress and strain-increment, 
and that the plastic potential is identical with the yield function. If, 
therefore, the yield function can be calculated for any particular metal 
the relationships between the ratios of the stress components and the 
ratios of the strain-increment components follow immediately from (3). 

We consider first what general properties of the yield function follow 
from the assumption of deformation by slip. It is evident that the 
hydrostatic part of the applied stress has no influence on yielding ; that is, 
the yield surface is cylindrical (though not, of course, necessarily circular). 
This follows since the superposition of a uniform hydrostatic stress 
throughout a plastic aggregate does not disturb the equilibrium of the 
microscopic stresses, nor does it alter the active slip-directions at any 
point, when the influence of elasticity is disregarded. In reality, however, 
a certain effect of hydrostatic stress would be expected in an aggregate 
of elastically anisotropic grains. 
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If the critical shear stress does not depend on the sense of slip, and if 

the aggregate is free from.internal stress when the applied loads are 
‘removed after cold-work, the yield surface for that state is symmetrical 
about the origin. That is, if S,, produces yielding when the aggregate is 
reloaded under constant stress-ratios so also will —S,,. For the elastic 
compatibility equations are linear in the stresses, and hence if a 
microscopic distribution o;, corresponding to §;; can be established by 
monotonic loading from a stress-free state so also can a distribution —o,, 
corresponding to —S;;._ Furthermore, if o,,; satisfies the yield conditions 
at any point so does —o;,, and if de,, is the increment of strain associated 
with o,;, —de,; is the increment associated with —c,,._ Plastic deformation 
would therefore be initiated under the stress —S,,. 

When the aggregate is macroscopically isotropic the yield surface can 
be shown to possess a six-fold symmetry in principal-stress space (Hill 
1950, p. 18). Since this symmetry does not depend on any particular 
mechanism of plastic distortion, we do not need to discuss it here. 

The exact calculation of the yield surface appears to be a matter of 
some difficulty. In the present paper we shall merely show that it is 
possible to calculate, with comparative ease, two cylindrical surfaces 
between which the yield surface must certainly lie when there is no 
Bauschinger effect. For this purpose we use the inequalities established 
in the last section. In one case we consider an aggregate in which the 
stress is uniform but displacement continuity is violated, and in the other 
case we consider an aggregate in which the strain is uniform but 
equilibiium is violated. 

Take first the uniform stress distribution 

ONT, Ty=Sil(Si:5 
where, for given stress-ratios r,,, A* is to be chosen so that the critical 
shear stress corresponding to o;, is just attained in one slip-direction 
at the “ weakest” point of the aggregate. It is then certain that oj; 
does not violate the yield conditions anywhere in the aggregate. oj; is 
therefore an equilibrium distribution such as is envisaged in (16) (the 
displacement function being continuous) which leads to 


Spon) oy | Samet tere ee beter 9o(22) 
Hence, for each direction r;; in stress hyperspace, (22) supplies a lower 
limit to the length of the corresponding “ radius ” to the yield surface. 
Since we know that the surface is cylindrical we need only take stresses 8,; 
with zero hydrostatic part in (21). 
Next take any wniform strain distribution de;,=dE;; with zero 
hydrostatic part, and such that r,; dK;,>0. Let dy* be any of the sets 
of shears which are equivalent to the strain dE; dy* is a function of 
position since the lattice orientation varies through the aggregate. Let a 
quantity .* be defined such that _ 


ur (ae | dy*|) dV |r ;; dk;,, ° e ° e e {23) 


vies (21) 
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where r;; is defined in (21) and the integration extends through the unit 
cube. Then, by (20), 
| Si; | ibs ae . . ° ° . . ° (24) 


According to (10) the least possible value of 4* for any choice of dE;; is 
to be obtained by taking dy* to be a physically possible set of shears for 
each orientation, the ‘“‘crystallites’’ being supposed completely 
independent. Then 

= (dH fou dV)/(ry dy), « - «ay » 29 (25) 
where o;, is a stress ae ae produce the strain dH;, in a free crystal. 
An advantageous choice of dE}, appears to be dE;;=r;; , (an arbitrary scale 


factor can obviously be omitted for convenience). Then, since 7,7 ,;=1, 


pt=rj,fo,dv. oo eke eM eta ea eete (26) 


This value of * is equal to the work which would be done if a unit cube 
of aggregate were split into its constituent regions of uniform orientation 
and each were separately given the strain 7;;. 


§6. DEPENDENCE OF YIELD SURFACE ON STRAIN-HIsTORY. 


We have not yet considered how the shape and size of the yield surface 
depend on the previous strain-history. For simplicity let us disregard 
possible Bauschinger effects both in a crystal grain and in the aggregate 
as a whole, and suppose also that the critical shear stress is the same in 
all slip directions. If there were no work-hardening it is evident that the 
shape and size of the yield surface would be preserved during continued 
deformation, provided the distribution of orientations did not vary. (This 
last condition will not be satisfied whenever a preferred orientation is 
developed during cold-working). When work-hardening occurs the 
grains harden by different amounts, due partly to the non-uniform 
distortion and partly to their different orientations. However, to a first 
approximation, it seems reasonable to expect that the yield surface will 
simply increase in size without changing shape. The size will be 
proportional to some mean value of the critical shear stress. 

There is some evidence (Taylor 1927) that the shear yield stress of 
a crystal is, very roughly, a function only of the total shear. That is, 


T=F(2y), 
where the function F is the hardening curve in simple slip. Now we have 
seen that the work done on the aggregate unit cube during an increment 
of strain is 
dW = [ (72 |dy|) dV 
=72' |dy|, 


say, where 7 and y now denote mean values over the aggregate. Hence, 
approximately 


aW=7 dz/F" 
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where F” is the slope of the shear-hardening curve. Since F’ is a function 
of + only, so also is W. Conversely, the mean critical shear stress is a 
function only of the total plastic work. 


These qualitative considerations lead, therefore, to the conclusion that 
the size of the yield surface is mainly a function of the total plastic work, 
and that other factors are secondary (in the absence of Bauschinger 


effects). This, indeed, is observed for face-centred cubic metals (see 
Hill 1950, pp. 27-32). 
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1. PROBLEMS SOLVED. 


Tue solution to this problem is important in estimating the accuracy 
of the capacitance method of measuring the irregularity of cotton yarn 
(Boyd 1949 and Walker 1950). The dielectric properties of cotton fibres 
are such that (Balls 1946) a singles cotton yarn may be represented by 
a uniformly twisted cylinder which, in the untwisted state, has one 
dielectric constant value along its length but a different, though axially 
symmetric, transverse value. 

We first find the field when the cylinder is placed in a general type of 
field F(R) [cos ndi—sin nj] (see equation (1)). This calculation is carried 
out by using a molecular model and mathematical process akin to that of 
Lorentz’s molecular dielectric theory, thus obtaining a pair of simple 
simultaneous integral equations (see (3)). If the twist is not too great 
a series solution is obtained when F(R)=R”. With this series the field when 
the cylinder is placed parallel to and between two infinite conducting 
planes at different potentials is found by means of multipliers which can be 
determined to any required accuracy. By assuming that finite plates do 
not alter the field, the charge on the plates is calculated and finally the 
capacity. For non-midway positions the calculated capacity is not quite 
exact but is a good approximation if the plates are wide compared with their 
distance apart. An approximate solution valid for most midway cases is 
given (see (10)). 

Acknowledgments are due to Mr. C. Illingworth of Manchester University 
for many helpful discussions and to Miss M. Castle of the Shirley Institute 
who has checked much of the mathematics and carried out the numerical 


work. 
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2. RESULTS. 


We consider an infinitely long circular cylinder of radius W which has a 
dielectric constant K, parallel to its axis but an axially symmetric trans- 
verse value K,. It is then twisted so that a line of molecules originally 
parallel to OZ, the axis, lies along the helix ~=R, cos bt, y=R sin bt, 
z=t, where b is a constant, O0<R<W and — o<t<-+ o (fig. 1). 

We then prove that an externally applied field 


F(R) [cos Nét—sin N@j] - . 2... (1) 


Fig. 1. 


Ni 


mnere. N==0:-1, 2°22... 3 , «=R cos d and y= sin 4, produces a distribution 
of molecular dipoles of density 


[A’ cos N¢é+ B’ cos (N+-2)¢] i+[—A’ sin N¢6+B’ sin (N+-2)d]7 
+C’sin (N+1)¢h . . . (2) 
where A’=k,A+bRJ, B’=k,B—bRJ, C’=k,C—23 
Amko(Ky+2)=3(K.—1), 4rk,(K,+2)=3(K,—1) 
2J(1-+b?R?)=(k,—k,) [DR(A—B)—C], 3C=4r0’ 
Ww 
A—F(R)=—27 | A(R)/3+B()—20N+1)RP {B's)ss as | 
R 
B=—27 | B(R)/3-+-4(B)—2(N pel) Bue de) { ‘ A/(S)S* a2 as | : 
(3) 
Outside the cylinder the field resulting from (2) is 
Ww 
D'[cos(N-++-2)di—sin(N+2)j]/R**?, D’=4(N-+-1)z E A’(S)SXt1d8 
(4) 
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Straightforward expansions solve the simultaneous integral equations 
(3) when F(R)=27(bR)X. Writing b|R=u, bW=U 


A’=u% Xa wm—Ax, say; Bae nem =Br, say; 


m 


Ces Ze, em=Cn, say; D’= 2n(N-+1)W5** Za, U2"/( Nek, 
=D x, say 
6K 14) =F a(Ky+2) [2Ky-+(N+1) (K,—1)/m], 3¢,=—2Iy,_-1(K +2), 
6K by) = —In-1(Ky +2) [2K,—(N+1) (K,—1)/(m+N+1)], 
Hi ove aes tp Ul s1; 
Jin|Sm-1= —1+(Ke—K,) (N+1)?/[4K ym(m-+-N-+1)], 
3(K,—1)Jp=2(K,+2)ao 
ao(K,+1)—(N-+1) (Ky—1) [b,0?+6,04/2+6,U%/3+ ...J=Ky—1... 


(5) 
Fig. 2. 
\, 
le 
K d a d > a= 


All quantities are first expressed in terms of a) then from the last 
equation the value of a) is found. 

When the cylinder is placed between two infinite parallel conducting 
planes kept at different potentials we assume the dipole density to be 


¥ pxRX[{At cos N$-+ By cos (N-+2)6i+{—A‘ sin Nd 
N=0 


+ By sin x (N+2)¢}j-+-Cx sin (N+1)¢k] Re cad Mt tacte aa Gd 
At each image of the cylinder we imagine a cylinder with the same 
distribution (6) except that, for odd N, the wy change sign. At the images 
of these images we assume (6) with unchanged sign. At the images of 
these, odd xy change sign and so on. We then calculate the field at the 


original cylinder due to the original field and these images and we expand 
it thus 


2 FyR* [cos N.di—sinN¢j] «7... (7) 
N=0 | 
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where the Fy are infinite linear combinations of the fy. From (5) the 
dipole density (7) would generate is calculated and must equal (6). By 
neglecting wy for N>M we obtain a finite number of linear equations for 
the wy. Numerical results indicate that by taking M large enough any 
desired accuracy can be achieved. 

The case of the cylinder midway between the plates is considered in 
detail (Fig. 2). Provided the potential difference between the plates is 
2d; we have, writing 


5 (1/s)"=S, 


1 


: oO ‘ oo 
Tg 3+ 2 Ham DemSam+a(2d) Ont, TDL = 22M 39 tom amSam+a(2d)—2"—4 
0 


Oh SOC aban Dim@2me(2d)e* Cl6. | Pps... 0... (8) 


The charge per unit length induced on the width AB (fig. 2) is 


L z (— if he? Limon wen (5) |e tanh B 
* m=0 (2m+1)! 2d dpe” B=nl,(2d). 


The charge induced on AB’ is the same with /, replaced by /,. In the 
case of finite plates we assume that the charge induced is not materially 
altered. Dividing by the potential difference we obtain the increase in 
capacity per unit length 6C. For fig. 2 a good approximation is 


C= 3Do {tanh [7l,/(2d)]+-tanh [7l,/(2d)]} 
- 4 (247d?—7?D5) ; 
A method of successive approximations is also given which may be 
applicable to other problems. 
The solution, in the form of an expansion, to an applied field 
R™ [cos Ndi—sin Ndj], MAN 


is also found (hence the more general case XfyR™ [cos Ndi—sin Njd] 
can be solved). 


(9) 


(10) 


REFERENCES. 
Batis 1946, Dielectric Properties of Raw Cotton, Nature, 1946, 158, 9-11. 
Boyp 1949, An Electronic Instrument for Measuring Weight Variations in 
Silvers, Rovings and Yarn, Journal of the Textile Institute (Transactions), 
Waker 1950, The Electronic Measurement of Sliver, Roving and Yarn 
Irregularity, Annual Conference of the Textile Institute, 1950. 


[ 432 ] 


XLVIII ‘CORRESPONDENCE. 


The Resistance-Minimum in Gold. 


By D. K. C. MacDonatp and I. M. TEMPLETON, 
The Clarendon Laboratory, Oxford*. 


[Received February 1, 1951.] 


Tue definite observation of a minimum in the electrical resistance of gold 
at low temperatures was first made by de Haas, de Boer and van den Berg 
(1933) and has so far received no satisfactory theoretical explanation. 
Later measurements by de Haas and van den Berg (1936, 1937) and 
Giauque, Stout, Clark (1937); Stout and Barieau (1939) have given 
evidence that the temperature of the minimum was a unique function 
of the relative resistance at the minimum (usually measured as (Ryyin/ Ro: ¢.,)) 
for annealed specimens and that for an ideally pure specimen the minimum 
would occur at absolute zero. It has been suggested (cf. Gorter 1938) 
that the resistance of any specimen would rise towards infinity if the 
temperature were lowered sufficiently far, but this suggestion hardly 
seems feasible particularly in view of the remarkable behaviour then to 
be expected at extremely low temperatures in a specimen of vanishingly 
small impurity. 

Two very tentative theories by one ‘of us (D. K. C. McD.) suggest that 
the resistance would pass through a point of inflexion below the minimum 
and then flatten off ultimately on approach to absolute zero at a higher 
value than the minimum. Both thebdries suggest, however, that the 
nature of the particular impurity atoms present would affect the situation | 
of the minimum. Since it appears likely that the major impurity present 
in the specimens available to the Dutch and American workers would 
generally be silver (the latter workers in fact worked specifically in one 
case with a range of Au—Ag alloys and report silver as the major impurity 
in the later work) it seemed profitable to investigate some prepared 
alloys containing other specific impurities. At the same time specimens 
with no intentionally inserted impurity were examined by us. 

Alloys containing respectively 0-05, per cent (by weight) of Cu and Ni 
and 0-2 per cent of Ni were made for us by Messrs. Johnson, Matthey 
and Co, Ltd., and we also experimented with a sample of fine gold which 
we drew down to a suitable diameter and a sample of “ spectrographically 
pure’ gold. As may be seen from the accompanying graph which also 
displays the data obtained elsewhere, the latter two specimens agree 
rather well with the curve of R,,;, against T,,;,, suggested by the other 
workers. (It should be noted that the minima derived in America are only 
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approximate since they do not have experimental results above ~4° K. 
or below ~10° K.). However the results for ‘our “ alloys” do not agree 
at all and this suggests very strongly that the presumption of a simple 
variation with residual resistance alone cannot be accepted. In fact, 
it seems tempting to assume that there is a characteristic curve for each 
type of impurity atom. 


Fig. 1. 
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[]: Data of : de Haas and v. d. Berg. A: Data of Giauque, Stout et al. 
©: Data of: Present work—‘ Specpure”’ and “ fine” gold. 
1. “ Hot” annealed. 


<>: Data of: Present work—0-05% Cu/Au < 2. “Cold” annealed. 
3. Unannealed. 


4, 0:05% ‘‘ Hot” annealed. 


F 5. 0:05294, Unannealed. 
He Data of: Present work—Ni/Au 6. one B Hotelannealed: 


7. 0:2% “Cold” annealed. 


Since it was thought that the precise method of annealing might be 
critical (e.g. Stout and Barieau annealed their wire for three hours at 
450° C.) we have annealed specimens at ~400° C. for 2-3 hours and also 
at ~850° C. (onset of bright red heat) for about the he period. There 
appears to be little difference in the cases of the “ pure” gold specimens 
and Cu/Au alloy but in the 0-2 per cent Ni/Au alloy a considerable 


2G3 
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improvement in residual resistance occurred with the “ hot” annealing 
over the “cold” process with a slight shift of the minimum to a lower 
temperature. We may also mention that in the 0-05 per cent Ni and Cu 
alloys annealing by either method produced a considerable minimum 
shift to a lower temperature—the points for the unannealed specimens 
are shown by dotted characters on the graph. 

Further experiments are in progress and it is hoped to publish the work 
in more detail later. 
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The Range-Energy Relation for Slow Alpha-Particles in Avr. 


By Joan M. Freeman and W. E. BurcHam, 
Cavendish Laboratory, Cambridge. 


[Received February 13, 1951.] 


DuRING some recent work on (pa) reactions (Burcham and Freeman 
1949 a, b) the ranges in air of a number of groups of alpha-particles with 
energies between 1 and 2-5 MeV. were compared. The measured range 
differences were used to obtain the energy differences for the corresponding 
groups of alpha-particles, with the aid of the range-energy relation given 
by Holloway and Livingston (1938). This range-energy relation was 
constructed on the assumption that the number of ions formed by an 
alpha-particle stopping in air was strictly proportional to its energy ; 
this assumption has recently been shown to be incorrect (Jesse and 
Sadauskis 1950), and a new range-energy relation has been given by 
these authors and also by Bethe (1950). The difference between the 
new and old range-energy curves is most significant for energies below 
about 2°5MeV., and the new relation is in good agreement with 
observations on the ®Li(n«) and B(n«) reactions, for which the alpha- 
particle ranges have been measured, and the energy releases are also 
accurately known as a result of magnetic deflection experiments (see 
Bethe 1950). The number of points on the range-energy curve for 
energies less than 2-5 MeV., for which direct measurements of range 
and energy have been made is however small. It was therefore thought 


* Or see: v. D. Bera, 1938, Thesis, Leiden, ‘‘De Electrische Weerstand van 
zuivere Metalen bij lage ...'Temperaturen ”’. 
+ Communicated by the Authors. 


Correspondence 435 


to be of some interest to use the range measurements made during the 
work on (pa) reactions to check the slope of the range-energy curve in 
the region from 2-5 down to 1-0 MeV. 

The measurements were made by a method previously described 
(Burcham and Freeman 1949 a, § 3 (a)). The reactions studied, together 
with the energy releases deduced from recent magnetic deflection 
measurements are shown below (Table I.). 


Reaction 


*Be(po)*Li 


10B (px) 7Be 


TABLE I. 


Energy Releases in (px) reactions. 


Q (MeV.) 


Reference 


2-121+0-012 
2-142 + 0-006 


Tollestrup, Fowler and Lauritsen (1949) 
Strait et al., quoted by Hornyak, Lauritsen, 
Morrison and Fowler (1950) 


1-148 + 0-006 
1-152 + 0-004 


1-147+0-010 


Brown, Chao, Fowler and Lauritsen (1950) 

Van Patter, Sperduto, Strait and Buechner 
(1950) 

Burcham and Freeman (1950) 


19F'( pa, )8O* 


1-977 + 0-008 
1-979 + 0-009 


Chao, Tollestrup, Fowler and Lauritsen 
(1950) 
Freeman (1950) 


19F'( pa.) 6O* 


1-204 + 0-008 


Chao, Tollestrup, Fowler and Lauritsen 


(1950) 


1-200+0-010 | Freeman (1950) 


Chao, Tollestrup, Fowler and Lauritsen 
(1950) 
Freeman (1950) 


19F(pas)!6O* | 1-002+0-008 


0-996 + 0-010 


At the time at which the measurements were made the Q-value for 
®Be(pa) seemed the best known, and all range comparisons were made 
with the alpha-particles from this reaction taken as standard. We have 
now chosen the 1!F(px,) reaction with E,—873-5 keV. and angle of 
observation 83°, as the standard, giving alpha-particles of energy 
2-272+0-010 MeV.; the corresponding mean range, from the curve 
given by Bethe (1950), is 1:157cem. The observed range differences 
then lead, after small corrections for differences of straggling and target 
thickness, to the mean ranges given in Table II. The energies shown in 
this Table were calculated, for the given conditions of bombardment, 
from the Q-values listed in Table I.; for the !F(p«) reactions the Q-values 
given by Chao, Tollestrup, Fowler and Lauritsen were used, since the 
values obtained by Freeman, although in good agreement, were based on 
a spectrometer calibration which involved the Q-value of the *Be(pz«) 
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reaction. This was assumed at the time to be 2:121 MeV., but the more 
recent measurement gives an appreciably higher value (see Table L.). 
The weighted mean of the two Q-values was used to calculate the energy 
given in Table II. for the, alpha-particles from the *Be(pa) reaction. 


TaBue II. 
Range Energy Relation. 


E,, (MeV.) Mean range 
Reaction P (Angle of em. of air at 15°C., 
hi observation 83°) 760 mm. 


19F(pa,)6O* 


2-:113+0-010 


1-068 + 0-008 


*Be(pa)*Li 


1-774+4 0-015 


0-895 +.0-008 


19F( po) 6O* 


1-649+ 0-010 


0-812 + 0-008 


19F(pag)18O* 


1-326 +0-010 


0-651 + 0-008 


0-559 + 0-011 


LOB(pa)7Li 1-113+0-007 


Fig. 1. 


ALPHA-PARTICLE ENERGY- MEV. 


Bethe 


Jesse & Sadauskis 


0s 0:7 0-9 4-1 1-3 


RANGE = CM.OF AIR AT 15°C AND 760 MM. 


These results are plotted in fig. 1, together with the curves proposed 
by Bethe (1950) and Jesse and Sadauskis (1950). The standard point 
for *F(pa,) at E,,=873-5 keV, has been placed on the former curve, and 
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all the other points, except that for *Be(pa), lie within 10 keV. of the 
curve. The position of the point for *Be(pa) suggests that the higher 
Q-value of the two given in Table I. is the more nearly correct. 

The results confirm the slope of the new range-energy curves over the 
energy range 1-0-2-3 MeV.; they are in slightly better agreement with 
the curve given by Bethe (1950) than with the figures given by Jesse and 
Sadauskis in the region where these differ. 
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XLIX. Notices of New Books and Periodicals received. 


Weltsystem, Weltather, und die Relativitaetstheorie. By K. JeEvLINcK. Pp. 450. 
(Verlag, Basel: Wepf and Co., 1949.) Price 45 8.Fr. 


Verstandlische Elemente der Wellenmechanik. By K. Jeuiinck. 1 Teil. Pp. 304. 
(Verlag, Basel: Wepf and Co., 1950.) Price 34 S.Fr. 


PROFESSOR JELLINCK may be known to many readers as the author of an 
extensive work on Physical Chemistry. In the preface to the book on the 
Theory of Relativity he confides that he has always wanted to understand this 
theory, but not finding a suitable textbook, he decided to write one both for 
himself and for the great number of experimental physicists and chemists. 

As the title suggests the author believes that the Theory of Relativity can be 
understood only by postulating an absolute frame of reference with respect to 
which all galaxies are, on the average, at rest (Weltbezugs system). The General 
Theory of Relativity states that even if it could be found, it must in no way be 
distinguishable from any other inertial frame of reference. It therefore seems 
unnecessary to introduce the Weltbezugs system in the first place. But this 
concept has, apparently, helped Professor Jellinck to a better understanding of 
the theory, and it may likewise help many others. Apart from this special 
attitude the book is easy to read and abounds in detailed discussions of particular 
problems. The tensor calculus is made familiar by working out many examples. 
This is probably the best way of becoming used to a new method. A discussion 
of the electro-magnetic equations, which should appeal strongly to experimenta- 
lists by showing at once the power and simplicity of the Special Theory of 
Relativity is unfortunately not included. 

The other book by Professor Jellinck is the first of two volumes on the 
principles and methods of wave mechanics. The demands on the reader with 
regard to mathematical knowledge are very small, as all necessary concepts are 
developed in the book. An attractive feature of the book is the discussion of 
various examples, such as harmonic oscillator, spinning top, and Kepler problem, 
from the points of view of Classical Mechanics, of the old Quantum Theory, and 
of Wave Mechanics. The final section is devoted to the problem of potential 
barriers and its applications. 

The text is amply illustrated and very readable. 

One unsatisfactory feature of the book must, however, be mentioned. It 
concerns the interpretation of some results of Wave Mechanics. For instance, 
in the case of the harmonic oscillator, Professor Jellinck asks : How is it possible 
for a particle with a definite energy to be found with a finite probability outside 
the classically allowed range? The answer is, of course, that as long as the 
energy is known exactly the position of the particle is completely unknown, but 
as soon as the position of the particle has been determined by experiment its 
energy is no longer known precisely. ; 

Professor Jellinck, however, argues as follows: The particle apparently can 
get to a place where its potential energy is greater than its total energy ; but its 
kinetic energy cannot be negative. Therefore the total energy as calculated by 
Wave Mechanics is only the constant average of the rapidly fluctuating actual 
energy of the particle. This fluctuating energy comes from the ether, which 
Professor Jellinck introduces very early in the book, together with the absolute 
frame of reference (Weltbezugs system). The ether has the property not only to 
propagate the wave function but also to interact energetically with the dynamical 
system. 

This way of explaining the “ paradoxes ’of Wave Mechanics is contrary to 
the logical scheme of Wave Mechanics, and misleading for the student. 
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Ultrasonics. By P. Vicournux. Pp. 163. (Chapman and Hall, 1950.) 
Price 25s. 


Tx author has dealt with the technique and with the theory of ultrasonics in 
simple terms which brings the work within the scope of a wide circle of readers 
including undergraduates. | 

Readers may remember the standard of work of Bergmann to which the author 
refers and which appeared in 1937. The author gives a good account of the 
progress which has been made since then in the study of attenuation, and in the 
effects associated with the vibrational states of molecules both in gases and in 
liquids. 


The Mathematical Theory of Plasticity. By R. Hm. [Pp. 354.] (Oxford: 
University Press.) Price 35s. 


THE mathematical theory of plasticity has undergone great development in the 
past dozen years or so, not only in its extension to more complex problems of 
industrial interest, but in its fundamental concepts and techniques. Much has 
been done to remove the reproach that it was one of the unrealistic branches 
of applied mathematics. This book is principally concerned with the behaviour 
of a model substance with the following properties: (i.) it is and remains isotropic; 
(ii.) it has a yield point; (iii.) this depends on the kind of stress in accordance with 
von Mises’ hypothesis (the yield surface in “ principal stress space ”’ is a right 
circular cylinder) or, as an occasional alternative, with Tresca’s hypothesis of 
critical shear stress ; (iv.) work hardening occurs as a function of the irrecover- 
able work done, or of Odquist’s integral plastic strain, or not at all; (v.) the 
principal axes of the plastic strain increment are parallel to the principal axes of 
the stress, frequently with the simplification that the principal strain incre- 
ments are proportional to the deviatoric principal stresses ; (vi.) the superimposed 
effects of elasticity are consistently regarded throughout. 

Considerable care is given to showing that these are the best assumptions, 
within the range of those capable of mathematical treatment, for conformity 
with the requirements of metal physics: and they lead to results in very 
satisfactory conformity with experimental measurements. ‘The weakest of the 
assumptions is clearly (i.), and a final chapter on plastic anisotropy develops the 
attack against this weakness. The author, who has been intimately concerned 
with the modern developments of the subject, has been in the fortunate position 
of being able to write what is simultaneously a complete text-book and a report 
of recent advancements. Ko Gole 


Mathematical Snapshots. By Prof. Huco Steryuaus, Oxford University Press. 
(Pp. 265, with 295 illustrations.] (London: Geoffrey Cumberlege.) Price 


27s. 6d. net. 


Tus book attempts with notable success to present mathematical ideas and 
theorems attractively to the layman with the help principally of illustrations. 
These are many and imaginative, while the text is always clear and to the point, 
if occasionally perhaps too brief. The author ranges widely, without evident 
design, from pursuit paths to polyhedra, from longimeters to loxodromes, from 
tesselations to topology, from counterfeit coins to cartography. Really difficult 
theorems are included, for example ones due to Minkowski and Lebesgue. 
A few results are proved ; most are stated with a challenge to the reader, but a 
bibliography of original papers is thoughtfully provided. Even a mathematical 


specialist, whatever his field, will surely find much in this book to divert ie 
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Electromagnetic Theory. By Ottver HeavistpE, [Pp. xxx+386.] (New York, 
Dover Publications Inc. British Agents: Interscience Publishers, Ltd.). 
Price 60s. 


Tuts is a photographic reproduction in one volume (12 in. by 9 in.) of the 
three volumes which appeared in 1883, 1899 and 1912. Four pages of the 
original appear on one page of the reproduction. Despite the reduction in 
size the print is quite as readable as in the reprint of 1925. Besides the original 
text there are a portrait of the author, a biographical sketch and an appraisal 
of Heaviside’s contribution to vector analysis, electromagnetism and communi- 
cation theory. J.D.E. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


